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NCNR  2008  Accomplishments  and  Opportunities 


Foreword 


I  am  very  pleased  to  present  this 
year's  summary  of  the  activities 
and  accomphshments  of  the 
NIST  Center  for  Neutron 
Research.  This  report  contains 
scientific  highUghts  selected 
from  a  few  of  the  many  diverse 
areas  of  research  that  illustrate 
the  large  body  of  exciting  work 
being  done  at  the  NCNR. 
However,  no  sampling  of 
highlights  can  fijlly  capture  the 
breadth  and  diversity  of  the  research  performed  each  year  at 
the  facility.  Therefore,  I  hope  that  this  report  serves  to  whet 
your  appetite,  and  I  invite  you  to  visit  our  website  at 
www.ncinr.nist.gov  to  learn  more. 

This  past  year  was  remarkably  busy  and  productive.  The 
reactor  operated  253  days  with  99  %  reliability  (operating 
on  the  days  originally  scheduled).  This  remarkable  record  of 
performance  enabled  another  banner  year  ot  research.  The 
NCNR  received  a  record  number  of  proposals,  and  more 
than  2200  scientists  participated  in  our  research  activities. 
The  resulting  scientific  output  of  320  published  articles, 
many  of  which  appeared  in  the  highest  impact  scientific 
journals,  is  outstanding. 


These  numbers  attest  to  the  importance  and  impact  of  neutron 
scattering.  However,  they  fail  to  capture  the  enthusiasm  and 
excitement  I  see  of  those  who  come  to  do  research  at  the  NCNR. 
For  me  these  attributes  were  evidenced  by  the  tremendous 
excitement  generated  by  work  on  the  recently  discovered  iron- 
based  superconductors,  seeing  the  excellent  presentations  at 
the  2008  American  Conference  on  Neutron  Scattering,  and 
observing  the  enthusiasm  ot  the  students  and  interns  who  come 
to  visit  and  work  at  the  NCNR.  I  hope  that  the  highlights  of  this 
annual  report  will  also  convey  to  you  the  spirit  of  excitement  and 
enthusiasm  of  those  who  performed  the  work. 

As  I  write  this  foreword,  we  are  wrapping  up  the  second  year 
of  the  NCNR  expansion  project  and  are  just  about  to  begin 
construction  of  the  new  guide  hall  extension  and  the  technical 
support  buildings.  Tremendous  progress  has  been  made  on  the 
cold  source,  the  new  guide  system,  the  instruments,  and  on  the 
control  room  upgrade  for  the  reactor.  This  expansion  represents 
an  enormous  opportunity  that,  when  completed,  will  be  a 
significant  enhancement  of  our  cold  neutron  capability. 

Thanks  to  the  hard  work  and  dedication  of  our  staff  and  the 
users  of  the  NCNR  facility,  we  continue  to  enjoy  success.  Their 
continued  effort  during  the  next  phase  of  the  NCNR  expansion 
in  the  coming  years  is  the  key  to  a  bright  future.  To  them  I  extend 
my  deepest  gratitude.  This  is  indeed  an  exciting  time  for  neutron 
science  in  the  United  States. 


The  NIST  Center  for  Neutron  Research 


Neutrons  provide  a  uniquely  effective  probe  of  the  structure 
and  dynamics  of  materials  ranging  from  water  adhering 
to  clefts  within  proteins  to  ordering  magnetic  moments  in 
giant  magnetocaloric  materials.  The  properties  of  neutrons 
(outhned  below)  can  be  exploited  using  a  variety  of  measurement 
techniques  to  provide  information  not  otherwise  available.  The 
positions  of  atomic  nuclei  in  crystals,  especially  of  those  of 
light  atoms,  can  be  determined.  Atomic  motion  can  be  directly 
measured  and  monitored.  Neutrons  are  especially  sensitive 
to  hydrogen,  so  that  moving  water  can  be  imaged  in  fuel  cell 
channels,  hydrogen  motion  can  be  followed  in  new  H-storage 
materials  and  water  take-up  can  be  monitored  as  cement  sets. 
Residual  stresses  such  as  those  deep  within  oil  pipeline  welds  or 
in  highway  trusses  can  be  mapped.  Neutron-based  measurements 
contribute  to  a  broad  spectrum  of  activities  including  in 
engineering,  materials  development,  polymer  dynamics,  chemical 
technology,  medicine,  and  physics. 


Cold  neutron 
Guide  Hall  and 
instruments. 
A  second 
Guide  Hall 
is  under 
construction. 


The  NCNR's  neutron  source  provides  the  intense,  conditioned 
beams  of  neutrons  required  for  these  types  of  measurements.  In 
addition  to  the  thermal  neutron  beams  from  the  heavy  water  or 


graphite  moderators,  the  NCNR  has  a  large  area  liquid  hydrogen 
moderator,  or  cold  source,  that  provides  long  wavelength  guided 
neutron  beams  for  the  major  cold  neutron  facility  in  the  U.S. 

There  are  currently  29  experiment  stations:  four  provide  high 
neutron  flux  positions  for  irradiation,  and  25  are  beam  facilities 
most  of  which  are  used  for  neutron  scattering  research.  The 
subsequent  pages  provide  a  schematic  description  of  our 
instruments.  More  complete  descriptions  can  be  found  at  www. 
ncnr.nist.gov/instruments/.  The  NCNR  is  currently  adding  a 
second  guide  hall  with  five  new  instruments  in  a  major  expansion 
(see  pp.  4-5). 

The  Center  supports  important  NIST  measurement  needs,  but 
is  also  operated  as  a  major  national  user  facility  with  merit- 
based  access  made  available  to  the  entire  U.S.  technological 
community.  Each  year,  more  than  2000  research  participants 
from  government,  industry,  and  academia  from  all  areas  of  the 
country  use  the  facility.  Beam  time  for  research  to  be  published 
in  the  open  literature  is  without  cost  to  the  user,  but  fuU  operating 
costs  are  recovered  for  proprietary  research.  Access  is  gained 
mainly  through  a  web-based,  peer-reviewed  proposal  system  with 
user  time  allotted  by  a  beamtime  allocation  committee  twice 
a  year.  For  details  see  www.ncnr.nist.gov/beamtime.html.The 
National  Science  Foundation  and  NIST  co-fund  the  Center  for 
High  Resolution  Neutron  Scattering  (CHRNS)  that  operates  six 
of  the  world's  most  advanced  instruments.  Time  on  CHRNS 
instruments  is  made  available  through  the  proposal  system.  Some 
access  to  beam  time  for  collaborative  measurements  with  the 
NIST  science  staff  can  also  be  arranged  on  other  instruments. 


Why  Neutrons? 

Neutrons  reveal  properties  not  readily  probed  by  photons  or  electrons.  They  are  electrically  neutral  and  therefore  easily  penetrate  ordinary 
matter.  They  behave  like  microscopic  magnets,  propagate  as  waves,  can  set  particles  into  motion  losing  or  gaining  energy  and  momentum  in 
the  process,  and  they  can  be  absorbed  with  subsequent  emission  of  characteristic  radiation. 

WAVELENGTHS  -  in  practice  range  from  =  0.01  nm  (thermal)  to  =1 .5  nm  (cold)  [1  nm  =  1 0  A),  allowing  the  formation  of  observable 
interference  patterns  when  scattered  from  structures  as  small  as  atoms  to  as  large  as  cells. 

ENERGIES  -  of  miliielectronvolts,  the  same  magnitude  as  atoms  moving  in  solids  or  liquids,  waves  in  magnetic  materials,  and  vibrations  in 
molecules.  Exchanges  of  energy  between  neutrons  and  samples  as  small  as  nanoelectronvolts  and  as  large  as  tenths  of  electronvolts  can  be 
detected. 

SELECTIVITY  -  in  scattering  power  varies  from  nucleus  to  nucleus  somewhat  randomly.  Specific  isotopes  can  stand  out  from  other  isotopes  of - 
the  same  kind  of  atom.  Specific  light  atoms,  difficult  to  observe  with  x-rays,  are  revealed  by  neutrons.  Hydrogen,  especially,  can  be  distinguished 
from  chemically  equivalent  deuterium,  allowing  a  variety  of  powerful  contrast  techniques. 

MAGNETISM  -  makes  the  neutron  sensitive  to  the  magnetic  spins  of  both  nuclei  and  electrons,  allowing  the  behavior  of  ordinary  and  exotic 
magnetic  materials  to  be  detailed  precisely. 

NEUTRALITY  -  of  the  uncharged  neutrons  allows  them  to  penetrate  deeply  without  destroying  samples,  and  pass  through  walls  conditioning  a 
sample's  environment,  permitting  measurements  under  extreme  conditions  of  temperature  and  pressure. 

CAPTURE  -  characteristic  radiation  emanating  from  specific  nuclei  capturing  incident  neutrons  can  be  used  to  identify  and  quantify  minute 
amounts  of  elemental  isotopes  in  samples  as  diverse  as  ancient  pottery  shards  and  lake  water  pollutants. 
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[1]  BT-5  Perfect  Crystal  Ultra-Small  Angle 

Neutron  Scattering  (USANS)  Diffractometer 
for  microstructure  up  to  10"^  nm,  part  of 
CHRNS. 

[2]  BT-4  Filter  Analyzer  Neutron  Spectrometer 
with  cooled  Be/Graphite  filter  analyzer  for 
chemical  spectroscopy. 

[3]  BT-2  Neutron  Imaging  Facility  for  imaging 
hydrogenous  matter  in  large  components 
such  as  water  in  fuel  cells  and  lubricants  in 
engines,  in  partnership  with  General  Motors. 

[4]  BT-1  Powder  Diffractometer  with  32  detec- 
tors; incident  wavelengths  of  0.208  nm, 
0.154  nm,  and  0.159  nm,  with  highest 
resolution  of  §d/d  =  8  x  lO.'"* 


[5]  BT-9  Thermal  Triple  Axis  Spectrometer  for 
measurements  of  excitations  and  structure 
in  condensed  matter. 

[6]  BT-8  Residual  Stress  Diffractometer 
optimized  for  depth  profiling  of  residual 
stresses  in  large  components. 

[7]  VT-5  Thermal  Neutron  Capture  Prompt 
Gamma-ray  Activation  Analysis  Instrument 
with  a  neutron  fluence  rate  of  3  x  10^/cm^/s 
used  for  quantitative  elemental  analysis  of 
bulk  materials.  Generally  used  for  the 
analysis  of  highly  hydrogenous  materials 
(=1  %  H)  such  as  foods,  oils,  and  biological 
materials. 

[8]  BT-7  Thermal  Triple  Axis  Spectrometer  with 
large  double  focusing  monochromator,  and 
interchangeable  analyzer/detectors 
systems. 

[11] 

[10] 

[12] 


D  KIC 


The  Center  for  High  Resolution  Neutron  Scattering  (CHRNS)  is  a  partner- 
ship between  NIST  and  the  National  Science  Foundation  that  develops 
and  operates  neutron  scattering  instrumentation  for  use  by  the  scientific 
community.  The  following  instruments  are  part  of  the  Center:  1  (USANS), 
13  (HFBS),  14  (NG-3  SANS),  15  (DCS),  22  (NSE),  and  23  (SPINS), 
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[9]  NG-0  MACS  cold  neutron  triple  axis  crystal 
spectrometer  (under  construction)  with 
double  focusing  monochromator  and 
multiple  crystal  analyzer/detectors  that  can 
be  flexibly  configured  for  several  energies 
simultaneously  or  for  high  throughput  at  one 
energy 

[10]  NG-1  Advanced  Neutron  Diffractometer  / 
Reflectometer  (AND/R)  a  vertical  sample 
reflectometer  with  polarization  analysis  and 
off-specular  reflection  capabilities  for  mea- 
suring reflectivities  down  to  10"^. 

[11]  NG-1  Vertical  Sample  Reflectometer  instru- 
ment with  polarization  analysis  capability  for 
measuring  reflectivities  down  to  10"^ to 
determine  subsurface  structure. 

[12]  NG-1  Cold  Neutron  Depth  Profiling  instru- 
ment for  quantitative  profiling  of  subsurface 
impurities. 


[14]  NG-3  30  m  SANS  for  microstructure  mea- 
surements, part  of  CHRNS. 

[1 5]  NG-4  Disk  Chopper  Time-of-Flight  Spec- 
trometer for  studies  of  diffusive  motions  and 
low  energy  dynamics  of  materials.  Wave- 
lengths from  =  0.18  nm  to  2.0  nm  give 
corresponding  energy  resolutions  from 
=  2  meVto  <  10  i.ieV,  part  of  CHRNS. 

[16-18]  NG-6  Neutron  Physics  Station  offering  three 
cold  neutron  beams  having  wavelengths  of 
0.5  nm,  0.9  nm,  and  "white"  that  are  avail- 
able for  fundamental  neutron  physics  experi- 
ments. 

[19]  NG-7  30  m  SANS  for  microstructure  mea- 
surements, in  partnership  with  NIST,  and 
ExxonMobil. 

[20]  NG-6  Neutron  Physics  Test  Bed  for  develop- 
ing measurement  techniques  and  detectors. 


[13]  NG-2  Backscattering  Spectrometer  (HFBS) 
high  intensity  inelastic  scattering  instrument 
with  energy  resolution  <  1  |ieV,  for  studies  of 
motion  in  molecular  and  biological  systems, 
part  of  CHRNS. 


NG-7  Prompt  Gamma  Activation  Analysis, 
cold  neutron  fluxes  allow  detection  limit  for 
[21]  H  of  1  jjg  to  10  jig,  depending  on  the  matrix. 
Focused  beams  are  available  for  profiling. 

NG-5  Neutron  Spin-Echo  Spectrometer 
(NSE)  for  measuring  dynamics  from  100  ns 
[22]  to  10  ps,  in  partnership  with  ExxonMobil, 
part  of  CHRNS. 

NG-5  Spin-Polarized  Triple  Axis  Spectrom- 
eter (SPINS)  using  cold  neutrons  with 
[23]  position  sensitive  detector  capability  for 
high-resolution  studies,  part  of  CHRNS. 

NG-7  Neutron  Interferometry  and  Optics 
Station  with  perfect  crystal  silicon  interfer- 
[24]  ometer.  A  vibration  isolation  system  provides 
exceptional  phase  stability  and  fringe  visibil- 
ity. 

NG-7  Horizontal  Sample  Reflectometer 
allows  reflectivity  measurements  of  free 
[25]  surfaces,  liquid/vapor  interfaces,  as  well  as 
polymer  coatings. 
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-a  multi-year  plan  to  meet  strong 
U.S.  demand  for  cold  neutron 
measurement  capability  by 
creating  new  beamlines  and 
instruments. 


New  guides  come  from  the  present  cold  source 
and  pass  through  confinement  wall  segments 


Control  Room  of  the  NCNR  Neutron  Source 
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MACS 


New  cold  source  chamber 
to  be  placed  at  the  MACS 
beamport 


4 


Present  Guide  Hall 
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2008  Activities 


Construction  contracts 
awarded,  ground  breaking 
scheduled 

Instrumentation  placement 
and  detailed  design 
underway 

Guide  design  completed, 
procurement  undenA/ay 

New  Cold  Source  design 
and  several  tests  completed, 
fabrication  and  procurement 
underway 

Control  Room  upgrade 
with  Brookhaven 
National  Laboratory 
carried  out,  bids  under 
evaluation 


New  cold  source  "thimble".  The  cold 
source  chamber  resides  within  the 
left  end  of  this  thimble,  and  cold 
neutrons  pass  out  to  the  right. 
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McMaster's  Pat  Clancy  and  Mehmet  Ramazanoglu 
prepare  a  sample  for  an  NG-3  SANS  measurement. 


NCNR's  William  Ratcliff  (left)  at  BT-9  with  Rutgers'  Valery 
Kiryukhin  and  Seongsu  Lee 


Helen  Hanson  (Brown  U.)  prepares  a  sample  during  the 
14th  IMCNR  Summer  School. 


James  Mclver  and  Bill  Clow  provide  entertainment 
at  the  holiday  party. 


NCNR's  Bulent  Akgun  [center]  at  the  NG-7  reflectometer 
with  U  Akron's  Gokge  Ugur  and  Wumin  Yu 


UCSD's  John  Goodkind  prepares  a  dilution  refrigerator  insert 
for  a  run  at  DCS 


Summer  School  students  focus  on  refiectometry. 
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NCNR's  David  Mildner  discussing  USANS  results  with  U.  Del.'s 
Sheng  Zhong  [left]  and  Congqi  Yan 


UCSD's  peripatetic  Sonny  Sinha  (2nd  from  left)  alights  for  an 
experiment  and  nearly  comes  to  rest  to  pose  with  NCNR's  John 
Copley,  Jack  Rush  and  Pat  Gallagher 


NCNR's  Dan  Neumann  and  Pat  Gallagher  flank 
Summer  Undergraduate  Research  Fellows  (SURF). 


Polymer-coated  magnetic  nanofibers  bring  together 
VTech's  Will  Miles  (left)  and  Rick  Davis  (middle  right] 
with  U  of  W.  Aus.'s  Matt  Carroll  (middle  left]  and  Rob 
Woodward  at  the  NG-7  SANS. 


NCNR's  Craig  Brown  guides  DCS  tutorial  participants. 


Tutorial  participants  Rolando  Valdez  (UMd)  and  Khaled 
Aamer  (NIST  Polymers  Division)  confer. 
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Magnetic  Order  versus  Superconductivity  in  the  Iron-based 
La(0^  /^)FeAs  Systems 

C.  de  la  Cruz  ^-2,  Q.  Huang^  J.  W.  Lynn^  J.  Li^^  W.  Ratcliff  IP,  J.L.  Zarestky^, 
H.A.  Mook2,  G.R  Chen^  J.L  Luo^  N.L  Wang^  and  P.  Dai^  ^ 


The  recent  discovery  of  superconductivity  in  the  rare- 
earth  (R)  iron-based  oxide  systems  [i^O^^F^FeAs  and 
(Ba,  Sr,  Ca)Fe,As  J  opens  a  new  field  because  these  materials 
are  the  first  oxide  superconductors  not  containing  copper 
with  transition  temperatures  (T)  exceeding  50  K.  The 
parent  (non-superconducting)  LaOFeAs  material  is  metallic 
but  shows  anomalies  near  150  K  in  resistivity,  dc  magnetic 
susceptibility,  and  specific  heat,  while  optical  conductivity 
and  theoretical  calculations  suggest  that  these  anomalies  are 
caused  by  the  development  of  a  spin-density-wave  magnetic 
order.  In  the  high-transition  temperature  (high-T)  copper 
oxides  it  is  generally  believed  that  antiferromagnetism  plays  a 
fi.mdamental  role  in  the  superconducting  mechanism  because 
superconductivity  occurs  when  mobile  'electrons'  or  'holes' 
are  doped  into  the  parent  antiferromagnetic  Mott  insulator 
compounds,  and  the  question  naturally  arises  whether  the 
superconductivity  in  these  two  classes  of  materials  have  a 
common  origin. 


The  present  neutron  measurements  demonstrate  that  LaOFeAs 
undergoes  an  abrupt  structural  transition  below  ~  150  K, 
distorting  the  tetragonal  crystal  symmetry,  and  this  is  the  origin 
of  the  transport  anomalies  in  these  systems.  However,  the 
structural  transition  is  closely  followed  by  the  development  of 
long  range  spin-density-wave-type  antiferromagnetic  order  at 
»  137  K.  The  ordered  moment  is  quite  small,  indicating  that  the 
magnetism  is  itinerant  in  nature,  but  with  a  simple  commensurate 
magnetic  structure.  Doping  the  system  with  fluorine  suppresses 
both  the  magnetic  order  and  structural  distortion  in  favor  of 
superconductivity.  Therefore,  much  like  high-  T  copper  oxides, 
the  superconducting  regime  in  these  Fe-based  materials  occurs  in 
close  proximity  to  a  long-range  ordered  antiferromagnetic  ground 
state. 

Neutron  diffraction  was  employed  to  study  samples  of 
polycrystalline  non-superconducting  LaOFeAs  and 
superconducting  La(0^  FJFeAs  with  x  =  0.08  {T  =  26  K)  [1]. 
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FIGURE  1 :  Temperature  dependence  of  the  [2,2,0]  structure  reflection.  Left:  The  diffraction  intensity  profile  in  the  inset  shows  that  the  angle 
dependence  of  the  (2,2,0].^  reflection  abruptly  splits  into  two  peaks  [(4,0,0)^  and  (0,4, 0)^,  {O  denotes  orthorhombic)]  as  a  function  of  temperature. 
The  graph  shows  the  data  at  175  K,  and  in  the  ground  state  at  4  K.  A  structural  transition  has  occurred  from  the  high  temperature  structure  with 
tetragonal  symmetry  to  orthorhombic  symmetry  [Cmma].  The  right  plot  shows  the  intensity  (red  symbols)  at  the  peak  position  for  the  tetragonal  peak, 
which  decreases  rapidly  below  the  transition  at  1 55  K  as  the  peak  positions  move  apart.  The  fitted  peak  positions  (blue  and  green  symbols]  are  also 
shown,  and  indicate  the  abrupt  nature  of  the  transition. 
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FIGURE  2:  Temperature  dependence  of  the  order  parameter  for  the  strongest  (1,0,3)  magnetic  reflection.  Data  were  obtained  in  diffraction  mode  using  the 
position-sensitive  detector.  The  solid  curve  is  a  simple  fit  to  mean  field  theory,  w/hich  indicates  that  the  order  develops  at  a  Neel  temperature  of  7^  =  1 37(3) 
K.  The  figure  on  the  right  shows  the  antiferromagnetic  structure  of  the  system,  giving  a  V2o^  xV2o.^  x2c^  magnetic  unit  cell.  All  the  magnetic  structures  for 
this  class  of  materials  have  antiparallel  spins  along  the  a  direction  and  parallel  spins  along  the  b  direction,  with  a>  b  and  the  spin  direction  along  a.  Both  the 
magnetic  and  structural  transitions  are  absent  in  superconducting  La(0,  FJFeAs  [T  =  26  K). 


At  170  K  detailed  refinements  of  the  data  obtained  on  BT-1 
show  that  the  crystal  structure  belongs  to  the  expected  tetragonal 
P4/nmm  space  group.  Upon  cooling  to  4  K,  the  (2,2,0)  reflection 
that  has  a  single  peak  at  170  K  spUts  into  two  peaks  as  shown 
in  Rg.  1.  This  immediately  demonstrates  that  a  structural  phase 
transition  has  occurred.  The  symmetry  of  the  low  temperature 
phase  can  be  described  by  an  orthorhombic  Cmma  space  group 
Avith  (approximate)  unit  cell  ^la-^  x^^la^  xc^  (the  subscript  T 
refers  to  the  tetragonal  parent  cell). The  detailed  temperature 
dependent  data  shown  in  Fig.  1  reveal  that  there  is  an  abrupt 
splitting  of  the  (2,2,0)  peak  that  occurs  at  ~  155  K.  These 
results  thus  indicate  that  the  non-superconducting  system  has  a 
structural  phase  transition  and  this  phase  transition  is  associated 
with  the  observed  resistivity  and  specific  heat  anomaUes.  Since  a 
similar  spUtting  of  the  (2,2,0)  peak  is  absent  in  superconducting 
La(OQ^2Foo8)-f'^-^^'       '^^'^  safely  assume  that  this  transition  is 
suppressed  with  the  appearance  of  superconductivity  in 
La(0^  ^FJFeAs  via  fluorine  doping. 

We  now  turn  to  the  search  for  magnetic  order  in  the  undoped 
system.  The  magnetic  scattering  is  expected  to  be  quite  weak, 
and  we  therefore  carried  high-intensity/coarse  resolution 
measurements  on  BT-7  to  search  for  spin-density-wave  magnetic 
order  in  LaOFeAs.  The  data  revealed  that  there  are  extra  peaks 
in  the  low  temperature  diffraction  pattern  that  are  simply  related 
to  the  chemical  unit  cell.  Indexing  these  reflections  indicates 
that  these  magnetic  peaks  arise  from  a  simple  stripe-type 
antiferromagnetic  structure  of  Fe  moments  with  a  magnetic  cell 
V2fl^  xVlfl^  x2c^  as  indicated  in  Fig.  2.  To  see  if  the  observed 
magnetic  scattering  at  low  temperature  in  LaOFeAs  is  associated 
with  the  150  K  phase  transition,  we  carried  out  a  measurement 


of  the  strongest  magnetic  peak.  Figure  2  shows  the  temperature 
dependence  of  the  square  of  the  ordered  magnetic  moment 
(proportional  to  the  magnetic  intensity),  which  vanishes  at 
~  137  K,  about  ~  18  K  lower  than  the  structural  phase  transition 
(Fig.  1).  That  the  magnetic  order  is  established  at  lower 
temperatures  than  the  structural  distortion  is  quite  surprising. 
Normalizing  the  magnetic  intensity  to  the  nuclear  scattering, 
we  find  an  ordered  Fe  moment  in  the  ground  state  of  0.36(5) 
The  small  value  of  the  moment  suggests  that  this  is  an  itinerant 
electron  (band)  magnetic.  Finally,  we  note  that  an  identical 
scan  in  the  superconducting  La(Og  ^^Fg  os^-^^-^^  shows  that  this 
magnetic  peak  is  absent.  Therefore,  the  resistivity  anomaly 
in  the  parent  compound  is  caused  by  structural  distortion, 
not  spin-density-wave  ordering  as  originally  suggested,  while 
both  transitions  are  suppressed  when  the  system  is  doped  into 
the  superconducting  regime.  The  disappearance  of  the  static 
antiferromagnetic  order  and  the  lattice  distortion  in  the  doped 
superconducting  materials  suggests  that  the  underlying  physical 
properties  of  this  class  of  superconductors  may  have  important 
similarities  to  the  high- 7^  copper  oxides.  In  any  event,  there  is 
no  doubt  that  this  new  class  of  materials  will  open  new  avenues 
of  research  regardless  of  the  origin  for  the  electron  pairing  and 
superconductivity. 
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Competing  Magnetic  Interactions,  Frustration,  and  the 
Structural  Phase  Transition  in  LaOFeAs 


T.  Yildirim^ 


The  layered  rare-earth  iron-pnictide  quaternary 
oxides  /?OFeAs  {R  =  La,  Ce,  Sm,  etc.)  exhibit 
superconductivity  up  to  a  transition  temperature  (T^)  of  55  K 
when  doped  with  either  electrons  (Oj  F  )  or  holes  (La^  Sr  ) 
[IJ.This  marks  the  first  non-copper  based  oxides  that,  like 
the  cuprates,  superconduct  at  relatively  high  temperatures 
upon  electron  and  hole  doping  of  their  non-superconducting 
parent  compounds.  Clearly,  understanding  the  electronic, 
magnetic,  and  structural  properties  of  the  parent  compound 
LaOFeAs  is  key  to  determining  the  underlying  mechanism 
that  makes  these  materials  superconduct  upon  electron/hole 
doping. 

Experimental  studies  including  resistivity  and  magnetic 
susceptibility  show  an  anomaly  near  150  K  in  LaOFeAs 
[IJ.The  origin  of  this  anomaly  has  been  very  recently 
determined  by  neutron  scattering  studies  [2].  It  has  been 
found  that  LaOFeAs  exhibits  spin-density-wave  (SDW) 
antiferromagnetic  long-range  ordering  with  a  small  moment 
(0.35  \Lg  per  Fe)  followed  by  a  small  structural  distortion  [2]. 
Even  though  the  SDW  ordering  was  predicted  from  Fermi 
surface  nesting,  there  is  no  proposed  microscopic  theory  that 
explains  the  origin  of  the  observed  structural  distortion.  It  is 
also  not  clear  if  the  magnetic  and  structural  phase  transitions 
are  related  to  each  other.  Finally,  given  the  tact  that  both  the 
cuprates  and  LaOFeAs  exhibit  antiferromagnetic  ordering, 
one  wonders  how  strong  and  what  kind  of  magnetic  spin- 
fluctuations  are  present  in  the  two  dimensional  (2D)  Fe- 
square  lattice  of  LaOFeAs. 

Very  recently,  from  accurate  all-electron  density  functional 
calculations  we  have  answered  some  of  these  questions 
[3].  We  considered  a  V2xV2-supercell  of  the  primitive  cell 
of  LaOFeAs  (see  Fig.l).  In  order  to  determine  the  true 
ground  state,  we  have  considered  four  different  cases.  These 
are  non-magnetic  (NM),  ferromagnetic  (F)  and  the  two 
different  antiferromagnetic  spin  configurations  shown  in 
Fig.2.  The  first  one  of  the  antiferromagnetic  configurations 
is  AFl  where  the  nearest  neighbor  spins  are  antiparallel  to 
each  other.  The  second  antiferromagnetic  configuration,  AF2, 
is  shown  in  Fig.  2b.  In  AF2  the  Fe  spins  along  the  square 


diagonal  are  aligned  antiferromagnetically.  This  is  the  stripe- 
phase  recently  observed  [2].  The  AF2  spin  configuration  can  be 
considered  as  two  interpenetrating  simple  square  AF  sublattices 
(red  and  blue  sublattices  in  Fig.  2b).  We  note  that  since  each  Fe 
ion  is  at  the  middle  of  a  square  AF  lattice,  the  mean  field  at  each 
spin  site  is  zero.  Hence  one  sublattice  can  be  rotated  freely  with 
respect  to  the  other  sublattice  without  costing  any  energy.  For  this 
reason  the  AF2  spin-configuration  is  fully  frustrated. 


As        ^  ^ 
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A+B 


FIGURE  1 :  The  crystal  structure  of  LaOFeAs  (left)  and  the  top  view  of  the 
FeAs-plane  and  the  relations  between  primitive  and  / 2xV 2-5upercell  used 
in  our  calculations.  The  dark  and  light  shaded  areas  indicate  the  As  atoms 
below  and  above  the  Fe-square  lattice,  respectively. 


(a)AF1  {J,>J,) 


(b)  AF2  {J,<J,) 


FIGURE  2:  Two  antiferromagnetic  configurations  considered  in  this  study. 

In  order  to  determine  which  spin  configurations  among  NM, 
F,  AFl,  and  AF2,  is  the  ground  state,  we  have  carried  out  total 
energy  calculations  as  a  function  of  Fe-magnetic  moment.  Our 
results  are  summarized  in  Fig.  3.  The  zero  of  energy  is  taken  as 
the  M  =  0  case  (i.e.,  NM  calculation).  From  Fig.  3,  it  is  clear  that 
LaOFeAs  has  only  one  magnetic  ground  state  which  is  AF2.  The 
Ferro  spin-configuration  always  results  in  the  highest  energy, 
regardless  the  Fe-magnetic  moment.  Similarly,  AFl  ordering 
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always  \ields  energies  higher  than  the  NM  case.  For  the  AF2 
ordering,  we  see  that  the  energy  minimum  occurs  near  the  fixed 
moment  calculation  with  M  =        Repeating  calculations  where 
magnetization  is  not  fixed,  we  obtained  the  optimum  magnetic 
moment  as  M  =  0.87[^  per  Fe.  As  we  discuss  below  in  detail,  the 
Fe  magnetic  moment  is  fiirther  reduced  almost  by  a  half  when  the 
structure  is  allowed  to  distort  due  to  AF2  stripe  ordering. 

In  order  to  gain  a  better  insight  into  the  nature  of  the  magnetic 
interactions  present  in  Fe-square  lattice  of  the  LaOFeAs  system, 
we  map  the  calculated  total  energies  of  the  F,  AFl  and  AF2 
configurations  shown  in  Fig.  3  a  to  a  simple  Heisenberg  like 
model.  For  fully  locahzed  spin-systems  this  is  a  perfect  thing  to 
do  but  for  the  case  of  LaOFeAs  this  is  only  an  approximation. 
Nevertheless,  the  calculated  Js  should  be  a  good  indication  of  the 
magnetic  interactions  present  in  the  system.  Fig.  3b  shows  and 

obtained  from  the  energies  of  the  F,  AFl  and  AF2  at  given 
magnetic  moment.  It  is  clear  that  both     and    are  quite  large 
and  positive  (i.e.,  antiferromagnetic).  /,  is  always  larger  than 
and  therefore  the  AF2  structure  is  the  only  ground  state  for  any 
given  moment  of  the  Fe  ion.  By  looking  at  the  exchange  paths  for 

and  /,,  we  notice  that  the  Fe-As-Fe  angle  is  around  75°  and 
120°  for  nn  and  nnn  Fe-pairs,  respectively.  Hence  it  makes 
sense  that  the  2nd  nn  exchange  interaction  is  as  strong  as  the 
nn  exchange  because  the  angle  is  closer  to  the  optimum  value  of 
180°.  It  is  quite  surprising  and  also  very  interesting  that  there 
are  strong  and  competing  antiferromagnetic  interactions  in  the 
LaOFeAs  system  that  results  in  a  totally  frustrated  AF2  spin 
configuration. 
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FIGURE  3:  Left:  The  total  energy  per  Fe  atom  versus  magnetic  moment  for 
F,  AFl  and  AF2  spin-configuration,  indicating  AF2  is  the  only  ground  state 
of  the  system.  Right:  The  magnetic  interactions  for  nn  and  nnn  Fe  ions 
obtained  from  the  energies  of  F,  AF,  and  AF2  configurations. 

We  next  discuss  the  imphcation  of  the  magnetically  frustrated 
AF2  configuration  on  the  structural  distortion  recently  observed 
by  neutron  scattering  [2].  Fig.  4  shows  the  total  energy  of  the 
AF2  spin  configuration  as  a  fiinction  of  the  y  angle.  When 
7  =  90°,  we  have  the  original  tetragonal  cell.  Once  y  deviates  from 
90°,  the  original  V2xV2-structure  (shown  as  a  dashed  Une)  is  no 
longer  tetragonal  but  orthorhombic  (i.e.,  the  cell  lengths  along 
the  a  and  I?  axes  are  no  longer  equal).  The  total  energy  versus  y 
plot  shown  in  Fig.4  clearly  indicates  that  the  structure  is  indeed 
distorted  with  y  =  91.0°,  which  is  in  good  agreement  with  the 


experimental  value  of  90.3°.  From  the  all-electron  method  we  get 
M  =  0.48[Xg  which  is  in  excellent  agreement  with  the  experimental 
value  of  0.35 fig.  The  net  energy  gain  by  the  structural  distortion 
is  about  12  meV  per  cell,  which  is  of  the  same  order  as  the 
temperature  at  which  this  phase  transition  occurs.  We  also 
considered  two  types  of  AF2  where  the  spins  along  the  short  axis 
are  aligned  parallel  or  antiparallel.  These  two  configurations  are  no 
longer  equivalent.  According  to  our  calculations  the  configuration 
in  which  the  spins  are  ordered  parallel  along  the  short-axis  is  the 
ground  state.  We  note  that  it  should  be  possible  to  determine 
the  exact  spin-configuration  from  neutron  powder  diffraction. 
Finally,  we  note  that  the  other  spin  configurations  including  the 
non-magnetic  case  do  not  yield  structural  distortion.  Therefore, 
the  experimentally  observed  structural  distortion  is  due  to  only 
AF2  stripe  ordering. 
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FIGURE  4:  The  total  energy  per  cell  versus  the  angle  y  for  NM,  F,  AFl 
and  AF2  spin-configurations.  Note  that  only  the  AF2  spin  configuration 
yields  structural  distortion.  The  inset  shows  that  as  y  increases,  the 
ferromagnetically  aligned  Fe  ions  (i.e.,  Fe1-Fe2)  get  closer  while  the 
antiferromagnetically  aligned  ions  (i.e.,  Fe1-Fe3)  move  apart. 

In  conclusion,  from  first-principles  fixed  spin  moment  calculations 
we  show  that  ferromagnetic  and  checkerboard  antiferromagnetic 
(AFM)  ordering  in  LaFeAsO  are  not  stable  and  the  stripe  AFM 
configuration  with  M(Fe)  =  0A8[i^  is  the  only  stable  ground 
state.  The  main  exchange  interactions  between  Fe  ions  are  large, 
antiferromagnetic  and  frustrated.  The  stripe  AFM  phase  breaks 
the  tetragonal  symmetry,  removes  the  frustration,  and  causes 
a  structural  distortion.  These  results  successfiiUy  explain  the 
magnetic  and  structural  phase  transitions  in  LaOFeAs  recendy 
observed  by  neutron  scattering.  The  presence  of  competing  strong 
antiferromagnetic  exchange  interactions  suggest  that  magnetism 
and  superconductivity  in  doped  LaOFeAs  may  be  strongly 
coupled,  much  Hke  in  the  high-TJ,  cuprates. 
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structural  Origin  of  the  Enhanced  Magnetostrictive  Properties 
of  Fe^^Ga^  (Galfenol) 

H.  Cao\"  P.  M.  Gehring^;  C.  P.  Devreugd,  J-F.  Li,  D.  Viehland^;  J.  A.  Rodriguez'^ 


Galfenol  is  the  name  given  to  alloys  of  Ga  and  Fe 
created  by  the  Naval  Ordnance  Labs  in  1998  that 
exhibit  very  large  magnetostriction  and  a  high  relative 
permeability  coupled  with  robust  mechanical  properties. 
Magnetostrictive  materials  are  the  magnetic  analogue 
of  piezoelectrics  in  that  they  can  be  used  to  convert 
between  electrical  and  mechanical  forms  of  energy 
via  their  mechanical  response  to  an  applied  magnetic 
field.  Although  Galfenol  has  only  1/3^''  to  1/4*  the 
magnetostriction  of  Terfenol-D,  an  alloy  based  on  Tb, 
Dy,  and  Fe  that  exhibits  the  highest  room-temperature 
magnetostriction  of  any  known  material,  it  is  nonetheless 
an  extremely  attractive  choice  for  use  in  device 
applications  such  as  magneto-acoustic  sensors,  actuators, 
and  transducers  because  it  can  be  used  in  mechanically 
harsh  environments  with  minimal  shock  hardening. 

The  addition  of  Ga  into  the  body-centered  cubic  (bcc) 
a-Fe  (or  yl2)  phase  is  known  to  produce  a  diversity 
of  crystal  structures.  In  fact,  both  equilibrium  and 
metastable  Fe^  Ga^  phase  diagrams  have  been  reported 
[1].  In  the  metastable  case  a  chemically  ordered  bcc 
Fe,Ga  (or  DO,)  phase  coexists  with  the  A2  phase  for 
0.14  <  X  <  0.20  [4],  where  both  the  magnetostriction  along 
the  [100]  direction  (Aj^q)  and  the  stability  of  the  two- 
phase  region  are  strongly  affected  by  the  rate  of  cooling. 
The  presence  of  two  peaks  in  the  magnetostriction  as  a 
function  of  Ga  content  has  been  reported  for  alloys  that 
have  been  cooled  rapidly  from  high  temperatures.  The 
first  peak  at  x  ~  0.19  has  been  attributed  to  an  increase 
in  the  magneto-elastic  coupling  that  results  from  the 
formation  of  short-range  ordered  (SRO)  Ga  pairs  along 
the  [100]  axis  of  the  ^2  structure,  while  the  second  peak 
at  X  ~  0.27  has  been  attributed  to  a  softening  of  the  shear 
modulus  c'  =  (fjj-fj-,)/2.  The  history  dependence  of  both 
Ajgg  and  the  A2-D0^  two-phase  region  suggests  that  the 
enhanced  magnetostriction  may  be  due  to  an  underlying 
heterogeneity  rather  than  being  due  to  a  conventional 
homogeneous  ferromagnetic  phase. 
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FIGURE  1 :  Neutron  scattering  intensity  map  measured  at  300  K  on  a 
furnace-cooled  single  crystal  of  Fe^g^Ga^^g.  The  symmetry  of  the  underlying 
matrix  allows  no  (1,0,0)  reflection.  The  broad  diffuse  scattering  that  peaks 
near  (0.95,0,0)  thus  indicates  the  presence  of  a  short-range  ordered 
structure  having  a  different  symmetry  and  lattice  spacing. 


(0,  0,  L)  (r.l.u) 


FIGURE  2:  Neutron  scattering  intensity  map  measured  at  300  K  on  a 
slow-cooled  single  crystal  of  Fejjg.,GaQ^g.  Under  slow-cooled  conditions, 
Fepg^Ga^^ig  exhibits  a  long-range  ordered  structure  with  DO^  symmetry,  for 
which  a  (1,0,0)  reflection  is  allowed.  The  diffuse  scattering  observed  in  the 
quenched  crystal  is  absent.  The  diffuse  scattering  located  below  the  line 
H  =  1  and  shown  in  green  is  believed  to  be  size-effect  scattering. 

Figure  1  shows  a  color  contour  plot  of  the  neutron  diffuse 
scattering  intensity  measured  at  room  temperature  in  the 
vicinity  of  the  (1,0,0)  reciprocal  lattice,  position  from  a  furnace 
cooled  (10  °C/min)  single  crystal  of  Fegg^Ga^     a  composition 
corresponding  to  one  of  the  two  peaks  in  the  magnetostriction 
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[2].  The  presence  of  substantial  diffuse  scattering,  slightly 
offset  from  (1,0,0),  is  evident  and  is  shown  in  orange  and 
yellow.  This  broad  scattering  is  consistent  with  the  presence 
of  a  short-range  ordered  structure  having  a  symmetry  and 
lattice  spacing  different  from  that  of  the  underlying  ^2 
matrix,  which  admits  no  structural  peak  at  (1,0,0).  Identical 
neutron  scattering  measurements  made  on  a  slowly-cooled 
(2  °C/min)  single  crystal  of  the  same  composition,  for 
which  the  magnetostriction  is  markedly  lower,  are  shown  in 
Fig.  2.  In  this  case  a  strong  and  sharp  Bragg  peak  is  seen  at 
(1,0,0),  consistent  with  a  long-range  ordered  DO^  structure. 
More  importantly,  no  evidence  is  seen  of  the  broad  and 
nearly  isotropic  diffuse  scattering  shown  in  Fig.  1.  Instead 
anisotropic  diffuse  scattering  contours  are  observed  below  the 
line  H  =  1,  which  are  believed  to  arise  from  correlated  atomic 
displacements  that  result  from  the  atomic  size  mismatch 
between  Fe  and  Ga,  also  known  as  size  effect  diffuse  scattering. 
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FIGURE  3;  Diffuse  scattering  profiles  measured  at  300  K  near  (1,0,0]  along 
[1 ,0,0]  are  shown  for  furnace-cooled  Fe.,^Ga^  [x  =  0.1 0,  0,1 5,  0.1 9].  The 
peak  and  0-integrated  diffuse  scattering  intensities  grow  with  increasing  Ga 
content  and  track  the  magnetostriction,  which  reaches  a  maximum  value  for 
X  =  0.1 9.  The  sharp  peaks  at  (1 ,0,0)  and  (0.67,0,0)  come  from  neutrons 
with  wavelengths  A/2  and  A/3,  respectively,  that  scatter  from  the  (2,0,0) 
reflection,  and  vanish  after  inserting  a  third  HOPG  filter  before  the  sample 
as  shown  in  the  inset. 

To  ascertain  whether  or  not  the  diffuse  scattering  observed  in 
Fig.  1  is  related  to  the  enhanced  magnetostriction  observed 
in  quenched  alloys  of  Galfenol,  linear  Q-scans  through  the 
diffuse  scattering  intensity  along  the  [100]  direction  were 
made  on  three  different  compositions  of  Fe^^Ga^  with  very 
different  values  of  magnetostriction.  Data  were  obtained  on 
single  crystals  for  x  =  0.10,  0.15,  and  0.19,  which  are  shown  in 
Fig.  3.  It  is  immediately  apparent  that  the  diffuse  scattering 


grows  markedly  with  increasing  Ga  content  and  in  fact  tracks 
the  corresponding  increase  in  the  magnetostriction.  These 
data  are  consistent  with  a  structurally  heterogeneous  model  [3] 
that  was  proposed  to  explain  the  enhanced  magnetostriction 
in  Galfenol  alloys.  This  model  assumes  that  the  quenched 
condition  is  in  fact  a  structurally  and  chemically  heterogeneous 
state  consisting  of  coarsening-resistant  DOj  nanoprecipitates 
within  an  A2  matrix.  An  essential  element  of  this  model  is 
that  the  nanoprecipitates  possess  a  lower  than  cubic  symmetry. 
Close  examination  of  the  data  in  Fig.  3  reveals  an  asymmetry 
in  the  diffuse  scattering,  especially  for  x  =  0.19,  that  is  well 
described  by  tetragonally-distorted  DO,  regions  with  a 
tetragonality  c/a  =  1.2.  This  value  agrees  precisely  with  that 
predicted  by  the  structurally  heterogenous  model.  This  finding 
lends  strong  credence  to  idea  that  the  large  magnetostriction 
in  quenched  Galfenol  alloys  is  directly  related  to  the  structural 
heterogeneity  of  the  -DO,  nanoprecipitates,  i.e.,  that  the 
enhancement  ofl^^^^  results  from  the  redistribution  of  the 
orientation  of  lower  symmetry  DO,-like  nanoprecipitate 
particles  under  field  or  stress. 
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The  magnetocaloric  effect  (MCE) — the  expulsion 
of  heat  from  a  material  when  a  magnetic  field  is 
applied  and  its  absorption  when  the  field  is  removed — was 
discovered  in  1881.  In  1933  using  a  technique  based  on 
MCE  called  "adiabatic  demagnetization"  a  temperature 
well  below  1  K  was  achieved  for  the  first  time,  and  for 
most  of  its  history  the  effect  was  only  useful  for  producing 
very  low  temperatures.  Recently,  magnetic  refrigeration  at 
ambient  temperatures  has  attracted  renewed  interest  with 
the  discovery  of  new  materials  with  improved  efficiencies 
and  advantages,  as  potential  replacements  for  the  classical 
vapor  compression  systems  in  use  today.  To  be  useful 
the  magnetorefrigerant  needs  to  have  a  large  change  in 
magnetic  entropy  near  room  temperature  for  modest 
applied  fields  (<  2  T),  be  inexpensive,  and  environmentally 
friendly,  and  Mnj^^Fej^Pj_^Ge^  is  the  first  system  to  fulfill 
all  these  requirements. 

We  have  used  neutron  diffraction  and  magnetization 
measurements  to  reveal  that  the  application  of  a  magnetic 
field  drives  the  material  through  a  first  order  phase 
transition  from  paramagnetic  (PM)  to  ferromagnetic 
(FM),  and  that  these  phases  correspond  to  two  very 
distinct  crystal  structures.  The  latent  heat  of  this  transition, 
deriving  both  from  the  sudden  onset  of  magnetic  order 
and  change  in  crystal  structure,  is  the  source  of  the  MCE 
[1].  This  entropic  control  of  combined  magnetic  and 
structural  transitions  circumvents  the  limitation  imposed  by 
the  magnetic  entropy  alone,  and  provides  a  new  pathway  to 
optimize  the  refrigerant  properties.  Indeed,  by  tuning  the 
physical  properties  of  this  system  we  now  have  achieved  an 
MCE  of  ~  75  J/Kg  K  on  both  increasing  and  decreasing 
field  cycles,  more  than  twice  the  value  of  the  previous 
record,  thrusting  Mn,  Fe,  P,  Ge  into  the  lead  as  a  viable 
commercial  refrigerant. 

Neutron  powder  diffraction  is  the  classic  technique 
to  solve  for  both  crystal  and  magnetic  structures, 
and  for  Mn^  jFe^  ^P^  gGe^  ^  we  find  it  is  single  phase 
and  paramagnetic  above  ~  260  K,  single  phase  and 


ferromagnetic  below  200  K,  and  in  between  undergoes  a 
strongly  first-order  phase  transition  as  a  function  of  temperature 
or  apphed  magnetic  field  (see  Fig.  1).  It  turns  out  that  both 
phases  possess  the  same  space  group  symmetry  {P62m)  but 
have  distinctly  different  structures,  with  the  a-  and  i^-axes  being 
-1.3  %  longer  while  the  f-axis  is  contracted  by  «  2.6  %  in  the 
FM  phase  compared  to  the  PM  phase.  There  is  no  significant 
change  in  the  unit  cell  volume  for  the  two  phases.  We  find  that 
one  site  (3^)  (see  upper  panel  of  Fig.  2)  is  completely  occupied 
by  Mn  atoms,  which  are  co-planar  with  P/Ge(l)  atoms  at  the 
lb  site  in  the  z  =  0.5  layer.  The  other  site  {3f)  has  ~  93  %  Fe, 
with  ~  7  %  Mn  distributed  randomly,  and  the  3/  Fe/Mn  site 
is  co-planar  with  the  P/Ge(2)  atoms  at  the  2c  site  in  the  z  =  0 
layer.  Ge  and  P  are  randomly  mixed,  although  Ge  atoms  prefer 
the  P/Ge(2)  {2c)  site  («  27  %  Ge  occupied)  to  the  P/Ge(l)  (U) 
site  (~  5  %  Ge).  For  the  ferromagnetic  phase  the  refinements 
give  Mn  and  Fe  moments  parallel  in  the  a-l>  plane,  with  ordered 
moments  (at  245  K)  of  2.9(1)     and  0.9(1)  ^i^  for  the  Mn  (3^)  ' 
site  and  Fe/Mn  {3f)  site,  respectively,  similar  to  what  is  seen  for 
other  compounds  with  the  Fe,P-type  structure. 

The  important  point  is  that  the  results  demonstrate  that 
magnetic  field,  or  temperature,  has  no  significant  effect  other 
than  to  convert  the  system  between  the  ferromagnetic  and 
paramagnetic  structural  phases. 

The  results  make  it  clear  that  the  large  MCE  originates  from 
the  applied  magnetic  field  driving  the  phase  conversion,  and  the 
measured  magnetic  entropy  changes  are  shown  in  Fig.  2.  For 
magnetic  field  increasing  and  decreasing  we  obtain 
74  J/Kg  K  and  78  J/Kg  K,  respectively,  for  a  field  change  of 
5  T  in  this  bulk  Mn^  jFe,|£,P,^gGeg2  compound.  These  values  are 
twice  the  previous  value  for  this  system  and  the  highest  MCE 
for  any  material  presently  known.  It  is  also  quite  unusual  to  be 
able  to  achieve  similar  values  in  both  field  directions,  as  typically 
the  field-decreasing  MCE  can  be  quite  drastically  reduced.  For 
the  present  material  the  field-decreasing  MCE  actually  exceeds 
the  field-increasing  values. 
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FIGURE  1:  Comparison  of  the  diffraction  and  magnetization  measurements. 
(A)  Temperature  dependence  of  the  magnetization  in  an  applied  field  of 
0.1  T.  [B)  Integrated  intensities  of  the  (001 )  reflections  for  the  PM-  and 
FM-  phases  as  a  function  of  temperature  on  cooling  and  warming.  (C) 
Field-dependent  magnetization  at  255  K.  The  magnetic  entropy  change 
|ASm|,  normalized  to  the  magnetization,  is  shown  for  comparison.  (□)  Field 
dependence  of  the  integrated  intensities  of  the  (001 )  reflections  for  the 
PM-phase  and  FM-phase  at  255  K,  showing  that  the  FM-phase  fraction 
tracks  the  magnetization  data  and  normalized  magnetic  entropy  change 
|ASm|.  (E)  Fraction  of  the  FM-phase  at  255  K  as  the  field  increases.  The  FM- 
phase  fraction  increases  smoothly  to  =  86  %.  Data  normalized  from  |ASm| 
are  shown  for  comparison.  (F)  |ASm|  increases  linearly  with  the  FM-phase 
fraction.  |ASm|  is  projected  to  be  =  103  J/Kg  K  if  the  transition  went  to 
completion  for  this  sample. 

The  present  results  directly  demonstrate  that  the  transition 
from  the  PM-phase  to  the  FM-phase  and  the  associated  huge 
magnetocaloric  effect  are  directly  controlled  by  the  first-order 
structural  phase  transition  between  these  two  phases,  and, 
moreover,  that  a  completed  phase  conversion  will  increase  the 
MCE  up  to  «  100  J/Kg  K  in  a  field  of  5  T  in  this  system.  The 
transition  can  be  driven  by  temperature  or  applied  magnetic 
field,  and  for  use  as  a  magnetic  refrigerant  the  field-dependent 
properties  are  critical.  The  improved  properties  and  overall 
advantages  of  this  material  open  the  possibility  for  its  use  in 
magnetic  refrigerant  applications. 
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FIGURE  2:  The  Magnetocaloric  Effect.  Temperature  dependence  of  the 
entropy  change  of  the  bulk  Mn^  i'^^os'^ob'^^os  compound  as  a  function 
of  applied  magnetic  field  up  to  5  T.  The  upper  panel  shows  the  crystal 
structure  for  the  FM-phase  which  differs  from  that  in  the  PM-phase  in  the 
ordering  of  the  spins  and  the  lengths  of  a,  b,  and  c. 

This  material  also  exhibits  a  large  magnetoelastic  effect — the 
change  in  the  lattice  constants  with  magnetic  field — and  this 
useflil  property  also  originates  from  the  structural  changes 
associated  with  the  first-order  transition.  Furthermore, 
our  investigations  suggest  a  number  of  ways  in  which  the 
magnetothermal  properties  can  be  tuned  by  chemical  doping 
or  mechanical  treatment.  For  example,  we  have  found  that  an 
applied  field  can  induce  substantial  preferred  crystallographic 
orientation,  which  results  in  improved  field-dependent  properties. 
We  also  expect  that  the  magnetic  and  structural  properties  can 
be  fiirther  optimized  by  selective  chemical  substitutions  and 
processing  techniques.  The  time  dependence  associated  with 
the  transition  does  tend  to  reduce  the  practical  MCE  available 
for  applications,  but  the  overall  improvements  and  prospects  for 
fiirther  advances  in  performance  make  this  material  the  magnetic 
refrigerant  of  choice  to  enable  a  wide  range  of  commercial 
magnetorefrigerant  applications. 
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Systems  that  are  simultaneously  magnetically  ordered 
and  ferroelectric — multiferroics — are  very  rare  in  nature. 
The  reason  is  that  typical  ferroelectrics  such  as  BaTiO^ 
develop  their  electric  polarization  via  hybridization  of 
empty  (j'-orbitals  on  the  Ti  ion  with  occupied  ^-orbitals  of 
the  octahedraUy  coordinated  oxygen  ions,  which  causes  the 
off-centering  displacement  of  the  Ti  ion  and  hence  the  net 
electric  polarization.  Magnetism  such  as  for  iron,  on  the 
other  hand,  requires  electronically  occupied  a'-orbitals,  and 
thus  magnetism  and  ferroelectricity  are  mutually  exclusive 
phenomena  through  the  standard  ferroelectric  mechanism. 
Nevertheless,  the  tremendous  technological  possibilities  that 
would  be  enabled  with  a  material  in  which  the  magnetic 
state  can  be  switched  with  an  electric  field  has  driven 
research  in  multiferroics,  and  the  recent  discovery  of  a  new 
class  of  multiferroic  materials  that  feature  strongly-coupled 
spontaneous  magnetic  and  electric  order  [1]  has  reignited 
this  dream  ol  materials  engineers. 

The  novelty  of  these  new  multiferroic  materials  is  that 
ferroelectricity  is  generated  by  the  magnetic  structure  itself, 
and  so  ferroelectricity  can  be  easily  switched  when  the 
magnetic  structure  is  tuned  by  magnetic  fields  [2,3].  More 
importantly,  here  we  show  that,  reciprocally,  electric  fields 
can  also  be  used  to  control  magnetic  properties,  and  that  the 
magnetic  chirality  or  handedness  can  be  directly  controlled 
with  the  electric  field — ^the  first  step  towards  electric  field 
control  of  the  macroscopic  magnetic  properties  for  device 
applications. 

The  model  multiferroic  system  that  was  chosen  to 
demonstrate  this  property  is  Ni^V^O^,  whose  magnetic 
and  electric  properties  we  have  extensively  characterized 
previously  [4-6].  The  magnetism  in  Ni^V^Og  arises  from 
weakly-coupled,  buckled  Kagome  planes  of  5  =  1  quantum 
spins  that  are  localized  on  the  Ni^*  sites,  while  the  V  ions 
are  nonmagnetic.  Ni^V^Og  is  an  insulating  material  (as 
all  ferroelectrics  must  be)  and  the  spin  interactions  are 


mediated  through  electronic  orbital  overlap  and  the  exchange 
of  electrons  in  quantum  mechanically  entangled  states.  The 
corner-shared  geometry  of  Ni-*  triangles  and  antiferromagnetic 
next-nearest  neighbor  interactions  leads  to  strong  magnetic 
fluctuations  even  below  the  temperature  where  simpler  magnets 
would  order,  and  the  competition  of  these  interactions  stabilizes 
a  series  of  distinctly  different  ordered  structures  with  decreasing 
temperature. 


Buckled  Kagome  sheet 

Electric  Polarization 


FIGURE  1:  (top)  Single  layers  of  a  Ni^*  Kagome  lattice  that  are  buckled 
along  the  c-axis,  and  that  in  NijV^Og  are  stacked  along  the  b-axis.  (Center/ 
Bottom)  spiral  magnetic  structure  that  is  only  possible  because  of  a  crystal 
distortion  that  also  leads  to  macroscopic  electric  polarization,  either  up  or 
down  depending  on  the  sense  of  rotation  of  the  spiral. 

Qiiite  generally,  all  magnetic  materials  attempt  to  minimize 
magnetic  disorder  or  entropy  at  low  temperature.  This  is  usually 
achieved  with  a  static  long-range  order  of  the  quantum  degrees 
of  freedom  such  as  the  spin.  For  incommensurate  magnetic 
order,  this  is  best  done  in  a  spin  spiral  structure  with  a  maximally 
ordered  moment  on  each  site.  However,  in  Ni3V,Og,the  crystal 
structure  and  the  interactions  favor  incommensurate  magnetic 
order  for  which  the  size  of  the  ordered  magnetic  moment  varies 
sinusoidally  from  one  Ni-*  site  to  next,  leading  to  a  great  amount 
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of  magnetic  disorder  near  the  nodes  of  the  wave.  One  way  to 
lower  the  magnetic  entropy  further  is  via  a  crystal  distortion, 
and  we  have  previously  shown  (see  Fig.  1)  that  Ni^V^O^  is  one 
of  the  few  knovra  materials  where  the  resulting  crystal  distortion 
involves  a  macroscopic  electric  polarization  that  is  magnetically 
driven. 

The  ferroelectric  polarization  in  NijV,Oj,  can  be  turned  on  and 
ofFby  simply  applying  a  magnetic  field,  because  a  magnetic  field 
induces  different  magnetic  structures,  some  of  which  allow  an 
electric  polarization.  The  more  interesting  question  is  whether 
the  inverse  coupling  can  also  be  achieved  in  Ni^V^Og,  namely,  can 
the  magnetic  structure  be  controlled  with  an  electric  field? 

To  test  such  a  scenario,  we  prepared  a  single  crystal  of  Ni,V,0;^ 
that  was  a  1  mm  thin  platelet  normal  to  the  ferroelectric  axis.  A 
capacitor  was  made  by  evaporating  an  electrode 
[(4  nm  Cr)/(40  nm  Au)  layer]  on  each  of  the  two  faces,  and 
then  attaching  Au  wires  using  silver  epoxy.  The  sample  was 
held  by  a  sapphire  plate  in  a  crv'ostat.  During  cool-down  an 
external  electric  field  was  applied  using  an  electrometer  to  obtain 
a  macroscopic  electrical  polarization.  We  also  surveyed  in-situ 
the  ferroelectric  polarization  during  the  neutron  diffraction 
measurements,  which  were  performed  on  BT-7  utilizing  large 
^He  polarizer  cells  to  achieve  the  neutron  polarization  and 
analysis.  The  neutron  polarization  P  was  chosen  to  be  parallel 
to  the  momentum  transfer  wave-vector  Q  under  investigation, 
[P|  1^].  In  this  case  the  neutron  polarization  is  reversed,  or 
"flipped",  when  the  neutrons  magnetically  scatter  from  the 
sample,  while  structural  Bragg  peak  scattering  never  flips  the 
neutron  spin.  Then  the  magnetic  scattering  is  observed  when 
the  incident  neutron  spin  is  either  parallel  (+)  and/or  antiparallel 
(-)  to  the  momentum  transfer  wave-vector  Q,  and  since  the  spin 
is  reversed  on  scattering  we  measure  either  the  (+  -)  or  (-  +) 
magnetic  cross  sections. 

For  typical  magnetic  materials  these  two  cross  sections 
would  have  equal  intensities,  but  in  the  ferroelectric  phase  of 
NijV^Og  the  incommensurate  magnetic  order  has  a  chirality, 
or  handedness.  The  modulation  wave-vector  is  along  the  a 
axis,  (^',0,0)  where  q  ~  0.265,  and  hence  the  magnetic  Bragg 
peaks  occur  at  satellite  positions  such  as  g  =  (l+y,l,l).  Our 
experimental  results  are  shown  in  Fig.  2,  where  we  see  that  the 
spin-flip  neutron  diffraction  intensity  of  the  magnetic 
Bragg  peak  depends  on  the  incoming  neutron  spin  polarization. 
For  an  applied  electric  field  of  +300  kV/m  the  Bragg  peak  is  only 
observed  for  the  (-  +)  configuration;  the  small  peak  observed  for 
the  (+  -)  configuration  is  simply  due  to  the  imperfect  polarization 
of  the  instrument.  More  importantly,  by  reversing  the  electric 
field  the  magnetic  Bragg  peak  is  only  observed  in  the  (+  -) 
configuration,  demonstrating  in  a  simple  and  elegant  way  that  the 
magnetic  structure  in  Ni^V^Og  can  be  tuned  with  an  electric  field. 


It  can  be  shown  that  the  reversal  of  the  neutron  diffraction 
intensity  with  electric  field  results  from  an  electrically  induced 
switching  of  the  magnetic  chiraUty — the  "handedness"  of  the 
spiral.  It  is  believed  that  the  microscopic  coupling  between 
the  magnetic  structure  and  the  ferroelectric  polarization  is 
mediated  by  the  Dzyaloshinskii-Moriya  antisymmetric  exchange 
interaction  of  the  form  D  x  (5,  x  S.),  where  D  changes  sign  when 
the  electric  polarization  is  switched.  Switching  the  electric 
polarization  of  the  ferroelectric  domain  thus  changes  the  sense  of 
rotation  of  the  spin  spiral. 
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FIGURE  2:  Full  polarized  neutron  magnetic  diffraction.  The  neutron 
polarization  Pis  parallel  to  Q,  [P|  Q],  so  that  all  the  magnetic  scattering  is 
spin-flip,  [a)  A  +300  kV/m  electric  field  selects  one  chiral  domain,  and  then 
the  magnetic  peak  is  only  observed  in  the  (-  +]  configuration;  the  small  peak 
In  the  (+  -)  configuration  is  due  to  the  imperfect  instrumental  polarization, 
(b)  With  the  ferroelectric  polarization  reversed  by  applying  -300  kV/m,  the 
magnetic  Bragg  peak  is  only  observed  in  the  (+  -]  configuration. 

Ni^V,Og  illustrates  the  great  potential  for  multifunctionality  that 
can  be  obtained  from  materials  with  competing  interactions. 
The  most  important  lesson  to  draw  from  studies  of  these 
magnetically-induced  ferroelectrics  is  that  the  competing 
magnetic  interactions  not  only  lead  to  novel  magnetic  ground 
states  of  interest  to  the  community  of  condensed  matter  scientists, 
but  that  they  also  have  the  potential  to  be  key  components  for  the 
design  of  next-generation  electronics,  data  storage  and  computing 
applications. 
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Killing  Cancer  Cells  with  Magnetic  Nanoparticles 


C.L.  Dennis\  AJ.  JacksonH  J.A.  Borchers^  R.  Ivkov^  A.R.  Foreman^ 
and  C.  Gruttner^ 


Fevers  offer  an  effective  system-wide  defense  against 
pathogens  like  bacteria  by  disrupting  their  biochemical 
processes,  rendering  the  pathogens  more  susceptible  to  attack 
by  the  immune  system.  However,  if  the  fever  gets  too  high 
(>  45  °C  in  humans),  then  damage  occurs  to  healthy  tissue, 
and  if  left  too  long  at  elevated  temperatures,  the  healthy 
tissue  will  die.  In  an  analogous  manner,  cancer  cells  can 
be  treated  with  heat  to  disrupt  their  biochemical  processes. 
A  key  difference  is  that  cancer  cells  start  to  experience 
damage  at  42  "C,  which  leaves  a  3  "C  temperature  window  in 
which  cancers  can  be  killed  with  little  or  no  damage  to  the 
surrounding  healthy  tissue.  Similar  to  a  fever,  nanoparticle 
mediated  hyperthermia  is  a  prospective  cancer  therapy  that 
destroys  tumors  by  locally  heating  tumor  cells. 

Although  hyperthermia  is  not  a  new  idea,  there  are 
two  distinguishing  aspects  of  nanoparticle  mediated 
hyperthermia  as  compared  to  "conventional"  hyperthermia 
(i.e.,  a  fever):  1)  the  heat  is  generated  locally  at  the  tumor 
cells  by  the  magnetic  nanoparticles  subjected  to  an 
alternating  magnetic  field,  and  2)  the  temperature  change 
can  be  correlated  with  the  energy  deposited  and  therefore 
to  the  amount  or  "dose"  of  nanoparticles  used.  The  primary 
benefit  of  this  technology  is  a  significant  reduction  in  the 
severity  and  duration  of  side  effects  that  are  common  with 
conventional  cancer  treatments  like  chemo-  and  radiation 
therapy  because  of  the  reduction  in  damage  to  surrounding 
healthy  tissue. 

In  designing  these  nanoparticle  systems,  the  primary 
requirement  is  a  biologically  stable  material  that  allows 
control  of  the  total  amount  of  heat  deposited  in  order  to 
kill  the  tumor  while  minimizing  damage  to  the  surrounding 
healthy  tissue.  To  design  more  effective  systems,  NIST 
researchers  in  collaboration  with  their  industrial  partners 
examined  two  candidates  with  significantly  different  heating 
characteristics.  Contrary  to  accepted  practice,  we  found 
that  strongly  interacting  systems  heat  significantly  better 
than  weakly  or  non-interacting  systems,  with  the  degree  of 
interaction  controlled  by  the  spacing  between  nanoparticles. 


FIGURE  1:  Schematic  of 
the  magnetic  nanoparticle 
developed  for  hyperthermia 
treatment  used  in  these 
studies:  The  iron  oxide  core 
generates  the  heat  while  the 
surrounding  dextran  shell 
stabilizes  the  nanoparticles  in 
solution. 


The  nanoparticles  studied  consist  of  a  magnetic  iron  oxide  core, 
surrounded  by  a  dextran  shell,  with  the  dextran  controlling  the 
spacing  between  of  nanoparticles  (see  Fig.  1).  The  difference 
between  the  two  systems  examined  lies  in  the  dextran  shell  (the 
cores  are  similar):  the  "single  dextran"  sample  was  coated  once 
with  dextran  whereas  the  "double  dextran"  sample  was  coated 
twice.  To  properly  compare  these  two  systems  their  physical  and 
magnetic  properties  were  extensively  characterized.  First,  the 
iron  oxide  core  of  both  systems  was  determined  by  Mossbauer 
spectroscopy  to  be  composed  of  magnetite  (Fe^O^),  except  for 
<  2  %  of  iron  hydroxide.  This  hydroxide  most  likely  originates 
from  the  storage  of  the  nanoparticles  in  water  and  is  assumed 
to  be  limited  to  the  surface  of  the  nanoparticle  core.  Analytical 
Ultra-Centrifiigation  (AUC)  yielded  a  size  distribution  of 
(44  ±13)  nm  for  the  nanoparticle  cores  for  the  double  dextran 
sample.  Transmission  electron  microscopy  images  also  confirm 
a  core  diameter  ~  50  nm  for  both  samples.  Photon  Correlation 
Spectroscopy  (PCS)  yielded  a  combined  size  and  size  distribution 
for  the  iron  oxide  core  and  dextran  together  of  (92  ±  14)  nm 
for  the  double  dextran  sample  and  (96  ±  32)  nm  for  the  single 
dextran  sample.  We  can  thus  estimate  the  dextran  thickness  to  be 
~  25  nm  for  both  samples  -  in  good  agreement  with  the  40,000 
Da  dextran  used.  This  may  indicate  that  the  addition  of  the 
second  dextran  layer  increases  the  density  of  the  dextran  coating. 

From  the  magnetization  vs.  field  (hysteresis)  loops  at  room 
temperature  of  the  two  colloidal  systems  (normalized  to  the  mass 
of  particles  present  in  the  colloid),  the  most  prominent  result  is 
that  the  saturation  magnetization  of  the  double  dextran  sample 
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is  (41.08  ±  0.03)  kA  m-/g,  33  %  less  than  that  of  the  single 
dextran  sample  which  is  (61.64  ±  0.03)  kA  m'/g.  Other  than  this 
difference  in  magnitude,  the  shapes  of  the  two  hysteresis  loops  are 
nearly  identical  (see  Fig.  2). 

The  Specific  Absorption  Rate  (SAR)  values  were  measured 
for  //=  86  kA/m  (1080  Oe)  and/=  150  kHz  using  colloids 
of  nominaUv  equal  concentrations  which  are  corrected  for  the 
thermal  properties  of  the  calorimeter,  coil,  and  water  and  then 
normalized  to  iron  concentration.  (The  iron  concentration 
is  different  for  the  two  samples,  but  the  variation  in  SAR 
with  concentration  for  each  sample  is  about  10  %  for  the 
concentrations  probed  here.)  Here  we  see  the  most  striking 
difference  -  the  double  dextran  sample  has  a  measured  SAR  of 
1075  W/(g  of  Fe)  while  the  single  dextran  has  a  measured  SAR 
of  537  W/(g  of  Fe)  -  a  factor  of  2  smaller. 
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FIGURE  2:  Magnetization  vs.  field  loops  show  the  same  shape  except  for  a 
33  %  difference  in  saturation  magnetization  between  the  double  and  single 
dextran  samples. 

Since  the  SAR  has  the  opposite  trend  from  what  we  would 
expect  based  on  the  saturation  magnetization  data,  the  question 
is:  if  the  nanoparticles  are  physically  and  magnetically  nearly 
identical,  what  is  the  origin  of  this  difference?  Only  Small  Angle 
Neutron  Scattering  (SANS)  and  Ultra- SANS  experiments  could 
directly  probe  the  separation  between  the  nanoparticles  to  show 
that  the  effect  of  the  dextran  layer  is  to  change  the  hard  sphere 
interaction  radius  (separation)  from  >  698  nm  for  the  single 
dextran  layer  to  157  nm  for  the  double  dextran  layer  -  a  factor  of 
four  smaller  (see  Fig.  3).  This  decrease  in  the  interaction  radius 
for  the  double  layer  is  not  due  to  a  change  in  average  diameter  (as 
determined  earlier  by  both  AUC  and  PCS),  nor  by  a  change  in 
volume  fraction  [1]. 

This  smaller  separation  in  the  double  dextran  sample  would 
have  a  two-fold  effect:  (1)  the  dipolar  interactions  would  be 
significantly  stronger,  enabling  the  nanoparticles  to  couple 
their  behavior  under  an  alternating  field  (thereby  amplifying 


the  heating),  and  (2)  the  smaller  interaction  radius  would  mean 
that  more  particles  are  grouped  closer  together,  enhancing  the 
local  heat  output  in  a  smaller  area.  This  increased  heating  effect 
shows  significant  enhancement  in  the  efficacy  of  these  particles 
in  pre-cHnical  (mouse)  trials  at  Dartmouth  College  [2],  where 
the  doubling  time  of  the  tumor  -  an  aggressive  mammary  tumor 
-  after  one  treatment  increases  to  an  average  of  60  d  from  an 
average  doubling  time  of  2  d  for  an  untreated  tumor. 

In  conclusion,  we  have  recently  determined  that  strongly 
interacting  coated  magnetic  nanoparticles  yield  large  increases  (by 
a  factor  of  two)  in  heat  output  as  compared  to  nominally  identical 
particles  which  are  only  weakly  interacting.  This  heating  enables 
larger  localized  temperature  changes  that  can  kiU  cancer  cells 
without  damaging  the  surrounding  healthy  tissue.  This  finding  is 
a  significant  departure  from  "common  knowledge"  in  the  medical 
community,  which  has  long  contended  that  non-interacting 
magnetic  nanoparticles  are  the  ideal  material  for  hyperthermia 
treatment.  This  interaction  strength  may  be  controlled  by  varying 
interparticle  distance  using  different  coatings,  and  varying  the 
magnetic  strength  of  the  nanoparticles  by  nanoparticle  chemistry. 
This  result  is  expected  to  provide  improved  guidance  to  industry 
for  the  development  of  new  nanoparticles  for  hyperthermia 
cancer  treatment. 


Q(A") 


FIGURE  3:  SANS/USAIMS  data  showing  the  interaction  peak  in  the  double 
dextran  sample  with  a  concentration  of  25.6  mg/ml  at  Q.  4H10'^  A'' 
(radius  =  1 57  nm),  and  the  possible  interaction  peak  in  the  single  dextran 
layer  sample  with  a  concentration  of  21.8  mg/ml  at  0  .  9H10^  A  '' 
(radius  =  698  nm). 
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Protein  Folding  in  Membranes 


X.  Han\  K.  Hristova\  W.  C.  Wimley^,  M.  Mihailescu^ 


In  any  organism,  roughly  a  quarter  of  all  proteins  assemble 
into  their  biologically  active  conformation  only  when 
they  are  embedded  in  a  lipid  bilayer  membrane  [1].  The 
basic  physical  principles  that  drive  membrane  protein 
folding  have  been  known  for  some  time;  the  bilayer  is  a 
two  dimensional  fluid  with  a  symmetric  normal  profile 
characterized  by  a  hydrophobic  (water-repelling)  core  and 
two  hydrophilic  bilayer-water  interfaces.  However  the 
influence  of  transmembrane  anisotropy,  dynamic  gradients 
and  chemical  heterogeneity  on  protein  folding  remains 
poorly  explored.  Furthermore,  the  degree  to  which  proteins 
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FIGURE  1:  A  lipid  film  containing  1  peptide  per  25  lipids  was  deposited  on  a  glass  substrate 
and  hydrated  through  the  vapor  phase  to  generate  oriented  fluid  multibilayers,  left.  The 
oriented  sample  was  placed  in  the  AND/R  beam,  and  Bragg  diffraction  peaks  were  collected  in 
a  6  -26  scan,  bottom.  The  bilayer  scattering  length  profiles,  upper  right,  are  generated  from 
the  diffraction  peaks  using  a  Fourier  transform.  Isomorphous  replacement  of  lipids  or  peptides 
in  the  bilayer  with  specifically  deuterated  variants  allows  for  the  localization  of  the  labels  in  the 
bilayer  using  difference  profiles,  top.  Bilayer  image  courtesy  of  Dr.  Scott  Feller. 


and  bilayers  must/ can  modulate  their  hydrophobicity 
profiles  to  achieve  hydrophobic  complementarity  is  not  well 


understood.  (Hydrophobic  complementarity  is  a  measure  of  how 
well  hydrophobic  regions  of  one  molecular  segment  line  up  with 
hydrophobic  regions  on  another  molecular  segment).  Recently, 
neutron  diffraction  was  used  to  address  these  questions  in  studies 
utilizing  peptide  model  systems  related  to  several  classes  of 
membrane  proteins  [2,3]. The  results  reveal  the  surprising  degree 
to  which  hydrophobic  complementarity  can  be  compromised  by 
protein-protein  interactions. 

The  studies  described  here  were  the  result  of  a  long-standing 
common  interest  in  protein  folding  in  membranes  between 
the  researchers  at  Johns  Hopkins  and 
the  Tulane  group  which  has  developed 
peptide  model  systems  for  studying  protein 
folding  in  membranes.  One  of  these  model 
systems,  designed  and  studied  at  the 
Tulane  lab,  is  a  hydrophobic  hexapeptide, 
Acetyl-Tryptophan-(Leucine)j  (AcWL^)  that 
spontaneously  and  reversibly  self-assembles 
in  membranes  into  highly  ordered  (i-sheet 
oligomers  [4] .  7\lthough  the  binding  and 
energetics  of  this  membrane  protein  model 
are  well  understood,  structural  details  of 
the  membrane-embedded  ^-sheets  were 
lacking  until  a  study  of  deuterium  labeled 
AcWL^  oligomers  in  oriented  bilayers  using 
the  Advanced  Neutron  Diffractometer/ 
Reflectometer  (AND/R)  at  the  NCNR  [3]. 

Previous  studies  of  the  assembly  of  AcWL^ 
into  (E-sheets  in  membranes  have  provided 
insights  into  membrane  protein  folding.  Yet, 
the  exact  structure  of  the  oligomer  in  the  lipid 
bilayer  is  still  unknown.  Based  on  numerous 
biophysical  experiments,  a  hypothetical 
model  was  proposed  in  which  AcWL^  forms 
transmembrane  antiparallel  |B-sheets  that  are 
centered  in  the  hydrocarbon  core  in  the  bilayer 
[5].  While  hydrogen  bonding  is  known  to  be 
the  main  driving  force  for  folding  of  AcWL,  in  membranes,  this 
model  also  contains  the  implicit  assumption  that  hydrophobic 
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FIGURE  2;  Transbilayer  scattering 
length  profiles  of  fluid  bilayers 
containing  either  the  monomeric, 
interfacial  peptide  AcWL^,  left,  or 
the  oligomeric,  transmembrane 
P-sheet  peptide  AcWL^,  right. 
Experiments  were  done  with  1 
peptide  per  25  lipids.  Isomorphous 
replacement  was  achieved  by 
substituting  protonated  peptide 
with  peptide  containing  two 
consecutive  deuterated  leucine 
residues.  The  difference  profiles, 
top,  show  the  average  location 
of  deuterium  labeled  leucines 
in  the  bilayer  +/-  experimental 
uncertainty.  While  AcWL^  shows 
the  expected  interfacial  localization, 
the  oligomeric  AcWL^  deuterons 
are  present  all  across  the 
membrane.  Thus  the  transbilayer 
disposition  of  the  peptide  p-sheets 
must  be  randomly  staggered  rather 
than  regular  and  symmetric  as  in 
existing  models,  bottom. 


complementarity  drives  a  regular,  symmetric  transbilayer 
disposition  of  the  peptide  (Fig.  2,  bottom  left.)  This  model  was 
tested  using  neutron  diffraction  of  fluid  lipid  bilayers  containing 
peptides  (Rg.  1).  Isomorphous  replacement  of  peptides  with 
variants  that  were  selectively  deuterated  on  two  consecutive 
leucine  residues,  either  the  L2L3  leucines  or  the  L5L6  leucines, 
was  used  to  study  the  disposition  of  the  peptides  in  bilayers  and 
compare  the  measured  distribution  with  that  predicted  by  the 
model. 


membrane  proteins  in  their  native  membrane  environment  and 
highlight  the  importance  of  the  interplay  between  hydrophobic 
complementarity  and  protein-protein  interactions  in  determining 
the  structure  of  membrane  proteins. 
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Unexpectedly,  the  pairs  of  adjacent  deuterium-labeled  leucines  in 
AcWL.  have  no  well-defined  peak  or  dip  in  the  transmembrane 
distribution  profiles,  but  instead  are  distributed  all  across 
the  membrane  (Fig.  2).  The  result  was  the  same  for  L2L3 
labeled  peptides  and  for  L5L6  labeled  peptides.  This  result 
indicated  that  there  is  heterogeneity  in  the  depth  of  membrane 
insertion,  implying  that  the  [^-strands  are  staggered  along  the 
transmembrane  direction.  At  the  same  time,  the  monomeric 
homolog,  AcWL^  exhibits  a  homogeneous,  well-defined, 
interfacial  location  in  neutron  diffraction  experiments  and 
is  strongly  excluded  from  the  hydrocarbon  core.  Thus,  while 
the  bilayer  location  of  monomeric  AcWL^  is  determined  by 
hydrophobicity  matching,  or  complementarity,  within  the  bilayer, 
the  AcWLj  molecules  in  the  oligomers  are  positioned  at  different 
depths  within  the  bilayer  because  they  assemble  into  a  staggered 
transmembrane  ji-sheet.  The  details  of  AcWL^  assembly  are 
thus  dominated  by  protein-protein  interactions  rather  than  by 
hydrophobic  complementarity  [3]. These  results  have  implications 
for  understanding  and  predicting  the  structure  and  folding  of 
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The  type  1  Human  Immunodeficiency  Virus  (HIV-1)  is 
responsible  for  the  AIDS  pandemic  and  has  infected 
over  30  million  people  worldwide.  Effective  public  health 
policy  can  play  an  important  role  in  drastically  reducing 
the  spread  of  AIDS.  However,  a  rise  in  the  number  of 
patients  infected  by  drug  resistant  strains  of  HIV  has  made 
treatment  more  difficult.  Patients  with  drug-resistant 
strains  are  often  out  of  options  and  have  greatly  reduced  life 
expectancies.  Although  treatments  often  involve  a  "cocktail" 
of  drugs  targeting  different  stages  of  the  viral  life-cycle 
simultaneously,  currently  there  are  no  therapeutic  agents  that 
target  the  assembly  stage. 

Formation  ot  HIV-1  is  mediated  by  the  viral  Gag  polyprotein. 
Expressed  in  the  cellular  cytoplasm,  Gag  eventually  targets 
the  inner  surface  of  the  cellular  membrane  of  the  infected  host 
cell  where  viral  assembly  occurs.  Molecular  insight  from  early 
cryo-electron  microscopy  data  showed  Gag  in  the  immature 
spherical  virus  as  elongated  rods  radiating  from  the  membrane 
with  one  end  tightly  bound  to  the  viral  genome  [1].  However, 
in  a  recent  study  using  small  angle  neutron  scattering  as 
well  as  other  techniques,  it  was  found  that  the  properties  of 
monomeric  Gag  in  solution  are  incompatible  with  an  extended 
structure  [2].  Rather,  Gag  likely  exists  in  several  compact 
conformations  in  solution,  most  likely  due  to  the  presence 
of  several  unstructured,  flexible  domains  in  the  protein. 
These  results  imply  that  the  protein  must  undergo  a  large 
conformational  change  when  it  assembles  into  a  vims  particle. 
Understanding  the  mechanism  of  this  conformational  change 
would  give  important  insights  into  retroviral  assembly. 

Our  work  has  elucidated  some  of  the  factors  that  contribute 
to  conformational  changes  of  the  HIV-1  Gag  protein  from  a 
compact  structure  in  solution  to  an  extended  structure  in  the 
immature  virus.  The  size  of  Gag  bound  to  the  bilayer  interface 
was  studied  in  a  well-defined  in  vitro  system.  The  bio-mimetic 
environment  for  observing  Gag  association  consisted  of  a 
supported  membrane  attached  to  a  gold  surface  via  a  polymer 
tether.  This  'sparsely  tethered' bilayer  lipid  membrane  (tBLM) 
provided  a  highly  fluid  lipid  system  decoupled  from  the  solid 
substrate  that  supports  it  [3].  Membranes  were  prepared  by 


the  rapid  solvent  exchange  method  [4].  Structural  characterization 
of  the  bilayer:protein  system  was  done  by  neutron  reflectivity  (NR) 
on  the  AND/R  instrument  at  the  NCNR  [5]. 

Viruses  acquire  their  membranes  by  budding  from  the  host  cell 
membrane.  Lipid  extracts  of  purified  HIV-1  particles  identified 
an  abundance  of  lipids  present  only  at  small  concentrations  in 
uninfected  cells  [6] .  Among  these  are  included  negatively  charged, 
anionic  lipids  that  are  thought  to  associate  with  the  n-terminal 
domain  of  the  Gag  protein  (known  as  the  Matrix  or  MA  domain). 
Our  experiments  involve  a  simplified  lipid  mixture  that  stUl 
captures  the  negative  surface  charge  density  of  the  native  viral 
membrane. 

To  test  if  the  membrane  targeting  domain  of  Gag  would  associate 
with  our  model  system,  we  first  used  NR  to  examine  just  the  MA 
construct.  Figure  la  shows  the  reflectivity  curves  of  the  neat  bilayer 
and  the  membrane  exposed  to  two  different  concentrations  of  the 
MA  protein.  Differences  between  the  pure  membrane  system  and 
the  membrane  with  MA  protein  are  clearly  seen  at  the  bottom 
of  the  panel.  In  Fig.  lb  the  reflectivity  data  were  fitted  using  a 
simple  box  model  to  represent  the  Au  substrate,  tBLM  and  the 
MA  domain.  In  this  model,  the  MA  protein  is  represented  by  a 
single  slab  of  fixed  neutron  scattering  length  density  (nSLD).The 
data  indicates  that  the  MA  domain  of  Gag  extends  40  A  from 
the  membrane  surface  and  that  this  extension  is  invariant  with 
concentration. 
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FIGURE  1 :  {a)  Fresnel  normalized  neutron  reflectivity  profile  of  the  pure 
tethered  bilayer  {black  circles),  bilayer  exposed  to  1  mmol/L  of  Gag  MA 
protein  (red  circles),  and  to  10  mmol/L  of  Gag  MA  (blue  circles).  Fits  to 
the  experimental  data  are  shown  as  solid  lines.  The  differences  in  the 
normalized  reflectivity  profile  between  pure  bilayer  and  bilayer  with  MA  are 
shown  at  the  bottom  of  the  graph,  (b)  Neutron  SLD  profile  that  best  fits  the 
reflectivity  data. 
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Because  MA  is  a  rigid  domain,  a  fixed  dimension  in  our  NR 
modeling  suggests  that  the  protein  binds  in  a  specific  orientation 
with  regard  to  the  bilayer  surface,  in  agreement  with  NMR  and 
x-ray  measurements  [7].  Using  an  atomic  resolution  structure  of 
the  MA  domain,  a  refined  nSLD  profile  was  used  in  place  of  the 
single  box  representation.  Figure  2a  shows  the  volumetric  and 
nSLD  profile  for  MA  using  the  suggested  binding  orientation 
from  previous  studies.  The  resulting  best  fit  model  is  shown 
in  Fig.  2b  where  only  the  volume  fraction  and  position  in  z  of 
the  protein  layer  are  free  parameters.  The  high  quality  of  fit  to 
the  NR  data  strongly  suggests  that  the  orientation  of  MA  on 
the  membrane  is  correct.  These  experiments  also  show  that 
the  volume  fraction  of  protein,  which  can  be  related  to  surface 
coverage,  doubles  from  15  %  to  30  %.  This  doubling  happens 
while  increasing  the  concentration  of  MA  in  solution  10-fold, 
hence  suggesting  the  surface  coverage  of  MA  is  saturated.  Finally, 
it  seems  that  MA  insertion  into  the  lipid  layer  is  slight 
(<  1  A)  v^dth  minimal  perturbation  to  the  membrane  structure. 
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FIGURE  2:  [a)  Matrix  (MA)  domain  oriented  with  membrane  binding  region 
on  the  left  hand  side.  Volumetric  profile  and  nSLD  profiles  were  calculated 
using  2.5  A  slabs  of  MA  along  z  taking  into  account  exchangeable  hydrogens 
on  the  protein,  (b)  Best  fit  nSLD  profile  of  the  tBLM  system  with  no  MA 
protein  (black  line),  and  bound  MA  protein  at  1  mmol/L  (red  line]  and  10 
mmol/L  (orange  line)  concentrations. 

Having  verified  MA  binding,  measurements  were  made  using  the 
fuU  length  Gag  protein.  A  series  of  measurements  were  made 
under  different  conditions  thought  to  affect  Gag  conformation. 
1)  Purified  Gag  protein  was  introduced  to  the  membrane  and 
allowed  to  associate.  2)  The  bound  protein  was  exposed  to  a  short 
DNA  strand  (TGx7)  known  to  bind  the  c-terminal  nucleocapsid 
(NC  )  domain  of  Gag  and  long  enough  to  span  approximately 
three  molecules.  The  NC  domain  is  normally  responsible  for 
associating  with  the  viral  genome.  3)  A  high  salt  buffer  rinse  of 
500  mmol/L  NaCl  was  flushed  through  the  system,  potentially 
screening  the  electrostatically  driven  interactions  between  the 
protein  and  DNA. 

Due  to  the  lack  of  a  well-defined  structure,  owing  to  the  protein's 
flexible  segments,  modeling  the  full-length  protein  provided 
a  difficult  challenge.  A  series  of  boxes  extending  240  A  from 
the  membrane  surface  was  used.  Because  of  the  large  number 
of  parameters  introduced  by  our  model,  a  Monte-Carlo  error 
analysis  was  performed  to  determine  the  statistical  variability 
in  the  nSLD  profiles.  Figure  3a  shows  the  nSLD  profile  of  the 
different  measurements  where  line  thickness  represents  95% 


confidence  intervals  in  our  fits.  The  neat  bilayer  is  seen  as  the 
black  solid  line.  Full-length  Gag  protein  bound  to  our  bio- 
mimetic  membrane  system,  in  red,  is  best  modeled  by  an  80  A 
structure  from  the  lipid  surface.  Upon  exposure  to  TGx7  the 
protein  layer  dramatically  extends  out  to  200  A.  Finally,  upon 
rinsing  with  the  high  salt  buffer  the  original  Gag  nSLD  profile  is 
nearly  recovered. 
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FIGURE  3:  (a)  nSLD  profile  of  full-length  HIV-I  Gag  protein  on  a  tBLM.  Neat 
lipid  bilayer  (black),  bound  Gag  protein  (red).  Gag  +  TGx7  DNA  strand  (blue), 
Gag:TGx7  500  mmol/L  NaCI  salt  rinse  (green).  The  inset  cartoon  illustrates 
how  the  charged  ends  of  the  Gag  cause  it  to  fold  toward  the  surface,  and 
then  how  the  viral  strands  attach  to  NC,  extending  and  crosslinking  the  Gag 
molecules,  (b)  Illustrative  models  of  folded  and  extended  conformations  of 
Gag  on  a  membrane  surface. 

This  extension  of  Gag  in  the  presence  of  both  the  membrane 
and  the  TGx7  nucleic  acid  strand  provides  a  model  for 
conformational  changes  in  Gag.  This  model,  illustrated  in  the 
inset  of  Fig.  3a,  shows  the  MA  and  the  NC  domains  associating 
with  the  negatively  charged  surface  of  the  lipid  membrane.  The 
NC  domain  also  contains  several  basic  residues  which  aid  in  its 
association  with  the  viral  RNA.  However  in  the  absence  of  RNA 
and  given  the  flexible  regions  of  the  protein,  the  NC  domain 
is  free  to  bend  down  and  associate  with  the  membrane.  Only 
when  nucleic  acid  of  significant  binding  strength  is  present  does 
the  NC  domain  disassociate  with  the  membrane  and  bind  to 
the  RNA  strand.  This  causes  extension  of  the  protein  molecule 
in  part  due  to  crosslinking  several  Gag  molecules  together.  This 
work  provides  highly  suggestive  insight  into  the  mechanism  by 
which  extension  occurs.  Our  model  unifies  earlier  in  vivo  and 
solution  studies  [1,2].  An  understanding  of  virus  assembly  may 
lead  to  the  development  of  therapeutics  that  inhibit  proper  virus 
formation. 
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Quantum  chromodynamics,  which  describes  the  strong 
interaction  between  quarks,  is  non-perturbative, 
making  rigorous  direct  calculations  at  low  energies 
intractable.  Instead,  complex,  multi-parameter  theoretical 
models  have  been  developed  to  tackle  nucleon-nucleon  (NN) 
interactions.  Combinations  of  computational  and  theoretical 
advances  in  recent  years  have  led  theoretical  models  to 
predict  properties  of  few  nucleon  systems  with  an  uncertainty 
of  less  than  one  percent.  In  systems  larger  than  two 
nucleons  less  phenomenological  and  more  poorly  understood 
three  nucleon  (3NI)  interactions  must  be  included  along 
with  NN  models.  Neutron  scattering  lengths,  which  describe 
a  neutron's  5-wave  interaction  with  a  target  nucleus,  are 
predicted  by  NN+3NI  models,  and  therefore  provide  crucial 
benchmarks  in  the  testing  of  various  theoretical  approaches. 
Scattering  lengths  also  play  an  important  role  in  effective 
field  theories  since  they  use  low-energy  observables  to 
constrain  mean-field  behavior.  In  this  highlight  our 
measurement  of  the  spin-dependent,  incoherent 
neutron-helium-3  (n-'He)  scattering  length,  b.,  is  described. 

The  NIST  neutron  interferometer  and  optics  facility  (NIOF) 
has  provided  measurements  of  neutron  scattering  lengths 
for  neutron-hydrogen  and  neutron-deuterium,  and  the 
spin-independent  scattering  length  for  n-^He  with  relative 
uncertainties  below  one  percent  [1].  This  is  only  the  second 
experimental  attempt  to  directly  measure  the  spin-dependent 
scattering  length  in  'He.  The  previous  measurement 
[2]  was  done  using  a  spin-echo  apparatus  at  the  Institut 


Laue-Langevin  which  achieved  a  result  inconsistent  with  any 
theoretical  prediction.  Using  an  interferometer  to  measure 
b,  introduces  a  completely  separate  set  of  systematics.  This 
experiment  marks  the  first  time  a  polarized  gas  sample  has  been 
used  in  a  neutron  interferometer  experimicnt. 

Upstream  of  the  interferometer  a  CoFe  supermirror  was  used  to 
polarize  the  neutrons  to  nearly  93  %  polarization.  The  neutron's 
spin  could  be  flipped  180°  after  the  supermirror  with  almost  100 
%  accuracy  using  a  pair  of  coils  tuned  to  the  neutron's  Larmor 
frequency.  After  the  supermirror,  the  neutron's  polarization  was 
preserved  by  a  guide  field  made  up  from  a  series  of  permanent 
magnets.  To  confirm  that  the  neutron  polarization  (PJ  was  not 
changing  during  the  experiment,  the  polarization  was  measured 
periodically.  This  was  done  by  replacing  the  interferometer  with 
an  optically  thick  'He  cell  which  could  provide  analyzing  power 
up  to  99  %.  Two  different  techniques  that  were  used  to  measure 
P^^  and  the  spin  flipper  efficiency  (s)  are  shown  in  Fig.  1. 

The  NIST  glass  shop  [3]  fabricated  four  boron-free  target  cells  for 
use  in  this  experiment.  Each  cyhndrical  cell  had  outer  dimensions 
0  25.4  mm  x  42  mm  and  was  sealed  with  approximately  1.5  bar 
of  ^He  gas.  The  cells  were  made  optically  thin  to  allow  for  some 
neutron  transmission  even  during  strong  absorption  when  the 
'He  and  neutron  spins  were  aligned  anti-paraUel.  The  'He  gas 
was  polarized  to  approximately  65  %  using  spin-exchange  optical 
pumping  at  a  separate  facility.  This  procedure  eliminated  the  added 
complexity  and  heat  loads  to  the  interferometer  setup  and  provided 
the  experiment  with  viable  target  samples.  Helmholtz  coils  placed 
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FIGURE  1 :  Polarization  analysis  using 
two  different  techniques,  (a)  The  results 
from  the  Asymmetry  Method  for  both  the 
neutron  polarization,  P^,  and  the  spin  flipper 
efficiency,  s  in  percent  versus  measurement 
instance,  (b)  The  polarization  results  using 
the  Normalized  Transmission  Method. 
Here  the  spin  flipper  efficiency  is  not 
directly  measured  but  instead  the  product 
(2s-1)*P  is.  In  both  cases  drawn  lines  are 
fit  to  the  data.  Both  techniques  agree  with 
each  other  with  a  relative  uncertainty  of 
less  than  0.001. 
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around  the  interferometer  provided  a  uniform  magnetic  field  of 
15xlO""'T  to  limit  the  loss  of  helium  polarization  due  to  magnetic 
gradients.  Cell  lifetimes  in  the  interferometer  were  up  to  175  h. 

The  interferometer  crystal  is  isolated  from  vibrations  and 
temperature  controlled  to  ±  5  K  which  enables  exceptional  phase 
stability  and  contrast  as  high  as  80  %.  The  skew  symmetric 
interferometer  consists  of  3  silicon  blades  on  a  common  base.  The 
first  blade  uses  Bragg  reflection  to  coherently  split  the  neutron's 
wavefunction  into  two  spatially  separated  paths,  I  and  II.  When  a 
neutron  passes  through  a  ^He  cell  in  Path  I  the  nuclear  potential 
of  the  helium  atoms  causes  a  change  in  phase  (zij)  of  the  neutron's 
wavefunction.  Beams  along  the  two  paths  are  brought  together 
in  the  final  blade  where  they  interfere  with  one  another.  Relative 
differences  in  the  neutron's  phase  {^Xr  ^Xn>  ^orig  each  path 
modify  the  count  rates  at  two  detectors,  labeled  O-  and  H-,  located 
behind  the  interferometer. 

A  thin  quartz  sample  called  a  phase  flag  is  rotated  to  controUably 
vary  the  phase  along  the  beam  paths  to  construct  an  interferogram 
(see  Fig.  2c)  that  is  used  to  determine  the  phase  of  the  target  cell. 
The  spin-dependent  scattering  length  is  proportional  to  the  phase 
difference  between  when  the  neutron  and  ^He  spins  are  parallel  to 
when  they  are  anti-parallel.  Boron-free  glass  was  placed  in  Path  II  in 
order  to  compensate  for  the  phase  shift  caused  by  the  cell  windows. 


FIGURE  2:  Experimental  setup,  (a)  Neutrons  enter  from  the  right  of  the 
picture.  One  can  see  the  compensation  glass,  the  target  cell,  and  the  quartz 
phase  flag,  (b)  One  of  the  target  cells  used.  Solid  rubidium  which  is  heated 
and  used  in  Spin-Exchange  Optical  Pumping  can  be  seen  in  the  picture 
as  dark  spots  along  the  cell  walls,  (c)  A  typical  interferogram.  The  blue 
=1     curve  (lines  are  fits)  is  the  0-beam  intensity  for  when  the  neutron  and  ^He 
spins  are  aligned  parallel.  The  red  curve  is  for  when  the  neutron  spin  has 
undergone  a  spin  flip.  Lower  intensity  In  the  red  curve  Is  due  to  increased 
neutron  absorption  in  that  case.  A  phase  shift  around  30°  can  be  seen  by 
comparing  the  two  curves. 


The  target  cell's  polarization  was  monitored  throughout  the 
experiment.  A  third  detector  labeled  C4  was  placed  after  the 
interferometer  and  directly  behind  the  -"He  cell.  C4  measured 
the  transmission  of  neutrons  through  the  ^He  for  both  neutron 
spin  states.  A  ratio  of  the  count  rates  for  both  cases  gave  the 
polarization  of  the  helium  gas. 

Our  (preliminary)  result,  b.  =  -2.410  (19)  fm,  is  in  good  agreement 
with  the  previous  measurement  (Fig.  3)  of  the  spin-dependent 
^He  scattering  length  [2]. The  present  result  and  the  previous  one 
are  systematically  limited  by  the  small  but  nonzero  spin-paraUel 
absorption  ((j  )  known  only  to  one  percent.  Better  experimental 
determination  of    is  needed  and  would  greatly  improve  the 
current  error  bars.  Known  NN+3NI  models  do  not  match  our 
results,  showing  the  need  for  improved  theoretical  work.  Four 
nucleon  interactions  are  not  included  into  the  models  due  to 
the  difficulty  in  handling  long-range  coulomb  forces,  but  should 
constitute  a  tiny  correction  to  any  NN+3NI  predictions.  This 
measurement  is  part  of  the  ongoing  exploration  into  few  body 
systems  by  the  NIOF. 
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FIGURES:  Experiential  results  found  in  [1],  [2],  and  [4]  compared  with 
predicted  values  of  some  N!\1+3NI  theoretical  models  [5].  The  coherent 
scattering  length  measured  in  Refs.  1  and  4  is  given  by  b^  =  %  b^  +  b^. 
The  results  are  shown  as  bands  representing  the  result  ±  one  standard 
deviation,  o. 
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Neutron  Diffraction  Measurement  of  the 
Load  Response  of  Concrete 
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Concrete  based  on  portland  cement  is  the  most  widely 
used  building  material  in  the  world.  Due  to  the 
importance  and  long  history  of  cement-based  materials,  a 
large  body  of  research  has  been  assembled  about  behavior 
and  properties  relevant  to  civil  engineering  applications.  This 
research  includes  the  macroscopic  response  to  mechanical 
and  thermal  loads  and  the  interaction  with  chemical 
agents.  However,  concrete  is  a  material  of  considerable 
complexity  consisting  of  several  phases  both  amorphous 
and  crystalline,  and  it  has  structural  heterogeneities  such 
as  capillary  pores  and  micro-cracks.  Currently,  there  is  no 
comprehensive  theory  that  explains  the  mechanical  and 
chemo-mechanical  responses  associated  with  cement- 
based  concretes.  The  advances  in  the  simulation  of  concrete 
behavior  so  far  are  hampered  by  insufficient  experimental 
data  for  model  calibration  and  validation.  This  work  is  part 
of  an  extensive  effort  to  develop  experimental  strategies 
aimed  at  gaining  insight  into  how  micro-  and  macro- 
scale  phenomena  are  Hnked  in  the  mechanical  responses 
of  cement-based  materials  under  mechanical  loading  and 
drying  environments. 

The  vast  majority  of  cement  and  concrete  experimental 
research  is  based  on  macro-level,  bulk  responses  to  stressors 
by  means  of  strain  gauges  or  other  similar  techniques  which 
do  not  observe  or  explain  the  internal  micromechanics  of 
the  system  which  is  where  failure  actually  originates  [1].  In 
order  to  observe  mechanical  responses  on  the  micro-  and 
macro-scale,  probing  techniques  need  to  be  developed  for 
obtaining  quantitative  measurements  within  these  complex 
composites.  There  is  much  work  in  the  fields  of  neutron  and 
x-ray  diffraction  that  demonstrates  that  these  techniques  can 
be  used  as  phase  resolving,  non-destructive  probes  [2,3].  The 
novelty  of  these  methods  is  that  diffraction  provides  a  means 
by  which  the  strains  within  native  crystalline  components 
of  the  cement  and  aggregate  can  be  used  as  internal  strain 
gauges,  thereby  providing  detailed  information  concerning 
the  load  transference  mechanism(s)  within  and  between 
individual  components  of  the  composite.  Since  the  major 
hydrated  portland  cement  reaction  product  (calcium  silicate 


hydrate)  has  an  amorphous  structure,  suitable  coherently  scattering 
crystalline  components  are  limited  to  the  minority  constituents 
calcium  hydroxide  (CH)  and  un-reacted  cement  phases  such  as 
tri-calcium  silicate  (C^S). 

In  recent  synchrotron  experiments  lattice  strains  were  measured 
within  the  native  CH  phase  of  hydrated  portland  cement.  Strains 
were  induced  through  compressive  uniaxial  loading  of  cement 
samples.  Proof  of  concept  was  demonstrated  with  a  strain  accuracy 
1x10''  and  sufficient  strain  sensitivity  for  resolving  the  anisotropic 
response  of  different  crystaUographic  directions  in  Ca-iOH)^ 
grains.  However,  the  scope  of  these  results  was  limited  by  the  low 
penetration  (<  50  \j.m)  of  x-rays  which  presents  difficulties  for 
comparing  surface  behavior  -  which  is  affected  by  wetting  and 
drying  effects  -  and  bulk  behavior  that  is  susceptible  to  the  same 
effects  but  on  a  much  longer  time  scale. 

In  the  research  presented  here,  neutrons  were  utilized  in  order  to 
investigate  bulk  behavior  for  comparison  with  synchrotron  x-ray 
results  as  part  of  the  multi-scale  hierarchy  of  novel  experimental 
results  needed  to  gain  new  insight  to  the  broader  spectrum  of 
micro  and  macro-scale  mechanical  responses  in  hydrated  portland 
cement  when  subjected  to  mechanical  loads.  7\lso,  due  to  the 
penetration  of  neutrons  in  the  range  of  10"^  m  (x-ray:  10'''  m) 
neutron  diffraction  provides  improved  volume  average  and  grain 
statistics.  In  order  to  improve  the  substantial  incoherent  scattering 
background  from  hydrogen,  aU  cement  paste  specimens  were 
cured  for  at  least  one  year,  then  machined  and  stored  in  deuterium 
oxide  D^O. 

Two  types  of  experiments  were  performed:  (1)  the  strain  response 
resulting  from  uniaxial,  compressive  mechanical  loading;  and  (2) 
the  strain  response  due  to  shrinkage  induced  stresses.  The  second 
was  needed  since  the  samples  were  kept  in  the  Ca(OD)2-saturated 
D^O  until  the  time  that  the  diffraction  experiment  began.  Thus, 
the  effects  of  drying  became  relevant  since  this  phenomenon  was 
occurring  simultaneously  during  aU  experiments.  The  measured 
strains  are  lattice  strains  that  characterize  the  response  of 
crystalline  Ca(OD),.  In  order  to  separate  the  strain  responses  from 
simultaneous  drying  and  loading,  the  experiment  was  performed 
as  a  series  of  loading-unloading  cycles  as  shown  in  Fig.  1. 
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FIGURE  1:  Compressive  tests  of  cement  paste  prisms  [10  mm  x10  mm 
x15  mm)  in  loading-unloading  cycles.  The  strain  follows  the  applied  load  in 
linear  dependence.  However,  with  increasing  time  and  water  evaporation 
the  strains  in  the  unloaded  states  (Q.43  MPa)  become  more  compressive. 
A  subsequent  drying  measurement  of  an  unloaded  sample  during  water 
evaporation  revealed  that  the  strains  in  the  unloaded  states  exhibit  a  similar 
rate  as  the  shrinkage  strains  originating  from  dehydration/evaporation  only 
(Fig.2]. 
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FIGURE  2:  Comparison  of  strains  from  a  drying  test  with  the  strains  of 
the  unloaded  states.  The  measured  CatOO)^  lattice  strains  characterize 
the  behavior  of  a  minority  constituent;  the  behavior  of  the  cement  paste 
aggregate  as  a  whole  is  better  captured  with  surface  attached  strain  gages. 
The  comparison  of  suitably  normalized  lattice  strains  and  macroscopic 
strains  is  shown  in  Fig.  3. 
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FIGURE  3:  Comparison  between  macro-strain  and  lattice  strain 
measurements  for  in  both  the  wetting  and  drying  environment.  The  lattice 
strains  from  wetting  are  for  a  dry  sample  being  equilibrated  in  room  air. 


The  rates  of  change  differ  substantially  between  macro-  and 
lattice  strains.  Lattice  strains  represent  elastic  changes  and  as 
such  they  are  a  fraction  of  the  macroscopic  strain  change  which 
consists  of,  among  other  contributions,  the  deformation  of  the 
pore  and  capillary  structure  and  the  opening  and  closing  of 
micro-cracks  driven  by  hydration  of  hygroscopic  constituents.  A 
notable  difference  between  lattice  strains  in  wetting  and  drying  is 
that  the  drying  strains  are  decreasing  and  they  qualitatively  follow 
the  macroscopic  strains.  The  wetting  strains  are  increasing  but 
see  a  drop  after  10  h,  indicative  of  the  very  low  tensile  strength 
(compared  to  its  compressive  strength)  of  cement  paste  and  the 
onset  of  micro-crack  propagation  and  inhomogeneous  stress 
relief.  The  macro-strain  data  in  Fig.  3  represent  surface  strains, 
and  their  changes  have  different  consequences  for  the  stresses 
present  in  the  bulk  (where  the  neutron  measurements  are  done). 
In  wetting,  the  near  surface  regions  experience  swelling  first 
which  creates  a  tensile  stress  in  the  interior,  some  of  which  is 
reflected  in  the  lattice  strains  from  wetting.  As  the  hydration/ 
wetting  gradient  progresses  deeper  into  the  bulk,  the  macro- 
strains  at  the  surface  begin  to  saturate  and,  with  the  swelling  now 
uniformly  present  in  the  bulk,  the  bulk  tensile  stresses  vanish.  This 
accounts  for  the  drop  in  the  lattice  strains  from  wetting  several 
hours  later.  The  rate  of  change  is  much  lower  for  the  drying 
strains,  hence  lattice  and  macro-strains  are  qualitatively  similar. 

Overall,  results  obtained  from  neutron  diffraction  contribute 
to  the  growing  base  of  experimental  data  that  measure  the 
mechanical  response  of  single  constituents  in  concrete.  Such  data 
is  essential  for  improved  understanding  of  the  material  and  for 
advanced  simulation  [4]  of  this  complex  heterogeneous  system. 
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Off-Shore  Oil  Pipelines:  Residual  Stresses  in  Girth  Welds 
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Offshore  oil  and  gas  extraction  plays  a  vital  role  for  the 
global  energy  supply.  It  is  expected  that  much  future 
development  and  exploration  of  major  reserves  will  take  place 
even  farther  offshore  and  in  greater  depths  than  today.  With 
extraction  and  production  systems  installed  at  depths  of 
2000  m  to  4000  m,  riser  pipeline  systems  that  are  used 
to  transfer  oil  or  gas  from  the  wells  on  the  seabed  to  the 
production  structure  have  to  withstand  motion  due  to 
floating  vessels  as  well  as  vortex  induced  vibrations  in  the 
marine  environment,  subjecting  them  to  the  axial  bending/ 
unbending  cycles  that  give  rise  to  accumulating  local  plastic 
deformation  (fatigue)  [1].  Because  of  residual  stresses 
and  geometric  stress  concentrations,  the  damaging  effects 
are  concentrated  in  the  regions  of  circumferential  girth 
welds  that  connect  individual  pipe  sections.  Therefore, 
the  determination  of  these  residual  stresses  is  essential 
for  understanding  and  estimating  fatigue  life  and  for 
maintaining  pipeline  integrity. 

This  highlight  shows  our  results  from  neutron  diffraction 
measurements  on  two  welds,  one  with  a  ferritic  weld  metal 
and  one  with  a  fully  austenitic  weld  where  different  stress 
profiles  are  formed  due  to  different  cooling  conditions,  phase 
transformation  and  deformation  behavior  upon  cooling  of 
the  weld  metal.  The  two  weld  metals  were  chosen  based  on 
preliminary  fatigue  testing  of  full-scale  pipes  showing  that 
girth  welds  with  nickel-based  weld  metal  exhibit  improved 
fatigue  performance  [2]. 


Failure  due  to  fatigue  has  not  been  of  much  importance  in  design 
and  safety  analysis  of  land-based  pipeline  systems  because  they  are 
mostly  subjected  to  static  loading.  However,  many  new  oil  and 
gas  reserves  are  found  offshore  at  increasing  depths  and  distance 
from  land  with  fatigue  in  riser  pipelines  arising  from  motion  as 
noted  above.  Major  factors  affecting  fatigue  life  include  the  yield 
strength  of  the  material  and,  with  respect  to  welds,  the  residual 
stress  state,  the  applied  stress,  and  the  number  of  loading/unloading 
cycles  to  which  they  are  subjected.  Interest  is  focused  on  the  girth 
welds  for  pipelines  in  general,  with  a  number  of  recent  neutron 
experiments  targeting  the  residual  stresses  in  the  girth  weld  region 
[3,4].  It  is  critical  to  characterize  the  residual  stress  state  of  welds 
in  order  to  develop  complete  mechanistic  understanding  of  fatigue 
performance  and  to  guide  the  fatigue  design  of  welded  joints.  We 
measured  residual  stresses  in  the  girth  weld  regions  of  two  sections 
of  X65  (minimum  yield  strength  451  MPa)  pipe:  one  with  ferritic 
weld  metal  and  another  with  a  fuUy  austenitic,  nickel-based  alloy 
weld  metal. 

Important  differences  between  the  two  weld  metals  are  in  the 
solidus-liquidus  ranges  and  in  the  microstructure.  The  temperature 
range  between  the  fuUy  liquid  state  and  the  fully  solidified  state  is 
between  1250  °C  and  1350  °C  for  the  austenitic  weld  metal  and 
between  1500  °C  and  1520  °C  for  the  ferritic  weld  metal.  The 
microstructural  differences  are  also  significant.  The  ferritic  welds 
exhibit  a  solid  state  phase  transition  of  y-Fe  (face  centered  cubic 
austenite)  to  a-Fe  (body  centered  cubic  ferrite)  at  ~  720°C.  There  is 
a  volume  expansion  associated  with  this  phase  transformation  that 
alters  the  residual  stress  generated. 
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FIGURE  1:  Neutron  diffraction 
geometry  for  a  large  pipe  section 
with  girth  weld  in  the  tangential 
measurement  direction.  The  right 
side  shows  the  cross  section 
of  the  fully  austenitic  girth  weld 
and  micrographs.  Note  the 
sharp  boundary  between  the 
crystallographically  different  base 
metal  and  the  weld  metal  (outlined 
by  the  dashed  line).  The  micrograph 
for  the  ferritic  weld  metal  is  not 
shown,  but  grains  sizes  are  of  the 
same  order  of  magnitude  as  the 
base  metal. 
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Also,  in  the  ferritic  welds  the  phase  transformation  refines 
the  as-sohdified  grain  structure  to  small,  often  needle  shaped 
crystaUites  with  typical  sizes  up  to  20  |am,  whereas  the  austenitic 
nickel-based  alloy  weld  metal  retains  the  as-solidified  columnar 
grain  structure  with  often  millimeter  size  grains  (lower  right. 
Fig.  l).'The  austenitic  grains  retain  the  face  centered  cubic  crystal 
structure  from  sohdification  to  room  temperature,  thus  allowing 
for  these  large  grains.  As  the  heat  transfer  conditions  change 
from  weld  pass  to  weld  pass,  grain  growth  direction  and  local 
texture  change  as  well,  which  makes  neutron  strain  measurements 
very  difficult  in  the  weld  region. 

A  fiandamental  underlying  requirement  of  neutron  diffraction 
strain  measurements  is  that  the  grains  contributing  to  a 
diffraction  peak  should  be  uniformly  distributed  within  the  gage 
volume.  With  gage  volume  sizes  of  3  mm  x  3  mm  x  3  mm,  in 
the  case  of  the  austenitic  weld  metal  with  large  grain  size  very 
few  or  no  grains  at  all  may  fulfill  the  diffraction  condition  at  a 
given  specimen  orientation.  Therefore,  the  grain  average  has  to  be 
improved  by  oscillating  the  sample  through  rotation  or  translation 
such  that  the  oscillation  does  not  introduce  additional  smoothing 
of  the  strain  gradient  to  be  measured.  This  requirement  represents 
a  considerable  challenge  when  confronted  with  a  pipe  section  of 
508  mm  diameter  and  25.4  mm  wall  thickness.  For  example,  the 
specimen  orientation  shown  in  Fig.  1  allows  a  possible  range  of 
oscillation  that  is  hmited  to  a  rotation  of  the  sample  table  (±  5  °) 
and  a  small  tangential  translation  (±  5  mm). 


Nickel  based  alloy  weld 


distance  from  weld  center  (mm) 

FIGURE  2:  Comparison  of  hoop  stresses  between  the  nickel-based  alloy  (top) 
and  the  ferritic  weld  (bottom).  Stress  levels  are  given  in  MPa.  The  dashed 
line  corresponds  to  the  girth  weld  region. 

The  principal  stresses  obtained  from  the  two  pipes  along  the  three 
principal  directions  (hoop,  axial,  radial)  show  the  effects  of  the 
aforementioned  differences  between  weld  metals.  Butt  welds  of 
the  above  shape  have  tensile  hoop  stresses  (Fig.  2)  that  peak  at 
the  level  of  the  yield  stress.  While  this  is  the  case  for  both  weld 
metals,  the  nickel-based  weld  with  the  lower  melting  range  shows 
considerably  less  hoop  stress  since  it  does  not  build  up  shrinkage 
stresses  to  the  magnitude  of  the  ferritic  weld  which  has  a  much 
higher  melting  range. 


However,  this  shrinkage  difference  is  not  reflected  to  the  same 
extent  in  the  axial  stresses  shown  in  Fig.  3.  The  near-surface 
tensile  stress  band  is  similar  for  both  welds  in  distribution  and 
magnitude.  Both  have  magnitudes  significantly  lower  than  the 
hoop  stress,  thus  offering  a  buffer  to  the  peak  tensile  stresses 
arising  from  bending.  However,  no  reduction  in  stress  levels 
comparable  to  the  hoop  direction  was  found  in  the  nickel  alloy 
weld  metal. 

These  are  among  the  first  measurements  of  residual  stresses  in 
full-scale  thick  section  pipes.  One  result  of  this  work  is  that  the 
improved  fatigue  life  in  nickel-based  weld  metal  compared  to 
ferritic  weld  metal  cannot  be  ascribed  to  a  reduction  in  axial 
stress.  It  remains  to  be  seen  whether  the  explanation  involves  the 
large  differences  in  hoop  stress.  Another  remaining  question  is 
the  role  of  microstructural  differences  in  resisting  fatigue  crack 
growth.  This  work  provides  high  quality  experimental  data  which 
wiU  aid  in  developing  and  calibrating  finite  element  models  for 
predicting  residual  stresses  in  weld  joints. 
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FIGURE  3:  Comparison  of  axial  stresses  between  the  nickel-based  alloy  (top) 
and  the  ferritic  weld  (bottom).  Stress  levels  are  given  in  MPa.  The  dashed 
line  corresponds  to  the  girth  weld  region. 
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Denser  than  Solid  Hydrogen:  Improving  Hydrogen  Storage 
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With  the  strong  dependence  on  other  countries  as 
sources  of  fossil  fuels,  the  high  prices  for  oil,  and 
concerns  about  environmental  impact  of  using  these  fuels, 
there  is  intense  interest  in  migrating  to  an  environmentally 
benign,  low  cost  alternative.  One  research  thrust  is 
investigating  the  use  of  clean-burning  hydrogen.  Specifically, 
one  of  the  key  engineering  challenges  to  building  a  clean, 
efficient,  hydrogen-powered  car  is  how  to  design  the  fuel 
tank.  Here  we  describe  a  candidate  storage  metal-organic- 
framework  (MOF)  material  that  displays  high  surface 
density  for  physisorbed  hydrogen. 

Storing  enough  hydrogen  to  allow  a  driving  range  equivalent 
to  a  gasoline  powered  car  requires  either  large  tanks  or  high 
pressures  of  gaseous  hydrogen.  While  there  are  promising 
mechanisms  for  storing  large  amounts  of  hydrogen  in 
solid  materials  that  take  up  much  less  space  than  the 
equivalent  amount  of  gas,  there  remain  several  roadblocks  to 
commercialization.  For  instance,  when  stored  on  high  surface 
area  materials,  hydrogen  is  weakly  bound  and  generally 
requires  low  temperatures  to  achieve  significant  storage 
capacities.  In  order  to  achieve  technologically  relevant  levels 
of  gravimetric  density,  the  density  of  adsorbed      must  be 
increased  beyond  levels  attained  for  typical  high  surface 
area  carbons.  In  this  highlight  we  report  our  work  showing 
a  strong  correlation  between  exposed  and  coordinatively 
unsaturated  metal  centers  (CUMCs)  and  enhanced  hydrogen 
surface  packing  density  in  many  framework  structures.  In 
particular,  we  show  that  the  MOF- 74  framework  structure 
with  open  Zn-*  sites  displays  the  highest  surface  density  for 
physisorbed  hydrogen  in  framework  structures. 

A  key  and  challenging  aspect  of  hydrogen  storage  research 
is  the  ability  to  control  the  H,  binding  energy  that  governs 
the  adsorption  in  materials  ranging  from  microporous  solids 
to  metal  hydrides.  In  this  regard,  a  very  attractive  property 
of  metal-organic  frameworks  (MOF),  comprised  of  metal 
ions  linked  via  organic  ligands  into  porous  three-dimensional 
solids,  is  that  their  pores  can  be  engineered  at  the  atomic 
scale,  thereby  enabling  some  control  over  the  H,  binding 
interaction.  Previously  we  have  shown  that  hydrogen  binds 


directly  to  the  unsaturated  metal  coordination  sites  within  certain 
classes  of  MOF  materials.  The  highest  observed  enthalpies  of 
hydrogen  adsorption  known  in  physisorption  systems  are  a  direct 
consequence  of  the  attraction  of  hydrogen  to  these  unsaturated, 
or  open,  metal  sites  [IJ.This  increase  in  adsorption  enthalpy  is 
necessary  to  increase  the  operating  temperature  of  the  storage 
medium  and  may  ultimately  result  in  a  room  temperature  storage 
capability  [2].  Despite  this  progress,  there  still  remains  the 
requirement  that  large  amounts  of  hydrogen  need  to  be  stored,  yet 
the  low  hydrogen  surface  packing  density  (SPD)  of  many  carbon- 
based  materials  has  limited  their  adsorption  capabilities. 

Tlie  metal  organic  framework  MOF- 74  resembles  a  bundle  of 
hexagonally  packed  drinking  straws  with  edge  sharing  ZnO^ 
octahedra  along  the  intersecting  walls  of  three  neighboring  tubes 
and  the  remaining  'straw'  composed  of  the  organic  ligand  linking 
the  metal  oxide  clusters.  Unlike  many  MOF  materials,  MOF- 74 
contains  exposed  Zn-*  metal  ions  that  can  coordinate  to  solvents 
or  even  hydrogen.  Gas  adsorption  studies  indicate  a  rather  large 
hydrogen  enthalpy  of  adsorption  of -8.8  kj/mol,  and  a  relatively 
large  excess  hydrogen  uptake  capacity  of  2.8  %  mass  fraction  for  a 
modest  surface  area  of  870  mVg. 
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FIGURE  1:  A  (001)  view  of  the  real-space  Fourier-difference  scattering- 
length  density  superimposed  with  hexagonal  pore  structure  of  MOF-74, 
indicating  the  location  of  the  first  adsorption  sites  (red-yellow  regions)  and 
the  second  adsorption  sites  (yellow  regions). 
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Neutron  powder  diffraction  experiments,  using  the  BT-1 
diffractometer,  were  performed  to  understand  the  local 
environment  at  and  around  the  exposed  Zn^*  sites  in  MOF-74 
[3].  A  diffraction  pattern  was  collected  for  the  desolvated  material 
prior  to  adding  D,.  Subsequent  Fourier  difference  mapping 
(Rg.  1)  and  Rietveld  analysis  indicated  the  existence  of  extra 
neutron  scattering  length  density  close  to  the  Zn'*  ions  in  the 
framework  skeleton,  as  expected  for  hydrogen  initially  adsorbed 
at  the  CUMC,  with  an  interaction  distance  of  2.6  A.  Upon 
additional  loading  of  deuterium  we  can  further  identify  two 
more  sites  that  generate  a  densely  packed  tube  of  hydrogen 
■within  the  pores  of  MOF-74  (Fig.  2).  A  more  detailed  study 
of  the  packing  geometry  of  the  deuterium  molecules  reveals 
that  the  second  site  lies  above  an  oxygen  triangle  while  the 
third  site  is  associated  with  the  benzene  ring  of  the  linker  with 
deuterium-framework  distances  larger  than  3  A,  as  is  typical 
for  the  weak  hydrogen  interactions.  Surprisingly,  the  second 
and  third  adsorbed  deuteriums  both  exhibit  tremendously  short 
distances  of  2.85  A  and  2.90  A  to  the  first  adsorbed  deuterium, 
respectively.  These  individual  distances  are  much  shorter  than  the 
3.6  A  intermolecular  distance  found  for  soUd  D,  under  similar 
measurement  conditions  and  without  the  application  of  pressure 
in  both  cases.  Although  the  diffraction  measurements  were 
performed  at  4  K,  the  D,-D^  intermolecular  spacings  between  the 
first  three  adsorption  sites  shed  light  on  the  relatively  large  SPD 
in  MOF-74  at  77  K. 


FIGURE  2:  One  dimensional  nanoscale  tube-like  structure  formed  from 
adsorbed  D^.     molecules  (large  green  spheres]  adsorbed  in  the  hexagonal 
pores  of  MOF-74  (partial  cross-section  shown  for  clarity.  Zn  :  blue,  C  :  grey, 
0  :  red,  H  :  pale  green)  form  a  one-dimensional  nanoscale  tube-like  structure 
from  the  first  three  adsorption  sites.  Deuterium  molecules  are  connected  if 
the  distance  between  D.,  molecules  is  less  than  3.65  A. 


the  SPD  as  the  saturated  excess  mass  adsorption  of  per 
unit  surface  area  measured  at  77  K  (e.g.,  the  nitrogen  derived 
Brunauer-Emmett-Teller  'BET' or  Langmuir  surface  areas). 
A  plot  of  the  SPD  value  for  MOF-74  and  a  selection  of  other 
framework  materials  is  given  in  Fig.  3.  Evidently,  the  presence 
of  unsaturated  metal  centers  in  MOFs  not  only  increases  the 
enthalpy  for  hydrogen  adsorption,  but  also  allows  for  hydrogen 
molecules  to  be  held  closer  together  when  adsorbed.  Surprisingly, 
in  the  case  of  MOF-74  the  hydrogen  molecules  are  even  closer 
together  at  77  K  than  in  solid  hydrogen  at  4  K,  in  agreement  with 
the  diffraction  data. 
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FIGURE  3:  The  surface  packing  density  [SPD]  of  hydrogen  in  metal-organic 
frameworks  is  typically  less  than  that  obtained  in  actiuated  carbon 
adsorbents  (dashed  black  line).  MOFs  with  SPD  values  above  the  line  contain 
CUMCs.  Excess  hydrogen  adsorption  capacities  are  indicated  by  the  dashed 
blue  and  solid  green  curves. 

In  this  study,  we  have  demonstrated  a  strong  correlation  between 
exposed  CUMCs  in  framework  structures  and  larger  associated 
SPD's.  These  observations  provide  an  avenue  to  improve  the 
gravimetric  H,  uptake  by  increasing  the  surface-packing  density 
of  hydrogen.  The  study  of  MOF-74  with  neutron  powder 
diffraction  alongside  isotherm  measurements  shows  that 
adsorbed  on  the  MOF-74  surface  has  a  higher  surface  density  at 
77  K  than  that  of  solid  H  at  ~  4  K  and  zero  pressure. 
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In  order  to  evaluate  the  hydrogen  surface  density  properties  of 
different  adsorbents,  we  use  the  experimental  surface  packing 
density  obtained  from  gas  adsorption  studies.  We  can  define 
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structure  of  Li^CaCNH)^  and  Hydrogen  Storage  Mechanisms  in 
the  Amide-Hydride  System 


Hui  Wu' 


Suitable  hydrogen  storage  materials  with  the  on-board 
operating  capabilities  for  fuel-cell  vehicular  applications 
is  one  of  the  major  challenges  to  widespread  use  of  hydrogen 
as  a  primary  fuel  and  to  reduce  dependence  on  fossil  fuel 
and  emissions  of  greenhouse  gases.  LiNH^  with  a  hydrogen 
storage  potential  up  to  10.4  %  mass  fraction  has  been  viewed 
as  a  promising  hydrogen  storage  system  [1].  However,  its 
practical  application  for  hydrogen  storage  is  limited  because 
of  its  low  hydrogen  equilibrium  pressure  (<  0.01  bar)  and 
high  desorption  temperature  (Tj^^  >  320  °C).  Another 
feature  of  this  system  that  complicates  the  decomposition 
mechanism  is  the  competing  release  of  ammonia  from 
LiNH,  at  high  temperature.  To  reduce  the  desorption/ 
absorption  temperatures  (T^J  T^^^)  of  pure  LiNH^/Li^NH, 
in  addition  to  the  use  of  catalysts  and  reduced  particle  size, 
several  investigations  have  been  conducted  on  the  effects  of 
hydride  additives  (LiH,  MgH^,  CaH^)  on  the  ratio  T^JT^^_^ 
and  the  amount  of  ammonia  released  [2,3].  In  some  cases  it 
has  been  claimed  that  the  ratio  T,  IT ,  of  the  mixed  amide/ 

des  abs 

hydride  systems  can  be  significantly  reduced  compared  to  the 
results  for  pure  LiNH,.  However,  the  hydrogen  storage  and 
release  mechanism  of  these  amide/hydride  mixed  systems 
has  not  been  established  mainly  due  to  the  lack  of  accurate 
crystaUographic  information  on  the  ternary  imides  formed 
after  dehydrogenation. 

In  our  recent  study  in  LiNH^-CaH,  system  [4],  we  have 
determined  the  crystal  structure  of  Li,Ca(ND),  by  neutron 
powder  diffraction  (NPD),  and  proposed  a  mechanism  for 
the  hydrogenation/ dehydrogenation  of  the  mixed  amide  and 
hydride  system.  This  finding  holds  the  key  to  understanding 
the  hydrogen  storage  mechanism  in  such  systems,  and  is 
critical  for  rational  development  of  better  candidates  for 
hydrogen  storage. 

Figure  1  shows  the  crystal  structure  of  Li.CaCND)^  with  D 
randomly  distributed  on  one  of  the  three  sites  around  each  N 
atom.  In  this  structure  each  Ca  is  coordinated  with  six  ND'" 
anions  forming  a  nearly  regular  Ca[ND]^  octahedron  with 
six  equal  Ca-N  bonds  (2.527  A);  each  Li  is  bonded  with  four 
ND^"  anions  with  Li-N  bond  lengths  in  a  range  of 


2.167  A  to  2.291  A.  The  large  difference  in  the  size  and 
coordination  preference  of  Li  and  Ca  drives  the  formation  of 
an  ordered  arrangement  of  alternating  Ca[NH]^-octahedra  and 
Li[NH]^-tetrahedra  along  the  c  axis.  Therefore,  Li2Ca(NH)2 
can  actually  be  viewed  as  a  combined-imide  structure  consisting 
of  two  different  imide  layers.  Alternatively,  Li2Ca(NH)^  can  be 
described  as  a  layered  structure  consisting  of  infinite  2D  slabs 
of  edge-shared  Ca[NH]^-octahedra,  which  are  separated  by  the 
motif  of  Li  cations  (Fig.  2),  as  in  many  "Li-intercalated"  layered 
oxides. 
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FIGURE  1:  Experimental  (circles],  calculated  (red  line),  and  difference  (grey 
line)  I\1PD  profiles  for  Li2Ca(IMD)g  at  1 5  K.  The  patterns  also  contain  peaks 
from  small  amounts  of  LilMDg  (1.09  %  mass  fraction]  and  CaiMD 
(2.38  %  mass  fraction].  Vertical  bars  indicate  the  calculated  positions  of 
Bragg  peaks  for  Li2Ca(l\ID]g,  LiND^  and  CaND  (from  the  top),  respectively. 
A  =  1 .5403  A.  Inset:  (Left)  off-[1 1 0]  view  of  the  trigonal  structure  of 
Li2Ca(NH]2  (P3m1).  Ca(NH)g-octahedra  are  in  green;  Li(NH]^-tetrahedra  are 
in  yellow.  Ca,  Li  and  N  atoms  are  represented  by  large  pink,  yellow,  and  blue 
spheres,  respectively.  H  atoms  are  randomly  distributed  at  one  of  the  three  . 
white  sites  around  each  IM  atom.  (Right)  Layered  structure  of  Lt2Ca{IMH)g  viewed 
as  a  "combined  imide"  consisting  of  ordered  CaNH-layer  and  Li^NH-layer. 

To  understand  the  formation  of  the  Li^Ca(NH)^,  intermediate 
products  at  different  dehydrogenation  stages  of  2LiNH^  +  CaH^ 
mixture  were  monitored  using  x-ray  diffraction  (XRD),  and 
the  hydrogen  contents  were  determined  by  the  prompt  gamma 
neutron  activation  (PGAA)  technique.  In  the  early  stage  of 
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desorption  the  sample  contains  a  multi-phase  mixture,  including 
CaNH,  Li,NH,  CaH„  and  LiNH,.  Upon  further  desorption, 
CaNH,  Li^NH  and  LiNH,  remain,  with  the  continued  presence 
of  Li^Ca(NH),  and  the  absence  of  CaH^.  After  complete 
dehydriding,  the  product  contains  single-phase  Li,Ca(NH)2. 
Rehydrogenating  this  product  results  in  a  mixture  of  CaNH, 
LiNH^  and  LiH,  differing  from  the  initial  mixture 
2LiNH2  +  CaH,.  From  these  results,  we  proposed  a  mechanism 
for  dehydrogenation  of  the  2LiNH,  +  CaH,  mixture, 
hydrogenation  of  Li,Ca(NH)2,  and  desorption  of  the  re- 
hydrogenated  product,  which  we  believe  best  explains  our 
structure  results  and  the  chemistry  of  the  system. 


FIGURE  2:  Structure  of  Li^CafNHJj  viewed  as  an  "intercalated"  layered 
structure  consisting  of  2D  slabs  of  Ca(NH)g-octahedra  separated  by  Li  ion 
motifs.  Ca,  Li,  N  and  H  atoms  are  colored  as  in  Fig.  1.  The  facile  motion  of 
Li  assists  the  hydrogenation  process  with  totally  1      absorbed  and  final 
products  of  CalMH.  LiNH^  and  LiH. 

Previous  study  in  LiNH^  indicated  the  presence  of  Frenkel 
defect  pairs  (a  charged  interstitial  [Li.LiNH^]*  and  a  lithium 
vacancy  [Dj^^NH^]")  and  the  resulting  motion  of  Li*  and  H*  [5]. 
CaHj  with  an  anti-PbCl^  structure  is  known  as  an  electrolyte 
material  with  highly  mobile  H'  ions  (also  created  by  the  Frenkel 
defect  pair  interstitial  [CaH^H^]"  and  vacancy  [CaHDj^]*). 
Therefore,  both  mobile  Li*  and  H*  in  LiNH^  and  H"  in  CaH^ 
are  involved  in  the  dehydrogenation  at  high  temperature.  Upon 
dehydrogenation  of  LiNH^/CaH^  (see  scheme  1  below),  the 
highly  mobile,  free  H"  anion  in  CaH^  and  protonic  H*  in  LiNH^ 
can  very  easily  combine  to  create  H,  gas  thanks  to  a  very  high 
reaction  enthalpy.  The  fast  combination  of  H*  and  H'  leaves 
[□^.NHJ"  and  [CaHD^J*  in  LiNH,  and  CaH,,  respectively,  and 
drives  the  reaction  between  these  two  defects  on  the  surface  of 
adjacent  particles.  When  these  processes  occur,  a  mixture  of 
2LiNH2  +  CaH^  will  release  2H.^  in  total,  consistent  with  the 
hydrogen  amount  observed.  More  importantly,  the  participation 
of  H'  enhances  interaction  between  H  and  H*  that  promotes 
hydrogen  to  be  released  at  lower  temperature,  and  suppresses  the 
interaction  between  H*  and  NH^'  that  produces  NH^. 

During  hydrogenation  of  the  layered  structure  Li2Ca(NH)2 
(Fig.  2),  in  the  layer  of  Li^NH,  mobile  Li*  reacts  rapidly  with  H^, 
forming  LiH  and  a  proton.  The  proton  produced  will  then  react 


with  the  negatively  charged  Li-vacancy  and  bond  with  NH'  to 
form  a  NH^"  group.  In  the  layers  of  CaNH  there  is  no  mobile 
species,  i.e.,  protons  in  NH'  groups  are  covalently  bonded  and 
Ca'*  cations  are  too  heavy  and  too  large  to  migrate.  Consequently, 
the  "inert"  CaNH  remains.  Continuous  ion  movement  and 
interactions  in  Li^NH-layer  will  finally  dissociate  Li2Ca(NH)2 
into  LiNH^,  LiH  and  unreacted  CaNH,  as  observed  in  our  XRD/ 
NPD  results  [4].  Li,Ca(NH)_,  could  absorb  H^  at  temperatures 
«  130  °C  lower  than  that  of  Li,NH  [3,4].  Considering  the 
structural  characteristics  of  the  ternary  imide,  we  believe  that  it 
is  the  special  layered  structure  with  "intercalated"  mobile  Li  ions 
that  is  mainly  responsible  for  the  dramatically  lowered  T  j^^.  The 
facile  mobility  of  Li  within  the  2D  channels  is  also  confirmed  by 
the  observed  increasing  nonstoichiometry  at  Li  lattice  sites  with 
elevated  temperature. 


Step  1 ;  C  reation  of  Frenkel  defect  pairs 
Li*  +  LiNH, 


LiNH;  +  LiNH^     [LiiLiNHj]*  +  [D^^NH^]-  h' +  K h,(-aE)  K^^HG  „]*+  [CaH^Hj"  <-*  CaHj  +  CaH^ 


Step  2:  tons  immigration  & 

Defect  reactions        (H*)+jLi  NH  I 


jCaNHj+H; 


Step  3:  Reaction  between  two  imides    LijNH  +  CaNH  ->  LijCa(NH); 

Scheme  1:  Dehydrogenation  of  2  LilMH^  +  CaH^  with  a  total  of  2H^  released 
and  products  of  CaNH,  Li^NH  and/or  Li^CafNHJg  observed  at  various 
desorption  steps. 

In  summary,  the  structure  of  ternary  imide  Li,Ca(NH)2  was 
determined  using  NPD  on  a  deuterated  sample.  Li,Ca(NH), 
crystallizes  in  a  layered  structure  consisting  of  infinite  2D  slabs 
of  edge-shared  Ca[NH]^-octahedra  separated  by  the  motif  of  Li 
cations.  The  movement  of  small  mobile  species  such  as  Li*,  H* 
and  H'  has  been  shown  to  play  a  key  role  in  both  decomposition 
and  hydrogenation  of  LiNH/CaH^.The  major  benefit  of  the 
participation  of  H"  is  the  enhanced  interaction  between  H"  and 
H*  that  promotes  hydrogen  release  at  a  lower  temperature,  and  a 
reduction  in  the  ammonia-producing  interaction  between  H*  and 
NH^  .  The  resultant  ternary  imide  with  a  layered  structure  allows 
a  rapid  Li*  movement,  leading  to  a  much  lower  T^^^  compared  to 
Li^NH.  Therefore,  we  conclude  that  the  differently  charged  small 
mobile  species  are  primarily  responsible  for  the  reduced 
T,  /T ,  and  for  the  reduced  ammonia  release  in  the  mixed 

dcs  :ibs 

amide  /  hydride  system. 
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Melilite  as  an  Interstitial  Oxygen  Conducting  Electrolyte 
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Jantsky\  M.J.  Rosseinsky^ 


Fuel  cells  date  to  the  mid-19th  century  when  the  concept 
was  first  proposed  by  the  Swiss  chemist,  Schonbein. 
They  have  been  widely  applied  through  the  years,  including 
providing  both  electricity  and  drinking  water  (as  a  by- 
product) in  NASA  space  programs.  As  opposed  to  batteries 
which  generate  electricity  from  stored  chemical  energy, 
fiiel  cells  generate  electricity  directly  by  electrochemically 
combining  fuel  and  oxygen.  Electricity  is  produced  as 
hydrogen  moves  from  the  anode  side  to  combine  with 
oxygen  moving  from  the  cathode  side.  The  reaction  product 
is  environmentally  benign  water.  Although  they  have  greatly 
improved  since  the  alkali-based  ceils  used  in  the  Apollo  space 
program,  fuels  cells  are  still  not  competitive  in  cost  or  power 
output  with  the  ubiquitous  internal  combustion  engine. 

The  development  of  more  efficient  systems  requires  the 
identification  of  new  materials.  The  ceramic  Solid  Oxide 
Fuel  Cell  (SOFC)  is  a  particularly  promising  variety 
and  is  characterized  by  very  favorable  power  output  but 
at  the  cost  of  high  operating  temperature  (~  850  °C),  a 
prerequisite  to  generating  sufficient  oxygen  mobility.  The 
time  presently  taken  to  heat  a  typical  commercial  fuel  cell 
to  its  operating  temperature  is  around  45  min  and  the 
US  Department  of  Energy  would  like  this  to  fall  to  less 
than  2  min.  Oxide  mobility  in  solids  is  associated  with 
defects.  Although  anion  vacancies  are  the  charge  carriers 
in  most  cases,  excess  (interstitial)  oxide  anions  give  high 
conductivities  in  isolated  polyhedral  anion  structures  such  as 
the  apatites.  The  development  of  new  families  of  interstitial 
oxide  conductors  with  less  restrictive  structural  constraints 
requires  an  understanding  of  the  mechanisms  enabling  both 
incorporation  and  mobility  of  the  excess  oxide. 

Neutron  diffraction  (ND)  allowed  us  to  show  how  the 
two-dimensionally  connected  tetrahedral  gallium  oxide 
network  in  the  melilite  structure,  La^    Sr^^^Ga^O^^y  (Fig- 1)> 
stabilizes  oxygen  interstitials  by  local  relaxation  around 
them,  affording  an  oxide  ion  conductivity  of  0.02  S  cm^^  to 
0.1  S  cm"'  over  the  600  °C  to  900  °C  temperature  range  [1]. 


To  identify  the  mobile  oxygen  species  in  La^^^  Sr^  ^^Ga30^  ^y)  we 
collected  ND  data  on  the  BT-1  diffractometer  at  the  NCNR. 
Refinement  showed  that  all  sites  were  fuUy  occupied.  Difference 
Fourier  map  calculations  showed  positive  scattering  density  within 
the  tetrahedral  layers  at  the  Gal-Ga2  level  between  the  La/Sr 
cations  along  the  c-axis  direction  (labeled  04  in  Figs.  1  and  2) 
the  occupancy  of  which  was  refined  to  0.136(5).  The  refined  total 
oxygen  content  of  7.27(1)  per  formula  unit  agrees  well  with  the 
oxygen  over-stoichiometry  obtained  through  other  measurements. 
Refinement  of  anisotropic  atomic  displacement  showed  a 
considerable  amount  of  positional  disorder,  presumably  associated 
with  the  interstitial  oxygen. 


La/Sr 


Gal 


Ga2 


structure.  The  04  positions  indicate  additional  interstitial  oxygen  ions,  not 
found  in  the  parent  structure. 

Further,  we  performed  maximum-entropy  method  (MEM) 
analysis,  known  to  minimize  a  bias  imposed  by  the  strucmral 
model  within  Rietveld  refinement,  which  confirmed  the  existence 
of  the  oxygen  interstitial  (see  Fig.  2).  The  additional  scattering 
density  around  all  the  atom  sites  is  assigned  to  positional  disorder 
associated  with  the  local  structural  relaxations  required  to 
accommodate  the  interstitial  oxygen. 
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In  order  to  address  the  possible  oxide  diffusion  pathway, 
a  comparison  can  be  made  between  the  ambient  and  high 
temperature  (800  °C)  ND  data.  Rien-eld  refinement  of  the  800  °C 
ND  data  re\'eals  the  interstitial  ox\'gen  is  essentially  constrained 
around  the  4e  sites  in  the  tunnels.  But  the  thermal  parameters 
for  04  are  now  extremeh^  flattened  within  the  laver,  reflecting 
predominant  motion  \\ithin  the  gallate  layer.  The  scattering 
density  of  the  framework  ox\'gen  atoms  shows  elongation  in 
the  MEM-derived  densitv^  compared  with  that  at  ambient 
temperature.  At  800  °C,  the  additional  scattering  densitv"  in  the 
ME^VI  analysis  is  due  to  the  gallate  framework  vibrational  motion 
with  no  apparent  additional  scattering  densit}-  between  the 
neighboring  04  sites  and  between  04  and  bulk  oxygen  sites.  This 
is  consistent  with  a  direct  interstitial  (single  particle)  mechanism 
for  the  ox)'gen  diffusion  rather  than  the  direct  involvement  of 
framework  oxygen  atoms. 


Ambient  temperature  800  =C 


families  of  materials  as  candidate  interstitial  oxide  conductors. 


FIGURE  3:  Structural  relaxation  around  the  interstitial  oxide.  The  bulk 
structure  (a)  and  local  defect  structure  (b)  around  the  oxygen  interstitial 
04.  [0^\  is,  unlike  (03)^^,  not  directly  bound  to  (Ga2)^  but  relaxes  away  from 
interstitial  04,  accommodated  by  a  change  in  bond  and  bridging  angles  at 
Ga2.  The  02  above  the  two  Ga2  centers  was  omitted  for  clarity. 
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FIGURE  2:  Comparison  of  the  nuclear  scattering  density  of 
La.i  ^Sr^^gGa^O^g?      obtained  from  maximum  entropy  analysis  maps. 
The  (001)  projections  [top]  and  corresponding  scattering  density  at  z  =  0 
[below]  are  shown  at  ambient  (left)  and  800  "C  (right). 

A  simple  model  for  relaxation  around  the  defect  that  involves 
those  Ga2,  La/Sr  01  and  03  species  neighboring  04  being 
displaced  to  {Ga2)^,  (La/Sr)^^,  (Ol)^^  and  (03)j_  sites  was  refined. 
The  resulting  bulk  structure  and  the  local  defect  structure  at  the 
04  interstitial  oxygen  are  shown  in  Figs.  3a  and  3b,  respectively. 
04  is  accommodated  by  the  synergy  of  the  change  in  bonding 
geometry  at  Ga2,  modification  of  polyhedral  internal  and 
bridging  angles,  and  counterion  displacement.  The  displacement 
of  (Ga2)j^  towards  04  shortens  this  gallium-oxygen  bond  from 
2.13(2)A  in  the  average  structure  to  1.81(2)A,  clearly  showing 
that  the  interstitial  O  enters  the  coordination  environment  of  one 
of  the  two  distinct  framework  gallium  centers.  This  relaxation 
process  shows  the  ability  of  the  melilite  structure  to  accommodate 
the  interstitial  oxygen  ion  via  local  lattice  relaxation  as  well  as 
to  sustain  oxygen  mobility,  which  opens  up  many  new  structural 
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Quasi-free  Methyl  Rotation  in  Zeolitic  Imidazolate 
Framework-8  (ZIF8) 
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In  general,  the  accurate  determination  of  interatomic 
potentials  in  solids  is  very  challenging.  However,  there 
are  cases  where  molecules  within  solids  undergo  sharp 
quantum  tunneling  transitions  which,  because  they  depend 
exponentially  on  these  potentials,  provide  a  stringent  test 
of  modern  computational  methods.  Direct  measurement 
of  these  transitions  by  neutron  spectroscopy  thus  provides 
a  very  powerful  probe  of  interatomic  potentials.  Especially 
rare  and  important  are  those  cases  where  a  low  barrier  to 
rotational  tunneling  provides  exquisite  sensitivity  to  the 
interatomic  potential.  Here  we  highlight  our  recent  work 
on  such  a  system,  ZIF8  [1]. 

Zeolitic  imidazolate  frameworks  (ZIFs)  [2,3]  are  a 
recently  developed  class  of  compounds  made  up  of  MN^ 
(M  =  Co,  Cu,  Zn,  etc.)  clusters  linked  together  with 
imidazolate  ligands  to  provide  tunable  nanosized  pores. 
ZIFS  (Zn[MeIM]„  MelM  =  2-methylimidazolate)  is  a 
prototypical  ZIF  compound  [4].  An  intriguing  feature  of 
the  ZIFS  structure  (Fig.  la)  is  its  methyl  groups,  bound  to 
the  imidazolate  ligands  (Fig.  lb).  The  relatively  isolated 
methyl  group  in  ZIFS  is  oriented  toward  the  large  cavity 
of  the  porous  ZIFS  structure  and  thus  exhibits  quasi-free 
rotation,  restricted  only  slightly  by  interatomic  barriers. 
Tunneling  of  the  methyl  groups  through  these  barriers  at 
low  temperatures  probes  the  interatomic  potential  in  detail. 


FIGURE  1:  (a)  ZIFS  crystal  structure  of  (space  group  /43m).  The 
3-D  open  framework  consists  of  ZnN^  clusters  [the  tetrahedra) 
connected  by  2-methylimidazolate  (MelM)  linkers.  The  central  space- 
filling sphere  denotes  the  pore  structure.  For  clarity,  only  the  H  atoms 
associated  with  the  methyl  groups  are  shown,  (b)  The  structure  of 
the  MelM  linker,  (c)  The  ellipsoids  of  the  MelM  linker  obtained  from 
ZIFS  neutron-diffraction  data  at  3.5  K. 


Neutron  diffraction  measurements  of  ZIFS  and  deuterated  ZIFS 
(D-ZIFS)  using  the  BT-1  diff tactometer  indicated  that  the  host 
lattice  maintains  the  same  structure  between  3.5  K  and  300  K. 
It  should  be  noted  that  the  torsional  amplitudes  of  the  methyl 
groups  are  quite  large  (Fig.  Ic),  but  they  suggest  a  mainly  3-fold 
character  for  the  rotational  potential.  Torsional  energies 
E     2.9  meV  for  CH^  and    1.7  meV  for  CD,  at  3.5  K) 
were  estimated  from  diffraction-derived  torsional  amplitudes 
within  the  harmonic  approximation,  (although  the  "harmonic 
assumption"  is  clearly  questionable  for  large  rotational 
amplitudes). 

For  1-D  rigid  rotation,  it  is  well  established  that  the  quantum 
rotational  levels  are  determined  by  the  stationary  Schrodinger 
equation.  For  a  free  CH^  rotor,  the  /=  0^1  transition  energy 
(splitting  of  the  ground  state  rotational  level)  is  equal  to 
0.655  meV.  For  a  non-zero  rotational  barrier,  this  transition 
energy  decreases  with  increasing  barrier  height.  To  quantitatively 
probe  the  methyl-group  rotations,  we  measured  the  O'^l 
rotational  transition  of  ZIFS  by  neutron  spectroscopy  at 
334  ±  1  peV  (Fig.  2).  The  lack  of  mode-splitting  indicates  little 
methyl-methyl  coupling.  The  corresponding  CD,  tunneling 
transition  in  D-ZIFS  was  measured  with  both  the  HFBS  and 
DCS  spectrometers  at  33  ±  1  [aeV  (at  much  lower  energy  due 
to  the  2x-higher  CD,  moment  of  inertia  compared  to  CH,). 
Additional  DCS  measurements  for  ZIFS  showed  a  broad  peak 
at  ^  2.7  ±0.1  meV,  which  we  tentatively  assigned  to  the /=  0^2 
level,  although  framework  modes  could  also  contribute  here.  If 
we  assume  a  pure  3-fold  barrier  to  rotation,  the  334  peV  tunnel 
splitting  for  ZIFS,  combined  with  a  second  rotational  level  of 
2.7  meV,  predicts  a  very  low  barrier  height  of  6  meV.  The  33  \ieV 
tunnel  splitting  for  D-ZIFS  suggests  a  slightly  higher  barrier  of 
about  S  meV  for  the  CD,  groups. 

At  higher  temperatures  (~  30  K  and  above),  the  tunneling  gives 
way  to  classical  jump  reorientation  as  the  peaks  shift  toward 
the  elastic  line  and  broaden  (Fig.  2  inset).  Figure  3  shows  the 
temperature  dependence  of  the  quasielastic  linewidths  {f)  for 
ZIFS  plotted  in  an  Arrhenius  fashion.  The  activation  energy 
for  3 -fold  jump  reorientation  can  be  estimated  from  the  slope 
under  the  assumption  that  the  methyl  groups  reorient  through 
a  jump  diffusion  process  by  thermal  activation  over  the  barrier. 
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The  transformation  from  3-fold  jump  to  rotational  diffusion  is 
manifested  by  the  nonlinear  nature  of  the  Arrhenius  plot  with 
increasing  temperature.  We  estimated  an  activation  energy  of 
«  3  meV  from  the  slope  associated  with  the  data  between  30  K 
and  40  K.  (At  higher  temperatures  (>  50  K),  this  model  loses 
its  validity  as  the  thermal  energy  is  comparable  to  the  barrier 
height.)  Further  assuming  that  this     was  a  measure  of  the 
distance  from  the  first-excited  rotational  (J  =  2)  level  (tentatively 
assigned  at  2.7  meV)  to  the  top  of  the  potential  barrier,  and 
estimating  a  zero-point  energy  of  %(2.7  meV),  led  to  a  barrier 
height  of  about  7  meV,  close  to  that  predicted  above  for  a  3-fold 
potential. 
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FIGURE  2:  ZIF8  tunneling  spectrum  at  1.4  K  using  the  DCS  spectrometer 
with  7  A  neutrons  (=  18  peV  resolution},  (inset)  Euolution  of  the  tunneling 
spectrum  with  temperature  using  the  PCS  spectrometer  with  6  A  neutrons 
(k  70  peV  resolution). 
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3:  A  plot  of  Inir)  vs.  T  '  for  ZIFB  between  30  K  and  90  K,  where  /' 
is  the  full-width-at-half-maximum  quasielastic  linewidth. 


To  obtain  further  insight,  we  calculated  the  rotational  barrier 
using  density  fiinctional  theory  (DFT).  From  the  optimized 
structure,  we  performed  total-energy  calculations  as  a  function 
of  the  rotational  angle  of  the  methyl  group.  For  simplicity,  we 
only  rotated  one  methyl  group  in  the  unit  cell.  Interestingly, 
the  calculated  results  (Fig.  4),  which  can  only  be  considered 
qualitative,  predict  a  shape  of  the  potential  that  includes 
both  3-fold  and  smaller  6-fold  terms  and  a  barrier  height  of 
~  6.6  meV,  consistent  with  the  value  determined  from  the  various 
neutron  scattering  measurements. 

Thus,  all  the  results  of  our  study  show  the  existence  of  an  unusual 
(for  a  solid)  quasi-free  rotational  potential  for  CH^  groups 
in  ZIF8.  The  high  methyl-group  tunneling  energies  and  low 
hindered  rotational  transitions  combined  with  the  results  of 
quasielastic  neutron  scattering  provide  a  rotational  barrier  of 
about  7  meV,  largely  3 -fold  and  consistent  with  first-principles 
calculations.  Two  previous  examples  of  quasi-free  methyl 
rotation  are  4-methyl  pyridine  [5]  and  Li  acetate  dehydrate  [6], 
which  have  high  tunneling  transitions  and  low  rotational  barriers. 
Yet,  both  show  methyl-methyl  coupling  effects,  which  result  in 
complex  tunneling  spectra.  In  contrast,  the  ZIF8  spectrum  shows 
no  evidence  of  sidebands,  which  indicates  "independent"  rotation 
of  CHj  groups  and  therefore  provides  a  "textbook"  example  of 
rotational  tunneling  in  a  low-barrier  system. 
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FIGURE  4:  Rotational  potential  of  the  methyl  group  in  ZIF8,  obtained  from 
DR  calculations.  Dotted  lines  are  the  fit  to  a  3-fold/6-fold  mixed  potential. 
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Recent  advances  in  the  synthesis  of  hydrogels  have 
resulted  in  materials  with  fracture  toughness 
comparable  to  that  of  articular  cartilage  [l].The 
toughness  arises  from  the  combination  of  a  stiff  but  brittle 
crosslinked  polyelectrolyte  (poly-2-acrylamido-2-methyl- 
1-propanesulfonic  acid,  PAMPS)  and  a  soft  and  pUable 
neutral  linear  polymer  (polyacrylamide,  PAAm).This  class 
of  hydrogels  has  been  named  double-network  hydrogels 
(DN-gels).  Although  many  polymer  pairs  have  been  tested, 
only  a  few  give  rise  to  hydrogels  with  superior  toughness. 
This  observation  suggests  that  some  specific  composition- 
dependent  associations  between  polymer  pairs  in  the 
DN  hydrogel  may  play  an  important  role  in  shaping  the 
toughness  of  DN-gels. 

In  our  previous  work,  small  angle  neutron  scattering  was 
used  to  measure  the  molecular  structure  of  DN-gels  under 
static  conditions  as  well  as  under  deformation  [2,  3].  The 
surprising  outcome  of  these  experiments  is  that  both  the 
PAMPS  and  the  PAAm  become  more  homogeneous  when  in 
the  presence  of  each  other  than  in  pure  water  by  themselves. 
This  result  suggested  that  there  could  be  associations  such 
as  hydrogen  bonding  between  the  charged  and  uncharged 
molecules  [4].  We  further  noticed  that  even  within  DN- 
gels  made  of  Pt^MPS  and  PAAm,  only  those  with  a  volume 
ratio  near  1:7  (PAMPS  :  PAAm)  exhibit  extraordinarily 
tough  mechanical  strength  [2].  This  result  suggests  that 
the  interchain  association  strength  depends  not  only  on  the 
chemical  structure  of  the  constituents  but  also  on  the  molar 
ratio  between  these  two  components. 

In  order  to  clarify  the  nature  and  the  origin  of  interchain 
associations  between  PAMPS  and  PAAm,  we  performed 
a  neutron  spin-echo  (NSE)  study  to  determine  the  chain 
dynamics  over  a  broad  scale  of  frequency  and  length  in  the 
PAMPS  and  PAAm  solution  blends.  The  solution  blends 
were  prepared  by  mixing  linear  polymer  solutions  of  PAMPS 
and  d^FAAm  at  various  volume  ratios.  Pure  PAMPS  and 
d^FAAm  solutions  were  prepared  by  UV-initiated  free-radical 
polymerization  of  their  monomers  in  D^O  with  2-oxoglutaric 
acid.  By  adjusting  the  ratio  between  D^O  and  H^O,  we 


brought  the  neutron  scattering  contrast  between  ^/SPAAm 
and  the  solvent  to  zero,  thus  the  dynamics  of  PAMPS  chains 
could  be  measured  using  NSE.  In  this  study,  the  total  polymer 
concentration  of  the  mixtures  is  kept  constant  (0.5  mol/L)  while 
the  concentration  of  PAMPS  and  PAAm  are  varied  (PAMPS: 
(0.375, 0.125,  0.0625,  0.03125)  mol/L).  NSE  measurements  were 
conducted  at  25  °C  on  the  NG5-NSE  beamline  at  the  NCNR 
using  a  6  A  wavelength  incident  beam. 

Figure  1  shows  the  intermediate  scattering  functions  for  the  pure 
0.5  mol/L  aqueous  PAMPS  solution.  Over  the  entire  ^'-range  the 
results  can  be  modeled  as  a  single  exponential  function,  i.e.,  the 
relaxation  time  (r^lq))  can  be  obtained  by  fitting  the  equation 


=  ^exp(-//T^,(g)) . 


FIGURE  1:  Partial  Intermediate  scattering  function  of  PAMPS  in  0.5  mol/L 

PAMPS  solution. 

In  all  solution  blend  samples,  the  interrrtediate  scattering 
functions  are  also  single  exponential  functions  similar  to  that  of 


^Polymers  Division,  National  Institute  of  Standards  and  Technology,  Gaithersburg,MD  20899 

^NIST  Center  for  Neutron  Research,  National  Institute  of  Standards  and  Technology,  Gaithersburg,  MD  20899 

^Indiana  University  Cyclotron  Facility,  Bloomington,  IN  47408 


the  pure  PAMPS  solution.  This  observation  indicates  that  the 
presence  of  neutral  polymers  does  not  induce  any  unexpected 
dynamics  of  the  PAMPS  chain  over  all  the  y-range  probed. 
The  5?-dependent  relaxation  spectrum  of  PAMPS  chain  in 
PAMPS/  PAAm  solution  blends  seems  to  follow  quaUtatively 
the  static  structure  S^Jiq)  measured  by  small  angle  neutron 
scattering  (SANS).  However,  it  is  interesting  to  note  that  the 
relaxation  time  measured  at  the  lowest  accessible  q  ~  0.05 
A"^  (Fig.  2)  exhibits  a  maximum  in  the  solution  blend  with 
0.0625  mol/L  PAMPS.  This  result  is  surprising  since  the  static 
structure  revealed  by  SANS  does  not  exhibit  an  anomaly  at  this 
composition;  the  normalized  low  q  intensity  of  SANS  simply 
increases  gradually  with  a  decrease  in  PAMPS  concentration  [2]. 
The  slower  relaxation  at  0.0625  mol/L  PAMPS  concentration 
is  significant  because  it  is  close  to  the  PAMPS:PAAm  ratio  at 
which  a  maximum  in  the  toughness  of  DN-gels  was  observed 
[2].  Using  the  experimentally  measured  ^-dependent  relaxation 
time,  the  static  structures  measured  by  SANS,  and  Akcasu  et  a/.'s 
theoretical  model  [5]  for  chain  dynamics  in  a  multi-component 
polymer  solution,  we  calculated  the  effective  solvent  viscosity. 
Again,  the  effective  solvent  viscosity  for  the  P/\MPS  chain  vdthin 
the  solution  blends  exhibits  a  maximum  at  the  0.0625  mol/L 
PAMPS  concentration  (Fig.  3).  It  is  noteworthy  that  the  high  q 
portion  of  the  data  of  the  relaxation  spectrum  dictates  the  value 
of  the  effective  viscosity,  whereas  the  result  given  in  Fig.  2  arises 
from  just  the  low  q  data.  The  above  observation  suggests  that  the 
interchain  association  between  PAMPS  and  PAAm  manifests 
itself  throughout  the  entire  q  range  probed  by  NSE. 
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FIGURE  2:  A  maximum  in  the  relaxation  time  of  PAMPS  (t^J  is  observed  to 
locate  at  a  PAMPS/  PAAm  molar  ratio  of  1  to  7.  (The  line  is  a  guide  to  the 
eye.) 

In  summary,  both  the  relaxation  time  at  the  large  length  scale  and 
the  segmental  diffusion  coefficient  /  the  effective  solvent  viscosity 
at  the  short  length  scale  measured  by  NSE  in  PAMPS/PAAm 
solution  blends  are  in  accord  with  the  theological  and  mechanical 
behavior  of  the  PAMPS/PAAm  double  network  hydrogels. 
To  the  best  of  our  knowledge,  this  work  represents  one  of  the 

NCNR 


first  NSE  measurements  of  polyelectrolyte  solutions  containing 
neutral  polymer.  The  NSE  results  can  be  modeled  adequately  with 
the  mean-field  theory  of  Akcasu  ef  al.  [5]  after  it  was  modified  to 
include  polyelectrolytes. 
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FIGURE  3:  Effective  solvent  viscosity  seen  by  PAMPS  segment  within 
PAMPS/PAAm  solution  blends.  A  maximum  is  also  located  at  a  PAMPS/ 
PAAm  molar  ratio  of  1/  7.  (The  line  is  a  guide  to  the  eye.) 
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Self-assembly  is  a  promising  route  for  creation  of 
high-fidelity  and  high-density  nanostructures  but 
many  questions  regarding  the  fundamental  behavior  of 
self-assembly  remain  unanswered.  Quantifying  the  three- 
dimensional  order  of  self-assembled  nanostructures, 
especially  in  the  interior  of  thin  films,  remains  challenging. 
We  present  a  simple  and  robust  strategy  for  directing  the 
self-assembly  of  block-copolymers.  The  technique  uses 
nanoparticles  with  tunable  surface  energy  (via  ultraviolet- 
ozone  (UVO)  oxidation),  cast  onto  flat  substrates,  where  the 
spin-casting  conditions  can  be  tuned  to  control  roughness. 
These  surfaces  can  induce  block-copolymer  lamellae  to  orient 
either  parallel  or  perpendicular  to  the  substrate.  In  order  to 
quantify  the  order  within  these  thin  films,  we  used  rotational 
small-angle  neutron  scattering  in  combination  with  neutron 
reflectometry  to  reconstruct  the  sample  reciprocal  space. 
We  demonstrate  that  by  carefully  modeling  the  data,  we 
can  extract  quantitative  information  regarding  the  angular 
distribution  of  lamellae. 

Block-copolymers  are  an  archetypal  self-assembling  system, 
and  are  being  actively  pursued  for  a  variety  of  applications, 
from  microporous  membranes  to  lithographic  resists  [1]. 
Crucial  to  many  applications  is  control  of  the  orientation  of 
the  well-ordered  phases  that  develop  in  block-copolymers 
during  annealing.  Hashimoto  and  collaborators  [2] 
demonstrated  the  use  of  rough  substrates  to  bias  the 
orientation  of  block-copolymer  assembly.  Building  upon 
this  work,  we  investigated  the  use  of  nanoparticle  surface 
treatments  as  a  means  of  generating  substrates  of  controlled 
roughness.  In  particular,  spin-casting  conditions  (solution 
concentration  and  spin  speed)  can  be  used  to  control  the 
roughness  of  the  resulting  layer.  The  organic-coated  silica 
nanoparticles  (20  nm  diameter)  are  amenable  to  UVO 
oxidation  treatment,  which  converts  the  hydrophobic 
particles  to  hydrophilic.Thus,  the  average  surface  energy  can 
be  tuned  over  a  wide  range  by  UVO  exposure  time.  This  dual 
control  of  roughness  and  surface  energy  allows  us  to  select 
the  substrate  conditions  most  conducive  to  a  desired  kind  of 
self-assembled  order. 


We  investigated  the  well-studied  lamellar-forming  block- 
copolymer  polystyrene-poly(methyl  methacrylate)  (PS-PMMA). 
By  casting  a  film  thickness  gradient  of  copolymer  over  the 
nanoparticle  substrate,  we  were  able  to  investigate  the  interplay 
between  film  thickness  and  substrate  interactions.  Remarkably, 
we  discovered  that  the  surface  orientation  of  lamellae  cast  on 
the  rough  nanoparticles  substrates  oscillated  with  film  thickness 
(see  Fig.  1).  Atomic  force  microscopy  indicated  that  the 
orientation  of  lamellae  at  the  surface  was  horizontal  for  some 
thicknesses  (lamellae  sheets  lying  parallel;  flat  along  substrate 
and  air  interfaces)  and  vertical  for  others  (sheets  perpendicular; 
'standing  up' with  respect  to  the  substrate).  This  theoretically- 
implied,  but  now  experimentally  verified,  oscillation  between 
orientations  points  to  the  key  role  of  film  confinement  and 
interfacial  interaction  in  determining  the  final  order.  However, 
atomic  force  microscopy  measurements  are  fundamentally  limited 
to  measuring  surface  morphology,  and  cannot  probe  within  a  thin 
film.  For  instance,  the  fingerprint  pattern  assigned  as  Vertical 
lamellae'  in  Fig.  1  could  conceivably  arise  from  a  section  through 
a  randomly-orientated  lamellar  phase. 
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FIGURE  1 :  As  a  function  of  film  thickness,  the  orientation  oscillates  between 
vertical  and  horizontal  states.  Typical  atomic-force  micrographs  of  'vertical' 
(top)  and  'horizontal'  (bottom)  lamellae  orientations  formed  on  nanoparticle- 
roughened  substrates. 

Neutron  scattering  measurements  can  provide  the  missing 
information  about  the  internal  arrangement  of  nanostructured 
thin  films.  We  used  semi-deuterated  block-copolymers,  where 
one  block  has  hydrogens  replaced  with  deuteriums,  which  scatter 
neutrons  strongly  and  yield  characteristic  patterns  in  a  small- 
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angle  neutron  scattering  (SANS)  geometry.  We  further  employed 
the  'rotational  SANS' variant,  whereby  we  acquire  a  scattering 
image  across  a  wide  range  of  sample  rotation  angles.  Each  angle 
measures  a  'sUce'  through  the  three-dimensional  reciprocal  space 
of  the  sample,  and  the  angular  sweep  thus  reconstructs  the  full 
reciprocal  space.  Since  the  scattering  pattern  at  each  angle  arises 
from  periodic  order  along  a  particular  plane  in  the  sample,  the 
reconstructed  reciprocal  space  in  principle  contains  information 
about  the  order  in  every  possible  orientation.  Figure  2  shows  a 
Q  -Q.  sUce  through  the  reconstructed  sample  reciprocal  space, 
for  three  representative  thin  films.  The  z-axis  points  along  the 
sample  normal,  with  the  sample  lying  in  the  xy-plane.  Scattering 
intensity  along  the  x-axis  thus  originates  from  in-plane  periodic 
structures,  whereas  intensity  along  the  z-axis  originates  from 
nanostructures  whose  periodicity  is  along  the  film  normal.  The 
Q  -Q^  slices  provide  information  about  the  range  of  lamellae 
orientations  from  vertical  to  horizontal  (note  that  the  j-direction 
is  essentially  identical  to  the  x-direction  since  the  samples  had  no 
preferred  in-plane  direction).  The  peak  at      =  0.16  nm"^  occurs 
due  to  vertical  lamellae,  whereas  the  scattering  intensity  near  the 

axis  arises  from  horizontal  lamellae.  Importantly,  the  data 
allow  us  to  rule  out  a  completely  random  internal  arrangement 
of  lamellae,  since  this  would  lead  to  uniform  isotropic  scattering 
for  all  5  =  ^joz+Q.'  =  0.16  nm"\  It  is  also  clear  that  the  extent 
of  vertical  order  is  affected  by  fUm  thickness.  A  smaU  variation 
in  the  thickness  of  the  film  causes  the  vertical  peak  to  nearly 
disappear. 
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FIGURE  2:  Rotational  SAMS  reciprocal-space  reconstructions  for  three 
different  film  thicknesses  (top  row)  show  varying  intensity  of  the  vertical 
peak.  The  associated  best-fit  models  (bottom  row)  provide  insight  into  the 
angular  distribution  of  the  lamellae. 

We  can  quantify  the  order  by  fitting  the  data  to  a  model  for 
the  scattering  from  lamellar  structures.  We  extend  a  literature 
model  [3],  where  the  angular  distribution  can  be  decoupled 
from  the  scattering  of  a  lamellar  'stack',  and  explicitly  include 
the  instrumental  broadening  [4].  By  minimizing  the  sum  of 
the  squared-residuals  between  the  model  and  the  experimental 
data,  we  can  determine  the  angular  distribution  of  lamellae  that 


distributions  for  the  three  samples  shown  in  Fig.  2,  along 
with  schematics  of  the  inferred  ordering.  Variations  in  film 
thickness  change  the  proportion  between  vertical  and  horizontal 
ordering.  The  rotational  SANS  technique  cannot  acquire 
data  at  near-grazing  angles  due  to  strong  reflection  from  the 
substrate  interface.  However,  neutron  reflectometry  is  perfectly 
complementary  in  this  regard,  since  it  is  sensitive  to  order  in 
the  film  normal  direction  (the  horizontal  lamellae  stacking). 
By  using  neutron  reflectometry,  we  were  able  to  corroborate 
the  best-fit  models  obtained  from  the  rotational  SANS  data: 
for  some  thicknesses,  the  horizontal  lamellae  form,  whereas  for 
other  thicknesses  their  formation  is  suppressed.  Moreover,  the 
reflectometry  makes  it  clear  that  the  formation  of  horizontal 
lamellae  is  disrupted  near  the  rough  substrate  interface. 
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FIGURE  3:  Best-fit  orientational  distribution  of  lamellae  for  the  samples 
shown  in  Figure  2  (top  row).  The  vertical  peaks  (0°)  are  broader  than  the 
horizontal  peaks  (±  90°).  Schematics  of  orientational  order  (bottom  row) 
can  be  constructed  from  the  distribution. 

This  analysis  of  the  three-dimensional  order  helps  elucidate  the 
origin  of  the  orientational  changes.  Block-copolymer  lamellae 
on  flat  substrates  always  form  horizontal  sheets  due  to  the  strong 
wetting  at  the  substrate  interface.  By  introducing  roughness, 
this  wetting  is  disrupting,  thereby  making  the  vertical  state 
energetically  accessible.  The  oscillation  as  a  function  of  thickness 
arises  due  to  commensurability.  When  the  film  thickness  is 
equal  to  a  whole  number  of  lamellae  sheets,  this  arrangement 
is  energetically  favored.  However,  when  the  film  thickness 
is  incommensurate,  and  a  whole  number  of  layers  cannot  be 
formed,  the  horizontal  state  is  unstable  and  the  vertical  state 
is  instead  the  energy  minimum.  Combining  the  atomic  force 
microscopy,  neutron  reflectometry,  and  rotational  SANS  thus 
gives  a  remarkably  thorough  picture  of  the  film  structure.  By 
modeling  the  rotational  SANS  data,  we  can  obtain  quantitative 
angular  distributions  for  the  nanostructures,  which  are  found 
to  be  consistent  with  the  atomic  force  microscopy  and  neutron 
reflectometry  data.  This  represents  a  general  strategy  that  can  be 
appUed  to  a  wide  variety  of  nanostructures  that  assemble  in  thin 
films. 
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hear  banding  is  a  flow-induced  phenomenon  observed 
in  a  wide  variety  of  materials,  including  highly  entangled 


polymers,  self-assembled  surfactants,  colloidal  suspensions, 
and  pastes.  Its  theological  and  macroscopic  signatures  have 
been  well-studied  for  over  a  decade,  yet  surprisingly  little  is 
known  about  the  underlying  microstructural  mechanism(s) 
that  give  rise  to  shear  banding  in  soft  matter.  In  this  report, 
we  demonstrate  the  capabilities  of  a  new  shear  cell  for 
spatially-resolved  small  angle  neutron  scattering  (SANS) 
measurements  in  the  flow-gradient  (1-2)  plane  to  elucidate 
the  microstructural  mechanisms  by  which  shear  banding 
occurs  in  a  model  wormlike  micellar  surfactant  solution. 

Perhaps  the  most  well-studied  fluids  that  exhibit  shear 
banding  are  viscoelastic  wormlike  micelle  (WLM)  solutions, 
comprised  of  long,  entangled  threadlike  aggregates  of 
amphiphilic  molecules  [1].  Shear  banding  in  these  systems 
was  originally  observed  visually  as  birefringent  bands  near  the 
rotating  wall  of  a  Couette  flow  geometry  [2].  The  signature  of 
banding  is  segregation  of  the  flow  field  into  two  fluid  layers, 
one  at  a  high  shear  rate  and  one  at  a  low  shear  rate,  that  span 
the  flow  geometry.  These  bands  coincide  with  a  stress  plateau 
in  steady  state  shear  rheological  measurements,  (Fig-  1). 


FIGURE  1:  Diagram  of  the  1-2  plane  flow-SANS  shear  cell.  A  plot  of 
the  neutron  transmission  demonstrates  the  gap  resolution  obtainable. 
The  illustration  shows  the  ability  to  probe  gap-resolved  structure  in  the 
two  different  shear  bands. 


Rheological  theories  that  explain  this  behavior  assume  a  non- 
monotonic constitutive  stress-rate  relationship,  similar  to  a  van 
der  Waals  loop  equation  of  state  for  first  order  phase  transitions 
[3].  This  has  led  many  investigators  to  propose  that  shear  banding 
coincides  with  an  underlying  shear-induced  phase  transition.  For 
example,  the  cationic  surfactant  CTAB  in  D^O  is  known  to  shear 
band  in  the  vicinity  of  an  equilibrium  isotropic-nematic  (I-N) 
transition  [4] .  However,  rigorous  validation  of  this  mechanism 
remains  elusive,  due  to  challenges  in  measuring  separately  the 
surfactant  aggregate  microstructure  in  both  the  high-shear  and 
low-shear  bands. 

Recently,  we  have  developed  the  ability  to  measure  spatially- 
resolved  microstructure  under  shear  via  SANS  by  using  a  short 
gap  Couette  cell  and  coUimating  the  incident  neutron  beam  down 
the  gap  (vorticity  axis).  A  slit  aperture  enables  collecting  SANS 
from  slices  in  the  flow-gradient  (1-2)  plane  of  shear  [5].  These 
1-2  plane  flow-SANS  measurements  provide  several  capabilities 
lacking  in  other  techniques  to  measure  structure  under  shear, 
namely:  (1)  the  size  ranges  accessible  by  SANS  cover  the  relevant 
microstructural  scales  of  WLM  solutions  (e.g.,  miceUe  radius, 
persistence  length,  mesh  size,  etc.);  (2)  measurements  in  the  1-2 
plane  allow  quantification  of  both  the  segmental  orientation  and 
degree  of  alignment,  whereas  more  common  measurements  in  the 
1-3  plane  measure  only  a  projection  of  the  segmental  alignment 
[4];  (3)  the  slit  aperture  provides  gap-resolved  measurements  of  the 
structure,  which  allows  discrimination  between  the  high-shear  and 
low-shear  bands.  Figure  1  shows  a  diagram  of  the  shear  cell  and 
demonstration  of  these  capabilities. 

In  this  highlight,  we  report  results  for  a  model  shear  banding  WLM 
of  the  cationic  surfactant  CTAB  at  0.49  mol/L  (490  mM)  and 
32  °C,  for  which  the  rheology  and  shear  banding  are  well- 
characterized  [6].  Flow-SANS  measurements  are  performed  at 
seven  positions  across  the  1.35  mm  Couette  gap  using  a  0.1  mm  slit. 
Figure  2  displays  a  visual  summary  of  the  results,  where  the  intensity 
ring  is  a  correlation  peak  due  to  segment-segment  interactions. 
Anisotropy  in  this  ring  indicates  segmental  flow  alignment,  with 
high  alignment  typical  for  a  nematic  phase.  A  significant  difference 
in  scattering  anisotropy  is  observed  between  positions  in  the  low- 
shear  and  high-shear  bands  (also  shown  in  Fig  1.) 
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FIGURE  2:  Two-dimensional  scattering  plots  for  nominal  applied  shear 
rates  and  normalized  gap  positions  spanning  the  shear  banding  transition 
for  the  CTAB  sample.  The  black  line  indicates  the  measured  location  of  the 
interface  between  the  high-shear  and  low-shear  bands. 

The  average  segmental  orientation,     of  the  micelles  relative  to  the 
flow  direction  is  defined  in  Fig.  1,  where  0°  <     <  45°.  Similarly, 
the  net  segmental  alignment  is  given  by  the  alignment  factor, 
which  characterizes  the  orientational  order  in  the  fluid  [6].  As  the 
average  shear  rate  at  each  gap  position  is  known  from  independent 
velocimetry  measurements  [6],  all  of  the  data  can  be  plotted  as 
master  curves  of  (p^  and     versus  the  local  shear  rate  (Fig.  3).  Doing 
so  demonstrates  a  clear  transition  at  ip^  ~  10°  3ind  Aj.~  0.15,  where 
a  jump  in  orientation  and  alignment  occurs  between  the  low-shear 
and  high-shear  bands.  Predictions  from  a  constitutive  model  that 
couples  the  fluid's  theology  and  micellar  orientation  order  are  in 
excellent  quantitative  agreement  in  both  low  and  high  shear  states. 


local  shear  rate  (s'  ) 


FIGURE  3:  Orientation  angle  [closed]  and  alignment  factor  (open)  versus 
local  shear  rate  measured  by  velocimetry  in  the  high-shear  (red]  and 
low-shear  (blue]  bands.  Half  filled  symbols  are  for  measurements  with 
contributions  from  both  bands.  Data  are  augmented  at  low  shear  rates 
using  flow-birefringence  measurements  [black].  Lines  give  predictions  from 
the  Giesekus  model  (solid). 

The  high  degree  of  microstructural  order  observed  at  the  highest 
shear  rates  is  consistent  with  a  flow-aligned  nematic  order, 


significantly  less  order  is  observed  for  non-banding  WLM  solutions 
at  comparable  shear  rates  [6]. This  is  further  confirmed  through 
sector-averaged  plots  of  the  intensity  in  the  flow-aligned  direction 
(Fig.  4),  which  show  that  the  low  shear  band  has  a  nearly  identical 
structure  to  the  fluid  at  rest,  whereas  the  high-shear  band  shows 
an  increase  in  sharpness  and  location  of  the  structure  peak  similar 
to  what  is  observed  for  an  equilibrium  nematic  phase  (shown  for 
reference.)  This  confirms  the  hypothesis  that  shear  banding  results 
from  a  shear-induced  I-N  transition  for  CTAB.  This  is  in  contrast 
to  some  other  surfactant  systems  that  show  a  transition  to  a  biphasic 
network  structure  [5].  Nonetheless,  in  both  systems,  the  underlying 
thermodynamic  phase  behavior  plays  a  critical  role  in  determining 
how  and  when  a  fluid  wiU  exhibit  shear  banding. 
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FIGURE  4:  Sector-averaged  intensity  in  the  flow  direction  for  the  shear 
rates  indicated.  Open  symbols  show  structure  in  the  high-shear  [red]  and 
low-shear  (blue)  bands.  The  sharpening  of  the  peak  at  high  rates  is  similar 
to  what  is  observed  for  a  flow-aligned  nematic  phase  at  rest  (black). 

In  conclusion,  spatially-resolved  1-2  plane  flow- SANS 
measurements  enable  direct  measurement  of  the  mesophase 
microstructure  and  orientational  order  of  WLMs  in  each  shear 
band  under  flow.  The  results  provide  critical  information  that  is 
being  used  to  test  and  refine  microstructure-based  models  of  shear 
banding  [6].  This  technique  is  unique  in  its  ability  to  resolve  local 
structure  in  non-homogeneous  flows  and  as  such  is  a  powerful 
tool  for  quantitative  interrogation  of  the  nonlinear  behavior  of  soft 
matter  under  shear. 

References 

[1]  H.  Rehage  and  H.  Hoffmann,  Molec.  Physics  74,  933  (1991). 

[2]  J.  P.  Decruppe,  R.  Cressely,  R.  Makhloufi,  E.  Cappelaere,  Coll.  and  Poly.  Sci. 

273,346  (1995). 

[3]  J.  K.  G.  Dhont  and  W.J.  Briels,  Rheol.  Acta  47, 257  (2008). 

[4]  J.  F.  Berret,  in  Molecular  Gels,  ed.  R.  G.  Weiss  and  P.  Terech  (Springer, 

Netherlands,  2005),  p.  663. 

[5]  M.  W.  Liberatore,  F.  Nettesheim,  N.  J.  Wagner,  L.  Porcar,  Phys.  Rev.  E  73, 
020504  (2006). 

[6]  M.  E.  Helgeson,  E.  W.  Kaler,  N.J.  Wagner,  R  A. Vasquez,  L.  P  Cook, J.  Rheol, 
Submitted  (2008). 


NCNR  2008  Accomplishments  and  Opportunities 


43 


A  Facile  Route  for  Creating  Reverse  Vesicles  in  Organic  Liquids 


S.-H.  Tung^  and  S.  R.  Raghavan 


Vesicles  are  nanoscale  containers  formed  in  water  and 
enclosed  by  a  bilayer  of  lipid  molecules.  They  are  of 
great  interest  in  drug  delivery  and  nanomedicine  due  to  their 
ability  to  encapsulate  hydrophilic  solutes  in  their  aqueous 
interior.  Much  like  such  "normal"  vesicles  in  v\^ater,  one  can 
also  envision  their  counterparts  in  organic,  nonpolar  solvents 
("oils"),  which  could  be  then  termed  "reverse"  vesicles.  If 
reverse  vesicles  could  be  formed  easily  and  reproducibly, 
it  is  easy  to  imagine  that  they  too  could  find  numerous 
applications  such  as  encapsulation  and  controlled  delivery 
of  hydrophobic  solutes  [1,2].  Here,  we  highlight  our  recent 
work  [3]  on  a  simple  route  to  forming  stable  unilamellar 
reverse  vesicles  in  nonpolar  solvents  such  as  cyclohexane.  Our 
method  also  offers  a  general  framework  for  tuning  reverse 
aggregate  geometry  from  reverse  spherical  micelles  to  reverse 
cylindrical  micelles  and  finally  to  reverse  vesicles. 

Figure  1  schematically  contrasts  normal  and  reverse  vesicles 
and  also  describes  our  method  for  forming  the  latter  [3]. 
As  is  well-known,  normal  vesicles  are  formed  in  water 
by  the  self-assembly  of  long-tailed  phospholipids  such  as 
L-a-phosphatidylcholine  (lecithin),  a  natural  two-tailed  lipid 
with  an  average  tail  length  of  17  carbons  and  an  unsaturation 
in  one  of  the  tails.  In  the  figure,  lecithin  is  shown  with  a 
blue  head  (hydrophilic)  and  two  red  tails  (hydrophobic). 
Note  the  structure  of  the  bilayer  in  normal  vesicles:  here, 
the  hydrophobic  tails  are  in  the  interior  of  the  bilayer  so 
that  they  are  shielded  from  water.  To  form  reverse  vesicles, 
we  combine  lecithin  with  a  short-chain  phospholipid, 
l,2-dibutyroyl-^77-glycero-3  phosphocholine  (C^-lecithin) 
-  a  lipid  with  two  four-carbon  saturated  tails.  When  added 


FIGURE  1:  Molecular 
structures  of  lecithin  and 
C^-lecithln,  and  schematics 
of  their  assemblies.  Lecithin 
alone  when  added  to  water 
forms  normal  vesicles,  which 
are  illustrated  schematically 
on  the  left.  The  present  study 
shows  that  mixtures  of  lecithin 
and  C^-lecithin  added  into 
oil  form  reverse  vesicles,  as 
depicted  on  the  right. 


"normal" 
vesicles 
in  water 


to  an  oil  (cyclohexane)  at  a  high  molar  ratio      of  C^-lecithin  to 
lecithin,  the  lipids  assembled  into  reverse  vesicles.  Note  that  the 
reverse  vesicles  have  an  oil  core,  and  their  shell  is  a  reverse  bilayer 
in  which  the  hydrophobic  portions  of  the  lipids  are  exposed  to 
the  oil  both  in  the  core  and  in  the  exterior. 


FIGURE  2:  Photograph  of  sample  vial  and  i  civi  n.ici  ogi  ap,i  o.  o.. .  es 
present  in  a  C^-lecithin  +  lecithin  sample  at     =  2.6  (total  lipids  =  20 
mmoi/L,  [NaCI]  =  3.5  mmol/L). 

Evidence  for  the  presence  of  reverse  vesicles  in  the  mixed  lipid 
solutions  was  obtained  by  a  combination  of  dynamic  light 
scattering  (DLS),  small-angle-neutron  scattering  (SANS) 
and  transmission  electron  microscopy  (TEM)  [3]. The  reverse 
vesicle  sample  (7?^  =  2.6)  has  a  clear,  bluish  tinge,  as  seen  from 
the  photograph  in  Fig.  2.  The  bluish  color  is  a  manifestation  of 
the  Tyndall  effect,  i.e.,  it  is  due  to  large  scatterers  in  solution, 
and  is  a  general  feature  of  vesicles  in  water  as  well.  DLS  gives  a 
hydrodynamic  radius  of  130  nm  for  the  structures  in  the  sample. 
Fig.  2  shows  a  TEM  micrograph  of  the  sample:  here  ammonium 
molybdate  was  added  as  a  positive  stain  that  would  bind  with 
the  headgroups  of  lipids  and  thus  clearly  reveal  the  bilayers.The 


+  "  "reverse 

vesicles 


water 


Lecithin 


C.  -Lecithin      in  oil 
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TEM  image  shows  a  number  of  spherical  structures  with  distinct 
shells,  much  like  micrographs  ot  normal  vesicles.  The  diameters 
of  these  structures  range  from  ~  300  nm  to  1  [.im.  While  these 
structures  are  larger  than  those  expected  from  DLS,  we  believe 
this  disagreement  is  due  to  the  coUapse  or  fusion  of  the  reverse 
vesicles  as  the  solvent  is  evaporated  from  the  TEM  grid. 

To  further  elucidate  the  structure,  we  chose  SANS,  which  was 
conducted  on  the  NG-3  beamline  at  the  NCNR  [3].  Samples 
were  made  in  deuterated  cyclohexane  to  achieve  the  needed 
contrast  between  scatterers  and  solvent.  SANS  spectra  (intensity 
/  vs  wave  vector  q)  for  20  mmol/L  lipid  solutions  are  shown  in 
Fig.  3  at  varying  R^.  The  data  for  Rq  =  0  and  R^^  =  0.5  are  both 
asymptotic  to  a  plateau  at  low  q  and  essentially  correspond  to 
micelles.  In  contrast,  there  is  no  plateau  at  R^  =  2.6,  with  the 
data  showing  a      decay  of  the  intensity  at  low  q.  Such  decay 
is  a  signature  of  scattering  from  vesicles.  The  SANS  data  thus 
indicate  a  transition  from  reverse  micelles  to  reverse  vesicles  with 
increasing  R^.  We  modeled  the  SANS  data  using  appropriate 
form  factors  and  the  fits  are  shown  as  solid  curves  through  the 
data.  For  lecithin  in  cyclohexane  {R^  =  0),  the  reverse  micelles 
are  ellipsoids  with  radii  of  2.2  nm  and  3.0  nm,  respectively,  for 
their  minor  and  major  axes.  Upon  the  addition  of  small  amounts 
of  C^-lecithin,  the  micelles  grow  axiaUy,  and  for  R^^  =  0.5,  they 
can  be  modeled  as  rigid  cylinders  (radius  of  2.2  nm  and  length 
of  19.8  nm).  For  R^  =  2.6,  the  data  are  fit  to  the  polydisperse 
unilamellar  vesicle  model.  An  average  reverse  vesicle  diameter 
of  234  nm  is  obtained,  along  with  a  bilayer  thickness  of  3.7  nm, 
and  a  polydispersity  of  0.22.  In  sum,  the  SANS  data  confirm 
the  evolution  of  self-assembled  reverse  structures  from  nearly 
spherical  micelles  to  cylindrical  micelles  to  vesicles  with 
increasing  molar  ratio  of  C^-lecithin  to  lecithin. 


FIGURE  3:  SANS  spectra  from  lipid  mixtures  in  deuterated  cyclohexane  for 
different  C^-lecithin:lecithin  molar  ratios  R^.  The  samples  contain  20  mmol/L 
total  lipid  and  3.5  mmol/L  NaCL 


Why  does  the  addition  of  C^-lecithin  induce  the  formation  of 
reverse  vesicles?  It  is  known  that  the  shape  of  self-assembled 
structures  formed  by  amphiphiles  is  governed  by  their  geometry 
[4].  This  connection  is  usually  expressed  in  terms  of  the 
critical  packing  parameter^  =  where  a^^.^^  and  a^^  are  the 

cross-sectional  areas  of  the  amphiphiles  tail  and  headgroup, 
respectively.  In  the  case  of  organic  solvents,  the  formation  of 
reverse  micelles  requires  a  packing  parameter  p  much  larger  than 
1,  and  spherical  reverse  micelles  correspond  to  an  inverse  cone 
shape.  For  these  spheres  to  transform  into  cylinders  or  vesicles, 
the  packing  parameter  p  has  to  decrease.  We  believe  that  such  a 
decrease  is  caused  by  the  binding  of  C^-lecithin  to  the  lecithin 
headgroups.  That  is,  since  C^-lecithin  has  the  same  headgroup  as 
lecithin  but  much  shorter  tails,  it  increases  the  headgroup  area  a^^, 
while  the  overall  tail  area  remains  about  the  same.  The  net  effect 
is  to  decrease  ^,  and  at  high  values  of  R^,  the  effective  geometry  is 
altered  to  a  cylinder  shape  {p  ~  1),  thus  leading  to  reverse  vesicles. 
Thus,  our  results  can  be  placed  within  a  general  framework  (or 
rules  of  thumb)  for  tuning  reverse  aggregate  geometry. 

One  aspect  that  we  have  not  discussed  thus  far  is  that,  in 
addition  to  the  two  lipids,  we  must  also  add  a  trace  amount  of 
salt  {e.g..,  NaCl)  to  ensure  the  stability  of  the  reverse  vesicles. 
Otherwise,  the  reverse  vesicles  become  unstable  and  phase- 
separate.  However,  in  the  presence  of  a  few  mmol/L  of  salt,  the 
reverse  vesicles  are  very  stable  and  robust,  retaining  their  size  and 
structure  over  a  period  of  many  months.  Why  is  salt  necessary? 
We  believe  that  the  salt  ions  (e.g.,  Na*  and  CI")  bind  with  the 
negative  and  positive  charges  on  the  lipid  headgroups  through 
short-range  electrostatic  interactions  and  thereby  serve  as  a  "glue" 
between  the  lipids  [3]. 

In  conclusion,  we  have  demonstrated  a  new  approach  for  creating 
reverse  unilamellar  vesicles  in  oils  -  it  involves  mixing  short  and 
long-chain  lecithins  along  with  a  small  amount  of  a  salt  such 
as  NaCl.  The  components  involved  are  commercially  available 
and  the  reverse  vesicles  thus  formed  are  stable  and  robust. 
Future  experiments  will  examine  the  ability  of  reverse  vesicles  to 
encapsulate  solutes  in  their  interior.  We  hope  that  our  approach 
wiU  stimulate  new  interest  into  this  fascinating  class  of  self- 
assembled  structures. 
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Dispersions  of  Carbon  Nanotubes  in  Polymers:  Hierarchical 
Structure  and  Viscoelastic  Properties 

T.  Chatterjee\  A.  Jackson^^  and  R.  Krishnamoorti^ 


Nanoparticles  demonstrate  a  number  of  advantages 
over  traditional  inorganic  fillers  used  in  the  polymer 
industry  as  a  result  of  their  huge  surface  area  per  mass 
and  often  very  high  aspect  ratio.  In  particular,  dramatic 
improvements  in  tensile  strength  and,  electrical  and  thermal 
conductivity  have  been  observed  at  much  lower  mass  fraction 
than  traditional  inorganic  fillers.  Understanding  the  origins 
of  these  structural  improvements  and  associated  theological 
phenomena  is  important  for  optimizing  the  properties  of 
polymer  nanocomposites  and  their  processing. 

At  relatively  high  concentrations,  dispersed  nanoparticles 
form  superstructures  and  networks  which  can  dominate  the 
mechanical  and  theological  properties  of  the  nanocomposites 
in  a  way  that  is  only  weakly  dependent  on  the  nature  of  the 
nanoparticles  or  of  the  polymer  matrix  [1,2].  Specifically, 
setni-dilute  dispersions  of  spherical  silica  nanoparticles,  fractal 
carbon  black,  rod-like  nanotubes,  and  disk-like  layered- 
silicates  show  remarkable  similarity  in  their  rheological 
properties.  These  dispersions  have  a  concentration  of 
nanoparticles  (/>)  significantly  larger  than  the  percolation 
threshold  or  jamming  transition  {p(^.  Such  dispersions  in 
oligomeric  or  low  molecular  weight  polymers  demonstrate 
linear  viscoelastic  properties  that  are  time-temperature- 
composition  superposable  and  are  gel-like. 

The  low-frequency  modulus  of  these  dispersions  displays 
a  concentration  scaling  ~  {p  -^  )\  which  corresponds  to 
the  network  properties,  independent  of  the  chemical  and 
topological  identity  of  the  primary  nanoparticles.  Using 
small  and  ultrasmall  angle  neutron  scattering  (SANS 
and  USANS  respectively),  and  focusing  on  the  nature  of 
the  network  structure  formed  by  dispersed  single  walled 
carbon  nanotubes  (SWNTs)  in  a  polymer  matrix,  we  have 
identified  the  relationship  between  the  network  structure 
and  rheological  properties  of  such  dispersions  [1],  which  we 
briefly  describe  in  this  highlight. 

Dispersions  of  SWNTs  in  three  matrices,  namely  the  epoxy 
resin  matrices  bisphenol  A  (BA),  bisphenol  F  (BF)  and  the 
higher  molecular  mass  poly(ethylene  oxide)  (PEO),were 


examined.  The  viscosities  for  these  matrices  are  comparable  and 
range  from  0.7  Pa  s  to  2.2  Pa  s.  The  concentrations  of  SWNTs 
in  each  of  the  dispersions  were  chosen  to  ensure  correspondence 
with  the  semi-dilute  region  (2  <  p/pc  <  15),  a  range  of  practical 
importance  as  well  as  largely  unexplored  in  terms  of  structure. 
The  sample  preparation,  dispersion  state,  and  linear  rheological 
behavior  of  some  of  these  systems  have  been  discussed  elsewhere 
[1,3]. 

On  the  basis  of  previous  studies  on  other  nanoparticulate 
polymer  systems,  we  hypothesize  that  on  a  macroscopic  length 
scale  a  matrix-spanning  network,  consisting  of  aggregated  clusters 
(floes)  of  dispersed  SWNTs,  is  established.  The  scattering  is 
dominated  by  the  large  floes,  characterized  by  the  cluster/floc 
size  {R).  Inside  the  floe,  individual  or  small  bundles  of  tubes 
(depending  on  the  dispersion  state)  overlap  and  the  average 
distance  between  two  adjacent  contacts  is  the  network  mesh  size 
(l).  A  schematic  diagram  of  this  hierarchical  structure  is  presented 
in  Fig;.  1. 


Mesh  size,  ^ 


Floe  size,  R 


FIGURE  1:  A  schematic  diagram  of  the  hierarchical  network  structure 
showing  two  different  length  scales;  we  note  that  there  can  be  quite  a  range 
of  mesh  sizes,  C  and  floe  sizes,  R. 

The  combination  of  SANS  and  USANS  offers  a  broad  range  of 
scattering  vector,  q,  ideal  for  the  investigation  of  the  hierarchy 
of  length  scales  that  are  expected  from  dispersed  SWTMTs  at 
concentrations  in  excess  of  the  percolation  threshold.  The  neutron 
scattering  measurements  were  performed  using  the  USANS 
instrument  at  BT-5  and  the  SANS  diff'ractometer  at  NG-7 
at  the  NCNR.  The  scattered  intensity,  when  scaled  by p{l-p) 
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(expected  for  disordered  systems),  superpose  to  form  a  master 
curve  (Fig.  2),  indicating  the  conservation  of  the  structure  over 
the  concentration  and  length  scales  (~  ^  ')  examined.  Scattering 
from  a  hierarchical  fractal  can  be  described  by  a  unified  equation  with 
multiple  structural  levels  [4] . 
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In  this  description,     ^  gives  a  measure  ot  the  floe  size  R,  and  R^ , 
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exponent corresponds  to  the  average  mass  fractal  dimension  i^of 
the  floe,  and  y,  to  the  fractal  dimension  d   .  of  the  nanotube  bundles. 
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FIGURE  2:  Concentration-effect  corrected  scattering  data  lead  to  a  master 
curve  indicating  a  hierarchical  network  structure  over  a  wide  range  of 
length  scales.  The  lines  show  slopes  in  two  regions. 

The  dfof  the  SWNTs  network  is  found  to  be  a  function  of  the 
surrounding  matrix.  Matrices  with  a  higher  value  of  pc  exhibit  a 
higher  value  of  ^^and  consequently  denser  aggregation  within 
the  floes.  The  other  relevant  length  scale  of  the  network,  the 
mesh  size  (l)  shows  a  dependence  on  nanotube  loading  (Fig.  3). 
Increased  loading  corresponds  to  decreased  I  and  exhibits  a  power 
law  (X~  p     dependence  with  a  values  ranging  between  0.3 
(±0.03)  and  0.4  (±0.05)  for  the  systems  studied  here.  These  weak 
crowding  values  of  a  are  consistent  with  diffusion-Umited  cluster 
formation  [5]. 

Interestingly,  the  average  floe  size  R  is  found  to  be  ~  4  ^m  and 
almost  independent  of  the  surrounding  matrix  and  particle 
concentration.  The  mesh  fractal  dimension,  d^^^^,  is  found  to 
be  ~  2.8  ±0.1  for  all  matrices,  consistent  with  the  formation  of 
nanotube  bundles  as  the  basic  unit  of  the  SWNT  network  within 
the  floes. 
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FIGURE  3;  Concentration  dependence  of  the  network  mesh  size  (0  and 
in  different  matrices.  The  scaling  of  Cis  consistent  with  diffusion  limited 
floe  formation  while  the  nearly  linear  dependence  of  G^,  roughly  measuring 
floe  number,  suggests  that  the  network  primarily  grows  through  new  floe 
formation.  Error  bars  are  ±1  standard  deviation. 

These  structural  results  help  in  understanding  the  viscoelastic 
properties  of  these  materials.  The  network  elastic  strength  might 
be  expected  to  be  proportional  to  the  number  of  stress-bearing 
junctions  between  nanotube  bundles.  But  the  high  values  of 
concentration  scaling  of  the  modulus  observed  for  the  three 
nanocomposite  series  (scaling  exponents  S  of  3.8  ±  0.3,  2.5  ±  0.2 
and  2.5  ±  0.4  for  PEO,  BA  and  BF  systems,  respectively)  are 
not  easily  reconciled  in  the  context  of  the  composition-invariant 
floe  size  and  the  weak  crowding  of  the  mesh  within  a  floe.  On 
the  other  hand,  Gj  demonstrates  a  nearly  linear  dependence 
with  nanotube  concentration  (~  ^'■^*°^).  Since  for  a  3D  isotropic 
system  the  number  of  inter-floc  contacts  can  be  estimated  to  go 
as  (~  ^'•^-°-')-\  we  suggest  instead  that  the  interactions  between 
floes  (either  directly  or  mediated  by  polymer),  control  the  scaling 
dependence  §  of  the  elastic  strength  of  the  network. 

This  study  shows  [6]  that  elastic  properties  in  a  variety  of 
polymer  nanocomposites  may  be  primarily  controlled  by  inter- 
floc  interactions.  The  specific  network  elastic  strength  can 
thus  be  used  as  a  quantitative  tool  to  understand  and  compare 
nanoparticle  dispersion  in  different  matrices  irrespective  of  the 
matrix,  chemical  and  topological  identity  of  the  nanoparticle,  and 
nanoparticle  loading. 
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Curved  Neutron  Guides:  Phase  Space  Tailorin 

J.  C  Cook^  and  D.  F.  R.  Mildner^ 


Spectrometers  that  demand  high  Q-resolution  or  have  long 
guide-free  pre-sample  drift  paths,  such  as  SANS  or  NSE 
instruments,  require  low  divergence  beams.  To  date  at  the  NCNR, 
cold  neutron  beams  have  been  delivered  to  these  instruments 
either  by  straight  guides  v^fith  a  cooled  crystal  filter,  or  by  a  neutron 
optical  filter  -  a  guide  arrangement  with  one  or  more  kinks  that 
excludes  direct  lines  of  sight  between  the  source  and  the  instrument. 
The  crystal  filter  is  typically  a  liquid  nitrogen-cooled,  Be/Bi 
combination,  with  the  Bi  required  because  of  the  gamma  sensitivity 
of  area  detectors.  The  main  disadvantages  of  the  crystal  filter  are 
that  the  minimum  transmitted  wavelength  is  set  by  the  Be  cutoff  at 
about  4  A,  the  filter  materials  attenuate  the  useful  1  >  4  A  neutrons, 
the  guide  break  for  the  filter  cryostat  generates  additional  losses,  and 
the  system  is  non-passive.  Some  of  these  problems  are  overcome 
using  an  optical  filter;  however,  the  transmission  only  tends  to  that 
of  an  unfiltered  straight  guide  ideally  in  the  long  wavelength  limit, 
and  spatial- angular  asymmetries  in  the  beam  may  be  present.  We 
have  developed  a  simple  "phase-space  tailoring"  alternative  whereby 
the  guide  transmission  can  ideally  match  that  of  an  unfiltered 
straight  guide  of  equivalent  length  with  equivalent  beam  divergence 
for  an  extended  wavelength  range  (often  the  entire  operating  range 
of  the  instrument)  [1]. 

Consider  first  a  two-dimensional  unfiltered  straight  guide  of  length 
L,  width  W,  and  side  coating  critical  angle  5 ,  delivering  a  beam  to 
an  instrument  whose  minimum  operating  wavelength  is  1^^^^^^.  For 
idealized  perfect  reflectivity  and  guide  entrance  illumination,  the 
emerging  beam  intensity  is  uniform  across  the  guide  width  and 
uniform  in  angle  between  the  limits  ±5 ,  with  unit  transmission  for 
the  entire  instrument  wavelength  range,  provided  that  L  >  W/?ck^.^, 
with  5  (2)  «  xl  where  a:  is  a  reflective  coating-dependent  constant. 
This  situation  would  be  ideal  for  the  instrument  except  that 
excessive  fast  neutron  and  gamma  transmission  usually  prohibits 
the  use  of  an  unfiltered  straight  guide.  We  propose,  therefore, 
to  precede  the  straight  guide  (whose  critical  angle  is  henceforth 
denoted  by  5    =  x^^X)  with  a  circular-arc  curved  guide  of  radius 
_p,  width  H^,  length  Z  ,  and  inner-  and  outer-curve  surface  coating 
critical  angles,  ^^.^  and  5     respectively. The  curved  guide  must 
function  so  as  to  (a)  remove  beam  contaminants  and  (b)  uniformly 
illuminate  the  straight  guide  within  the  required  limits  ±S  for 
all  wavelengths  greater  than  A    .  Satisfying  both  conditions  is 
achieved  by  eliminating  Unes  of  sight  within  the  curved  section, 
which  requires  that  L_  >  yliSWj} ).  Condition  (b)  cannot  be  achieved 


X  =  -w/2 


unless  a  beam  entering  at  the  inner  surface  just  reaches  the  inner 
surface  again  upon  exiting,  which  happens  when  9  >  yJ(2W/p)  =  f,  the 
"characteristic  angle".  The  minimum  wavelength,  1,  for  which  this  is 
possible  occurs  when  1^  =  f-^  /^„„,>  where  x^^^  is  the  critical  angle  per  unit 
wavelength  of  the  outer  surface  coating.  If  f{x)  denotes  the  angle  of  a 
trajectory  with  respect  to  the  guide  axis  at  a  transverse  position  x,  the 
angular  limits  of  neutrons  transmitted  though  the  curved  guide  for  1  > 
A.  are  described  by  the  parabolas  xp{xj  =MlN^r^  +ipl,d'-^^^^,)+ip^{x/W -111), 
where  5.  ,^  =  x^l.  From  this  equation  and  the  definition  of  2^,  it  can  be 
shown  that  when  5    <  9    ,  the  divergence  at  the  inner  radius  of  the 

c  in  c  our  O 

curved  guide  exit,  |  f{x  =  -  W2)  | ,  "saturates"  at  5    for  all  wavelengths 
exceeding  A'=  \{k^„„  I yjiK^,,  -^l)]  and  increases  parabolically  to 
|y/(.Y  =  +IV /2)\  =        + at  the  outer  radius  (see  dashed  curve  in  Fig. 
1).  Note  that  for  5  .  >  9     (i.e.,2''^oo)  or  for  2  <  2  <  2'  this  saturation 
is  never  achieved  because  \f{x  =  -W/2)\  is  always  less  than  5^    see  solid 
parabola  in  Fig.  1. 

FIGURE  1:The 
acceptance  diagram 
+W/2  (fansverse  position 
X  vs  trajectory  angle 
>  V       f]o^  B  curved/ 
straight  guide 
system  that  is 
constrained  by 

cm  cstr 

In  conclusion,  the  idealized  transmission  of  the  curved-straight 
combination  is  100  %  (that  of  the  unfiltered  straight  guide),  the 
spatial-angular  distribution  is  uniform,  and  the  divergence  is  hmited 
between  ±5    for  all  2  >  2^^^  when  the  following  four  conditions 
are  satisfied:  (i)     >  ^l{SWf),  (ii)  x^^^^>  x.^ ,  (iii)  W%^^2^,^,  and 

(iv)  2„,_,_>  {\lx.)^{2^lp)l^{{xjx.y-\).l}c^&  combination  guide 
acceptance  is  then  represented  by  the  solid  green  shaded  region  in  Rg.  1. 
The  unwanted,  divergent  part  rejected  by  the  straight  guide  is  represented 
by  the  cross-hatching. 

Example:  A  SANS  guide  requires  a  5  cm  wide  beam  with  a  divergence 
defined  by  the  critical  angle  of  natural  Ni  {x.^  =  x^^  =  1.73x10"^  radians 
A"^)  for  all  wavelengths  greater  than  2^^^,^  =  4  A.  Spatial  constraints 
impose  a  minimum  curved  guide  radius  of  1500  m.  Conditions  (i)  to  (iv) 
then  lead  to  L,  >  24.5  m,  xjx.^^  =  1.55  (i.e.,  5  ^^^^  is  55  %  greater  than  that 
of  natural  Ni),  and  L^^  >  7.23  m  (minimum  total  guide  length  is  31.7  m). 
This  is  technically  feasible. 
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A  Novel  Optical  Technique  for  Rapid  Detection  of  Neutrons 
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Almost  every  instrument  at  the  NCNR  and  other  neutron  scattering 
facihties  depends  on  ^He  proportional  tube  detectors  because  of 
their  high  efficiency,  good  background  rejection,  and  reliability.  New 
research  results  [1]  may  lead  to  a  new  detector  for  thermal  and  cold 
neutrons:  the  Lyman  Alpha  Neutron  Detector  (LAND).  This  detector, 
based  on  the  same  fundamental  nuclear  reaction  as  ^He  proportional 
tubes,  measures  ultraviolet  light  of  122  nm  wavelength  produced  by  the 
reaction  instead  of  amplifying  and  collecting  charge.  This  new  technique 
may  be  able  to  circumvent  limitations  of  ^He  proportional  tubes  while 
preserving  their  advantages  over  other  techniques. 

The  fundamental  nuclear  reaction  that  is  the  basis  of  the  LAND  is  n 
+  'He    t  +  p,  where  «  is  a  thermal  or  cold  neutron,  ^He  is  the  nucleus 
of  an  isotope  of  heUum,  and  /  and  p  are  a  tritium  nucleus  and  a  proton, 
respectively.  The  reaction  releases  763  keV  of  energy  that  is  shared 
between  t  and/>.  The  velocities  of  the  ^He  atomic  electrons  and  the 
outgoing  t  have  similar  magnitude,  which  makes  it  possible  to  have  a 
reaction  product  that  is  a  tritium  atom  (T)  instead  of  a  bare  tritium 
nucleus.  This  tritium  atom  is  hkely  to  be  formed  in  an  excited  state  that 
quickly  decays  to  its  ground  state,  emitting  one  or  more  photons  in  the 
process.  The  decay  from  the  first  excited  state  of  hydrogen  to  the  ground 
state  emits  a  characteristic  photon  at  121.6  nm,  the  "Lyman  Alpha 
Line"  (Ly  ).  Because  the  energy  needed  to  produce  Ly  light  is  10.2  eV, 
there  is  enough  energy  available  between  the  /  and  p  to  produce  tens  of 
thousands  of  photons. 

The  experimental  apparatus  (inset.  Fig.  1)  consisted  of  a  gas  cell  with  a 
neutron-transparent  window,  a  gas-handling  system  to  allow  evacuating 
and  filling  the  cell  with  ^He,  ■*He,  or  a  mixture  of  the  two,  and  a  high- 
efficiency  Ly  sensitive  photodetector  (a  Hamamatsu  R6835  "solar-bUnd" 
photomultipher,  which  detects  light  only  between  120  nm  and  180  nm). 
The  detector  was  installed  on  the  4.96  A  NG-6M  neutron  beam  hne. 

We  measured  photon  detection  rates  with  the  gas  cell  evacuated,  at 
various  pressures  of  "He,  and  at  various  pressures  of  ^He.  The  ''He  has 
essentially  no  interaction  with  neutrons,  and  provides  a  check  that  we 
were  not  seeing  signals  produced  by  beam-related  non-neutron  radiation. 
A  series  of  measurements  both  with  and  without  a  narrow  band-pass 
filter  demonstrated  that  the  signal  was  within  the  filter  bandwidth 
(8.7  nm  bandwidth  at  119.2  nm),  and  thus  almost  certainly  Ly^  hght. 

Rgure  1  shows  the  Ly  photon  yield  per  reacted  neutron.  If  the 
photon  production  were  a  result  of  the  n  +  ^He    t  +  p  reaction  alone. 


their  number  per  reacted  neutron  should  be  constant,  independent  of 
^He  pressure.  Instead,  the  data  fall  on  a  curve.  We  take  this  as  strong 
evidence  that  Ly^  photons  are  produced  in  atomic  interactions  with 
"spectator"  ^He  atoms  occurring  after  the  primary  nuclear  reaction. 
Preliminary  theoretical  calculations  suggest  that  most  of  the  radiation 
we  observe  comes  from  excitation  of  neutral  atoms  of  H  and  T  after  they 
have  been  slowed  to  below  1  keV.  At  93  kPa  (700  Terr),  46  photons  are 
produced  for  every  neutron  reacting  with^He.  This  high  yield  of  photons 
is  the  main  result  of  this  investigation. 
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FIGURE  1:  Lyman  alpha  photon  yield  per  reacted  neutron  as  a  function 
of  ^He  pressure.  The  Inset  states  the  nuclear  reaction  and  the  cartoon 
indicates  photon  production  and  the  experimental  setup. 

One  way  in  which  the  LAND  can  improve  upon  current  technology  is 
in  significantly  reducing  the  time  signal  per  detection  event.  Observed 
pulses  were  in  the  nanosecond  range,  compared  to  microseconds  for 
typical  ^He  proportional  tubes.  Moreover,  charge  collection  along  a 
high-voltage  anode  wire  requires  a  cylindrical  geometry  that  has  variable 
efficiency.  The  LAND  would  not  be  constrained  to  this  geometry. 

In  sum,  carefully  caHbrated  experiments  on  a  new  type  of  neutron 
detector,  LAND,  using  a  ^He  gas  detector  near  atmospheric  pressure, 
showed  that  tens  of  Ly^  photons  were  generated  per  neutron  absorbed. 
In  recognition  of  its  promise  for  a  transformational  approach  to  neutron 
detection,  this  result  has  garnered  a  2008  R&D  100  award  for  LAND  as 
one  of  the  most  significant  technologies  developed  during  the  previous 
year. 
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Visualizing  Nanoscale  Magnetism  in  3D  with  Polarized  SANS 
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Understanding  how  long-range  magnetic  correlations  between 
nanoparticles  form  and  evolve  as  a  function  of  temperature, 
shape  anisotropy,  and  near-neighbor  particle  proximity  takes  on 
heightened  importance  as  technology  increasingly  utilizes  nanoscale 
magnetism.  For  example,  novel  hyperthermia  cancer  treatment 
depends  critically  on  the  ability  of  magnetic  nanoparticles  to  couple 
strongly  for  localized,  high  heat  output.  At  the  other  extreme  as 
densities  of  =  TByte/cm^  (TByte/in-)  are  approached,  magnetic 
interparticle  cross-talk  in  hard  disk  nanopatterned  media  must  be 
minimized  to  avoid  unintentional  bit  reversal  during  reading  and 
writing.  Both  examples  motivate  the  development  of  techniques 
for  three  dimensional  (3D)  profiling  of  these  magnetic  nanoparticle 
assemblies  and  this  is  where  polarized  small  angle  neutron  scattering 
(SANS)  excels. 
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FIGURE  1:  Magnetic  and  nuclear  separation  in  a  remanent  magnetic 
field.  The  left  inset  is  an  asymmetric  spin-flip  scattering  SANS  pattern. 
The  right  inset  shows  nuclear  scattering  at  larger  Q.  In  the  main  plot 
nuclear  scattering  (black)  Is  distinct  from  magnetic  scattering  (red 
and  green:  perpendicular  to  neutron  polarization,  blue:  parallel). 

Selecting  the  neutron  polarization  before  small  angle  scattering  from 
a  sample  and  analyzing  it  after  allows  us  to  distinguish  between 
non-spin  flip  scattering  (involving  nuclei  and  magnetic  moments 
parallel  to  the  neutron  spin)  and  spin-flip  scattering  (magnetic 
moments  perpendicular  to  the  neutron  spin).  We  can  use  these 
measurements  to  extract  the  magnetism  and  separate  it  into  its 
average  3D  components.  This  is  important  in  systems  where  the 
nuclear  scattering  dominates  the  magnetic  signal  or  when  magnetic 
saturation  cannot  be  achieved  and  used  as  a  background  [1]. 

The  key  to  performing  polarization  analysis  on  a  divergent  scattered 
beam  is  to  use  a  ^He  spin  polarized  filter.  These  are  excellent  neutron 
spin  polarizers/analyzers  because  of  their  high  absorption  asymmetry 
for  parallel  and  anti-parallel  neutron  spins,  and  they  can  be  reversed 


with  negligible  loss  [2].  Their  polarizing  efficiency,  however,  decays  with  a 
lifetime  of  150  h  to  200  h  in  our  SANS  set-up.  Based  on  the  formalism  of 
Re£  3,  we  have  implemented  an  algorithm  to  extract  all  four  polarization- 
corrected  spin  scattering  cross-sections  by  explicitly  including  the  ^He  time 
dependence. 

Our  samples  consist  of  ferromagnetic,  monodispersed,  magnetite 
nanoparticles  7  nm  in  diameter  that  exhibit  strong  magnetic  coupling  with 
an  average  edge-to-edge  separation  of  2.5  nrn  [4].  The  nanoparticles  form 
a  powder  of  hexagonal  close  packed  thin  (several  monolayer)  crystallites  on 
the  order  of  a  micrometer  in  the  plane.  A  range  of  magnetic  behaviors  is 
experimentally  accessible  because  the  maximum  temperature  at  which  bulk 
ferromagnetism  is  observed  (the  blocking  temperature)  is  well  below  room 
temperature. 

To  test  our  algorithm,  we  collected  data  at  50  K  in  a  remanent  magnetic 
field  (2.5  mT  to  3.0  mT)  where  the  magnetic  scattering  should  be  nearly 
homogenous  along  the  directions  perpendicular  (Y  and  Z)  and  parallel  (X) 
to  the  neutron  polarization.  As  shown  in  Fig.  1,  the  polarization-corrected 
data  indeed  exhibits  this  pattern.  Even  though  the  magnetic  scattering  is 
much  smaller  than  the  nuclear  scattering,  it  can  be  fully  resolved  into  its 
vectorial  components  (which  happen  in  this  case  to  be  nearly  the  same). 

Utilizing  this  technique,  we  studied  the  dependence  of  the  magnetic 
correlation  length  on  temperature  from  50  K  to  300  K.  While  the 
average  structural  separation  remained  constant,  the  long-range  magnetic 
correlations  showed  marked  decrease  well  past  their  ferromagnetic 
blocking  temperature  at  65  K.  Quantitative  analysis  similar  to  Ref  5 
reveals  that  the  magnetic  domains  ranged  from  100  nm  at  50  K 
{-  10  nanoparticles)  to  10  nm  («  1  nanoparticle)  at  300  K. 

In  conclusion  we  have  shown  that  'He  spin  filters  with  the  appropriate 
time-dependent  polarization  correction  can  be  used  to  perform  full 
polarization  analysis  and  subsequent  3D  profiling  of  magnetic  systems. 
We  envision  a  suite  of  neutron  magnetic  techniques  usefiJ  for  investigating 
a  variety  of  3D  systems  such  as  biotagged  magnetic  sensors,  in  situ 
magnetically-manipulated  bioprobes,  controlled  viscosity  magnetic  gels, 
self-assembled  magnetic  templates,  and  patterned  thin-film  magnetic 
media. 
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Neutron  Phase  Imaging  with  Talbot  Optics 

D.S.  Hussey\  S-W.  Lee^  D.L.  Jacobson\  and  M.  Arif 


eutxon  phase  imaging  can  be  potentially  10  to  100  times  more 
sensitive  than  traditional  neutron  radiography,  see  Fig.  la,b. 


However,  it  has  seen  hmited  appHcation  because  the  quasi-coherence 
required  of  the  propagation  technique  reduces  the  incident  neutron 
intensity  by  about  IC,  which  prevents  realizing  the  improved  sensitivity. 
By  making  use  of  the  Talbot  effect  it  is  possible  to  form  phase  images 
with  only  a  factor  of  two  loss  in  intensity,  which  makes  neutron  phase 
imaging  feasible. 


(c)    source  grating  (GO)       phase  object   phase  grating  analyzer  grating 


FIGURE  1:  Comparison  of  the  transmission  image  (a)  to  the  neutron  phase 
gradient  image  (b]  of  a  U.S.  quarter.  The  features  of  the  quarter  are  not 
resolved  in  the  transmission  due  to  insufficient  contrast.  However,  the 
eagle  is  clearly  visible  in  the  neutron  phase  gradient  image.  The  setup  is 
shown  schematically  in  (c). 


The  Talbot  effect  is  that  a  coherendy  illuminated  phase  grating  produces 
a  self-image  at  a  specific  downstream  distance  called  the  Talbot  distance. 
Placing  an  object  to  be  imaged  in  front  of  this  phase  grating  (Fig.  Ic) 
results  in  a  phase  shift  of  the  self-image  that  is  Unearly  related  to  the 
neutron  phase  gradient  produced  in  the  object.  The  neutron  beam 
must  be  coherent  only  along  the  direction  parallel  to  the  phase  grating 
but  perpendicular  to  the  grating  lines.  In  our  case  the  beam  is  formed 
with  a  source  grating  (GO)  («  10  |im  thick  Gd  Hnes  with  a  ~  0.8  mm 
period,  60  %  duty  cycle)  which  reduces  the  incident  neutron  intensity 
by  60  %.  The  phase  grating  is  composed  of  lines  of  Si  that  are  30  |j.m 
thick  and  4  pm  wide  periodically  spaced  by  8  pm.  The  period  of  the 
self-image  is  smaller  than  current  neutron  detector  resolution.  However, 
the  self-image  phase-shift  can  be  measured  by  placing  an  absorbing 
analj^er  grating  at  the  Talbot  distance,  and  then  scanning  the  source 
grating  through  one  period.  The  analyzer  grating  is  composed  of  lines 
of  Gd  that  are  3  pm  thick  and  2  pm  wide  periodically  spaced  by  4  pm. 


The  gratings  were  fabricated  using  the  NIST  Nanofabrication  facility, 
and  installed  in  a  prototype  phase  imaging  facility  in  the  NG7  Neutron 
Interferometer  and  Optics  Facility. 

To  obtain  a  neutron  phase  image  with  Talbot  optics,  a  set  of  neutron 
images  is  acquired  as  the  source  grating  is  stepped  through  one  period, 
analogous  to  a  traditional  neutron  interferogram.  Making  a  pixel- 
by-pixel  comparison  of  the  mean,  amphtude,  and  phase  shift  of  the 
interference  patterns  of  the  object  to  the  open  beam,  one  obtains  the 
transmission  image,  the  dark-field  image,  and  the  phase  gradient 
image,  respectively,  as  shown  in  Fig.  2.  The  neutron  phase  image  can  be 
obtained  by  integrating  the  phase  gradient  image  along  the  dimension 
perpendicular  to  the  phase  grating  lines. 


Distance  (mm) 


FIGURE  2:  Talbot  optics  measurement,  (a]  A  sapphire  wedge  and 
single  crystal  aluminum  cylinder  are  a  simple  model  system,  (b)  The  raw 
Image  with  the  source  grating  (GO)  stepped  to  the  4th  position,  (c)  The 
interference  pattern  of  points  A,  B,  and  C  labeled  in  (b).  (d)  The  neutron 
phase  gradient  image,  and  a  plot  of  the  average  gradient  and  the  expected 
phase  gradient  for  the  wedge  and  cylinder,  (e)  The  transmission  radiograph, 
(f)  The  dark-field  image,  which  clearly  shows  the  rough  edges  of  the 
aluminum  cylinder. 

With  the  potentially  greater  sensitivity  than  radiography  it  may  be 
possible  to  study  residual  stress  in  materials  through  the  change  in 
number  density.  Since  the  dark-field  image  is  affected  by  the  porosity 
of  the  material,  Talbot  optics  would  facilitate  tomography  of  porous 
materials  such  as  metal  foams. 


'Physics  Laboratory,  National  Institute  of  Standards  and  Technology,  Gaithersburg,MD  20899-8461 
^Korea  Atomic  Energy  Research  Institute,  Daejeon,  South  Korea 


NCNR  2008  Accomplishments  and  Opportunities 


Neutron  Imaging  Straw  Detectors:  Getting  the  Efficiency  Right 

R.G.  Downing^  L.R.  Cao\  J.L.  Lacy^ 


Straw  detectors  are  a  new  concept  in  position  sensitive  neutron 
detectors  recently  developed  by  Proportional  Technologies, 
Inc.  [l].The  detectors  are  especially  valuable  in  large-area  imaging 
applications  and  scattering  instruments  [2]  and  are  being  installed 
at  the  Spallation  Neutron  Source  (SNS)  at  Oak  Ridge  National 
Laboratory.  The  atomic  thickness  of  boron  carbide  deposits  on  the 
individual  tubes  is  critical  to  the  efficiency  of  the  overall  detector. 
Neutron  Depth  Profiling  (NDP)  [3],  a  nondestructive  metrological 
technique,  has  played  a  critical  role  in  optimizing  the  performance 
of  these  novel  high  resolution  neutron  systems. 

Neutrons  are  not  easily  detected  with  high  resolution  and  good 
efficiency.  Consequently  high-quality  neutron  detectors  remain  a 
challenge  to  manufacture.  A  few  nuclides  such  as  ^He,  ^"B  and  ''Li 
are  commonly  used  to  convert  impinging  neutrons  into  an  electrical 
pulse  for  electronic  post-processing.  The  two  most  commonly 
used  neutron  position  sensitive  detectors  are  the  ^He  proportional 
counter  and  the  ZnSALiF  scintillator.  The  first  has,  perhaps,  the 
highest  detection  efficiency,  but  is  costly  to  make  and  the  spatial 
resolution  is  relatively  poor.  The  latter  has  better  spatial  resolution, 
but  suffers  as  an  optimal  solution  since  the  afterglow  in  a  scintillator 
can  last  for  up  to  10  |is,  making  it  unsuitable  for  any  high  count 
event  application  [4]. 


The  neutron-active  surface  of  a  straw  detector  is  sputter-coated  with 
enriched  boron  carbide  ("^B^C).  In  operation,  two  highly  energetic 
charged  particles  are  generated  by  neutron  capture  within  the  deposited 
width.  One  of  those  two  particles  must  escape  the  deposit  and  ionize  gas 
within  the  hollow  straw  to  generate  a  signal  pulse.  The  signal  is  read  out 
via  the  suspended  wire  running  through  the  center  of  the  straw.  Many 
straws  are  close-packed  into  a  2D  array  to  form  a  large  area  neutron 
position-sensitive  detector  (Fig.  la). 

The  ^'-^B  atomic  density  of  '''B_|C  coating  applied  to  the  straw  strongly 
influences  detector  efficiency  and  overall  performance.  The  ^°B^C 
coating  must  be  sufficiendy  thin  to  permit  efficient  escape  of  one  or  the 
other  reaction  product  emitted  upon  each  neutron  reaction  with  a  '"B 
atom.  Otherwise,  the  event  is  lost  in  self  absorption.  Furthermore,  if 
the  film  is  too  thick,  remaining  neutrons  cannot  continue  along  their 
trajectories  to  the  next  straw,  leading  to  further  inefficiency.  Obviously, 
the  accurate  determination  of  '"B  depth  distribution  (stoichiometry)  and 
total  mass  -  not  linear  thickness  -  of  this  thin  (0.1  |im)  lining  is  critical 
in  the  manufacture  of  straw  elements.  The  required  metrology  is  readily 
achieved  using  NDP. 

Sample  coupons  were  selected  representing  material  along  the  length 
and  width  of  the  "'B_|C  deposition  chamber.  NDP  spectra  were  then 
determined  from  each  sample;  three  spectra  are  overlaid  and  presented 
in  Fig.  lb  illustrating  the  variations  in  the  atomic  thickness  of  B  found. 
NDP  results  revealed  that  the  masks  used  in  applying  the  film  yielded 
a  smoothly  varying  thickness  in  coating  across  the  chamber,  see  Fig.  Ic. 
Using  the  metrology  provided  by  NDP,  a  practical  and  reliable  approach 
to  compensate  the  deposition  mask  was  provided  as  determined  in 
subsequent  production  runs. 
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FIGURE  1:  (a]  Illustration  of  a  3D  stack  of  straw  detectors.  The  balloon 
shows  the  neutron-''°B  reaction,  (b)  An  overlay  of  three  NDP  spectra 
for  samples  deposited  with  different  masses  of  ^°B^C.  Only  the 
1 .472  keV  alpha  particle  portion  of  the  spectrum  is  shown,  (c)  A  plot  of 
the  boron  carbide  thickness  versus  position  in  the  deposition  chamber 
as  determined  by  NDP  in  the  first  production  lots. 
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Neutron  Source  Operations 


The  NIST  reactor  (NBSR)  operated  for  253  full  power  (20  MW)  days  or  99  %  of  the  scheduled  time  during  FY2008.  A  typical  operating 
year  consists  of  seven  cycles.  A  cycle  has  38  d  of  continuous  full  power  operation  followed  by  11  d  of  shutdown  for  maintenance,  refueling, 
and  startup  preparations.  The  cold  source  continues  to  operate  at  a  remarkable  99.99  %  reliability.  This  outstanding  record  of  safe  and  reliable 
performance  of  the  reactor  and  cold  source  lies  at  the  foundation  of  the  success  of  the  neutron  scattering,  nuclear  methods,  and  neutron  physics 
programs.  Licensed  operations  continue  while  the  Nuclear  Regulatory  Commission  nears  completion  of  the  renewal  of  the  reactor  license, 
which  will  permit  another  20  years  of  operation.  Some  of  the  activities  of  the  Reactor  Operations  and  Engineering  Group  are  outlined  below. 


Facility  Improvements 

Shim  Arm  Radiography:  The  NIST  Reactor  is  controlled  by 
shim  arms,  which  consist  of  neutron  absorbing  cadmium,  lodged 
between  two  supporting  aluminum  tubular  structures.  The  cadmium 
gets  depleted  ("burned")  in  the  process  of  controlling  the  reactor. 
It  is  important  to  monitor  the  rate  of  Cd  burnup.  A  radiographic 
method  (Fig.  1)  has  been  developed  and  used  to  characterize  the  Cd 
distribution  of  four  new  shim  arms.  Radiography  wiU  be  used  again 
at  the  ends  of  their  life  cycles  to  determine  the  cadmium  depletion 
profile. 


Shim  arm  t>eing  inserted  into 
if'  the  neutron  radiography  box. 


FIGURE  1:  Shim  arm  Cadmium  burnup  monitored  by  neutron 
radiography.  As  an  example  of  resolution,  the  inset  radiogram  shows 
an  edge  of  the  Cd  layer. 


Thermal  Shield  Flow  Tests:  The  apparatus  for  detecting  leaks  in  pipes 
surrounded  by  a  vacuum  described  in  the  2007  report  has  been  tested  in 
a  number  of  critical  situations  involving  water  flow  in  copper  tubing  in 
the  presence  of  corrosive  media,  both  in  a  laboratory  setting  and  in  situ. 
These  tests  continue  to  provide  useful  information  to  improve  the  safety 
of  operations.  A  contract  to  study  the  feasibility  of  a  fuU  scale  system  has 
been  negotiated. 

Electric  Panel  Upgrade:  Reactor  Operations  is  collaborating  with  the 
NIST  Plant  Division  to  replace  original  equipment  motor  control  centers 
(MCC)  with  state  of  the  art  units.  The  new  units  allow  the  NCNR 
increased  reliability  and  several  special  features  including  relaying  the 
status  of  the  MCC  directly  to  the  reactor  control  room. 

Thermal  Shield  Coohng  Water  Filter:  During  testing  it  was  discovered 
that  the  dose  rate  emanating  from  the  thermal  shield  headers  is  originating 
from  the  suspension  of  solids  in  the  cooling  water.  A  design  for  a 
continuously  operating  suspended  particle  filtering  system  has  been  made 
and  is  ready  to  be  installed. 

Control  Room  Upgrade:  The  modernization  of  the  Reactor  Control 
Room  initiated  with  Brookhaven  National  Laboratory  last  year  has 
completed  the  analysis  stage  and  bids  for  proceeding  with  various  aspects 
of  the  upgrade  have  been  received  and  are  being  evaluated. 


New  Staff 

This  year  the  Reactor  Operations  and  Engineering  Group  hired  several  new  staff  members.  Two  senior  reactor  operator  candidates, 
Sam  MacDavid  and  Sam  Colvard,  are  on  pace  to  be  licensed  in  October  of  2008.  Electronic  Technician  Susan  Deeb  wiU  participate  in 
instrumentation  development  and  calibrations  as  well  as  general  electronics  trouble-shooting.  Joe  Reyenga  was  hired  as  an  Electronics  Engineer. 
His  main  responsibility  will  be  to  lead  the  modernization  project  of  the  NBSR  Control  Room.  Mechanical  Engineer  Mike  Middleton  will 
focus  on  thermo-mechanical  aspects  of  the  Cold  Source  and  its  operational  integration  with  the  reactor  control  systems.  We  are  very  pleased  to 
welcome  these  talented  new  staff  members. 


FIGURE  3:  New  Reactor  Operations 
and  Engineering  staff  members:  Senior 
operator  candidates;  Sam  MacDauid  and 
Sam  Colvard,  Electronics  Technician  Susan 
Deeb,  Electronics  Engineer  Joe  Reyenga, 
and  Mechanical  Engineer  Mike  Middleton. 
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Facility  Developments 


Developments  at  the  NCNR  facility  to  improve  the  quality  of  the 
user  experience  are  an  ongoing  commitment.  Recent  progress  in 
instruments,  software,  sample  environment  instrumentation,  and  user 
laboratories  is  described  below. 

MACS 

The  Multi-Axis  Crystal  Spectrometer  (MACS)  is  a  high  intensity  cold 
neutron  spectrometer  that  will  gready  expand  access  to  atomic  scale 
dynamic  processes  in  condensed  matter. The  fabrication  and  installation 
of  this  instrument  has  progressed  to  the  point  where  all  major 
components  have  now  been  assembled  on  the  experiment  floor. 

During  the  past  year,  work  has  focused  on  the  MACS  detector  system 
which  consists  of  20  channels,  each  of  which  subtends  a  solid  angle  of 
0.1  sr  at  the  sample.  The  central  functional  element  in  each  channel 
of  the  detector  system  is  the  double  crystal  analyzer  (DXAL).  A 
single  stepping  motor  positions  the  device  for  scattered  beam  energy 
analysis  in  the  range  from  2.5  meV  to  15  meV.  Transmission  through 
the  DXAL  is  improved  through  the  use  of  84'  graphite,  40  %  coarser 
than  the  60'  graphite  on  the  monochromator.  This  choice  ensures 
similar  elastic  energy  resolution  before  and  after  the  sample  despite 
the  different  geometries. The  nearly  complete  analyzer  assembly 
now  carries  4500  kg  of  polyethylene  and  boron  carbide  loaded 
polyurethane  to  create  a  quiet  environment  for  the  forty  'He  neutron 
detectors.  To  simplify  the  interfice  to  this  mobile  detector  system, 
compartments  below  the  analyzer  assemblies  accommodate  the  wide 
range  of  electronic  support  systems  needed  for  operation.  The  sample 
positioning  system  which  provides  fill  gimballing  capabilities  for  a  400 
kg  sample  environment  system  has  also  been  installed  this  year.  Testing 
of  the  complete  system  is  now  underway.  The  scientific  productivity 
of  both  of  the  NCNR's  previous  cold  neutron  spectrometers  SPINS 
and  DCS,  has  been  remarltable.  And  P^ACS  promises  to  usher  in  an 
exciting  new  era  for  cold  neutron  spectroscopy.  In  particular,  the  new 
instrument  will  allow  scientists  to  tackle  the  many  important  problems 
in  modern  condensed  matter  physics  in  which  only  sm-all  samples  are 
available  or  in  which  extreme  sample  environment  conditions  that  limit 
sample  size  are  required.  Thus  we  eagerly  anticipate  the  first  scattered 
neutrons  to  be  detected  by  MACS  in  late  2008. 


Magnetic  i-ield  Abatenieni 

Stray  fields  from  superconducting  magnets  used  as  sample 
environments  interfere  with  the  operation  of  neutron  measurement 
methods  that  depend  on  precession  of  the  neutron  spin.  The  NCNR 


has  controlled  stray  fields  administratively  by  restricting  the  use  of  magnets 
during  some  periods,  thereby  allowing  more  "sensitive"  experiments  to  be 
performed  on  the  NCNR's  neutron  spin-echo  spectrometer  (NSE).  This 
situation  not  only  limits  the  number  of  experiments  that  can  employ  high 
field  magnets,  it  also  reduces  the  eff"ectiveness  of  the  spin  echo  instrument 
during  those  periods  that  magnets  are  allowed  to  be  used.  Thus  techniques 
that  shield  or  reduce  stray  fields  are  necessary  in  order  to  achieve  the 
maximum  scientific  output  of  the  facility. 

Recently  the  NCNR  completed  a  magnetic  field  cancellation  demonstration 
project  for  SPINS  which  is  located  near  the  NSE  instrument. This  system, 
fuUy  automated  through  the  SPINS  instrument  control  software,  sends 
compensating  currents  through  two  rings  encircling  the  SPINS  instrument 
as  a  function  of  the  superconducting  magnet  field  applied  to  the  sample  and 
the  orientation  of  the  spectrometer.  As  a  result,  NSE  is  able  to  continue 
collecting  high-quality  data  even  when  superconducting  magnets  are  in  use 
on  the  adjacent  SPINS  instrument.  This  proof-of-concept  demonstration 
has  provided  valuable  practical  information  on  the  feasibility  of  magnetic 
field  cancellation  for  potential  deployment  on  other  instruments.  The 
NCNR  now  intends  to  develop  and  deploy  cancellation  coils  on  those 
instruments  which  use  high  magnetic  fields  in  the  NCNR's  guide  hall.  This, 
along  with  moving  the  NSE  instrument  to  a  more  magnetically  isolated 
location  at  the  new  guide  NGa  as  part  of  the  expansion  wall  reduce  the 
variability  of  stray  magnetic  fields  to  less  than  the  natural  variation  of  the 
Earth's  field. 

Standardization  of  Data  Acquisition  Electronics 

This  year  the  NCNR  completed  a  program  of  instrument  upgrades  intended 
to  standardize  the  data  acquisition  hardware  of  all  neutron  scattering 
instruments.  The  common  suite  of  motor  and  counter  controls  allows  more 
straightforward  expansion,  improves  the  reliability,  and  reduces  downtime 
due  to  greater  ease  in  troubleshooting.  Instrument  controls  on  the  NSE 
instrument  were  migrated  to  NCNR  standards  early  this  year  resulting  in 
increased  reliability. 

^He  Neutron  Spin  Filters 

The  development  and  application  of  ^He  neutron  spin  filters  (NSF)  has 
continued  this  past  year.  'He  NSFs  were  employed  in  20  experiments 
with  more  than  100  bar-liters  of  polarized  ^He  delivered  over  the  past  year 
alone.  These  devices  have  been  used  on  SANS,  reflectometry  and  triple-axis 
instruments.  Methods  and  software  have  been  developed  that  allow  users 
to  account  for  the  polarization  decay  of  the  cells  and  to  account  for  differing 
paths  through  the  ^He  cell.  (See  page  50)  We  are  actively  working  to 
develop  and  employ  'He  cells  on  even  more  instruments,  particularly  those 
with  detectors  covering  large  solid  angles. 

In  the  past,  precession  coil  flippers  have  been  used  to  control  the  neutron 
spin  of  the  scattered  beams.  However,  we  have  found  that  the  aluminum 
wires  in  such  flippers  produces  weak  small  angle  scattering,  hence  such  a 
flipper  would  be  unsuitable  after  the  sample  for  polarized  beam  experiments 
in  SANS  and  diffiise  reflectometry.  We  have  developed  a  method  that 
utilizes  the  adiabatic  fast  passage  nuclear  magnetic  resonance  technique 
to  invert  the  'He  polarization,  hence  flipping  neutron  polarization  in 
a  compact,  tight  space  on  the  beam  Une.  This  method  integrates  the 
spin  flipper  and  the  analyzer  into  a  single  neutron-polarizing  device. 


Furthermore,  flipping  the  neutron  spin  with  this  method  is  independent  of 
wavelength  and  very  efficient. 

One  of  the  primary  challenges  in  the  deployment  of  ^He  NSFs  is  the 
reduced  lifetime  of  the  polarization  in  the  field  due  to  inhomogeneities  in 
the  magnetic  field  environment.  To  extend  the  lifetime  to  that  observed 
in  the  lab,  we  are  pursuing  new  methods  to  optimize  the  magnetic  field 
environment  of  the  Ln-situ  ^He  cells.  Previous  magnetic  field  coUs  were 
large  in  order  to  provide  for  a  uniform  field  around  the  ^He  polarized  gas 
cell,  while  the  room  available  for  the  cells  on  most  instruments  is  limited. 
Detailed  calculations  of  the  fields  have  been  performed  for  a  smaller  device 
that  has  been  fabricated  and  employed  on  the  BT-7  triple  axis  instrument. 
This  development  not  only  is  an  important  enhancement  of  our  current 
capabilities,  but  is  also  a  key  development  for  accommodating  'He  cells  on  a 
wider  range  of  NCNR  instruments. 

Sample  Environment 

The  vast  majority  of  neutron  scattering  experiments  take  place  with  the 
sample  in  a  special  environment.  User  surveys  routinely  emphasize  the 
importance  of  high  quality,  reliable  sample  environment  equipment  for  the 
success  of  their  experiments.  To  this  end,  the  NCNR  currendy  maintains 
and  operates  over  120  different  pieces  of  sample  environment  equipment. 
In  addition,  the  NCNR  continues  to  expand  our  capabilities  through  the 
acquisition  of  new  equipment  for  NCNR  users.  During  the  past  year, 
these  include  a  15  T  superconducting  magnet  with  a  dilution  refrigerator 
insert  and  a  1.36  MPa  high-pressure  helium  gas  pressurization  rig  along 
vrith  a  variety  of  pressure  cells  including  a  high-pressure  cryogenic  cell  for 
the  SANS  instruments.  In  addition,  we  continue  to  improve  the  reliability 
and  performance  of  our  current  suite  of  equipment.  For  instance,  the  9  T 
horiTOntal  field  magnet  for  SANS  has  been  made  more  robust  through 
the  use  of  titanium  alloy  windows.  Fmally,  we  are  working  to  improve  our 
tracking  of  the  reliability  of  our  sample  environments  through  the  NCNR 
Information  Management  System.  This  should  allow  us  to  address  issues 
before  they  become  major  problems. 

Software  Development 

The  NCNR  released  a  new  version  of  our  popular  SANS/USANS  data 
reduction  and  analysis  suite  which  includes  many  new  features  addressing 
frequent  user  requests.  A  new,  more  intuitive  interface  for  analysis  has 
been  developed  that  allows  simpler  control  and  manipulation  of  data  and 
data  modeling.  The  performance  of  the  data  fitting  routines  has  been 
improved  by  as  much  as  a  factor  of  10  by  using  optimized  C  routines 
and  exploiting  multi-threading  capabilities  of  modern  computing 
platforms.  This  significant  performance  enhancement  has  had  a  direct 
impact  on  computation-intensive  analysis  operations  making  it  possible  to 
simultaneously  fit  SANS  and  uSANS  data  sets  and  even  parametric  studies 
with  multiple  temperatures,  concentrations,  etc.  Video  tutorials  covering 
experiment  planning,  data  reduction  and  analysis  now  supplement  the 
software  documentation  and  can  be  found  on  the  NCNR  website:  http:// 
www.ncnr.nist.gov/programs/sans/data/red  anal.html 

New  features  and  enhancements  continue  to  be  incorporated  into  the  DAVE 
software  suite  for  inelastic  neutron  scattering  data.  A  significant  redesign 
of  DAVE,  known  as  DAVE  2,  was  released  this  year  DAVE  2  includes 
substantial  improvements  such  as  easier  session  management,  simpler 
handling  of  multiple  datasets,  and  enhanced  usability.  Users  benefit 
from  the  new  functionality  in  a  much  more  integrated  application.  Two 
application  development  courses  were  offered  this  year  in  order  to  train 
NCNR  technical  staff  interested  in  developing  applications.  Twenty 
participants  from  NIST  and  elsewhere  attended  both  courses.  Such 


training  is  crucial  in  order  to  increase  the  number  of  scientists  that 
contribute  to  the  development  of  DAVE.  For  example,  developers  at  the 
Paul  Scherrer  Institut  in  Switzerland  now  contribute  data  reduction  support 
for  two  of  their  instruments.  In  the  past  year  use  of  DAVE  was  cited  in  more 
than  30  publications. 

NCNR  User  Accounts  for  Proposals  and  Facility  Access 
Substantial  progress  has  been  made  during  the  past  year  with  the 
development  of  a  number  of  different  software  applications  that  are  intended 
to  expedite  and  enhance  the  user  experience  at  the  NCNR.  Among  the 
most  visible  of  these  is  the  NCNR-IMS  (Information  Management 
System),  which  is  a  web-based,  database  application  that  enables  NCNR 
users  to  submit  beamtime  requests/proposals,  copy/edit  old  requests/ 
proposals  to  expedite  new  submissions  (this  avoids  having  to  retype 
collaborator  names  and  affiliations),  submit  visitor  registration  requests, 
request  NCNR  health  physics  training,  or  edit  their  account  information.  At 
the  same  time  the  NCNR-IMS  tracks  various  user  statistics  and  important 
facility  metrics  such  as  the  number  of  NCNR  research  participants, 
publications,  and  instrument  usage,  many  of  which  are  included  in  this 
Annual  Report  and  are  essential  to  effective  NCNR  management.  Presentiy 
all  NCNR  beam  stations  are  scheduled  through  the  NCNR-IMS,  save  for 
the  two  SANS  instruments,  and  all  schedules  are  published  on  the  web  for 
users  to  see.  The  NCNR-IMS  also  automatically  informs  users  when  their 
health  physics  training  is  due  to  expire,  automatically  requests  dosimetry  for 
users  when  they  register  to  visit  the  NCNR  for  the  purpose  of  running  an 
experiment,  and  allows  instrument  scientists  to  send  out  automated  email 
messages  to  inform  users  when  their  beamtime  has  been  scheduled.  The 
NCNR-IMS  also  offers  an  email  "support"  address  to  which  users  can  report 
problems,  ask  for  assistance,  or  request  specific  new  functionality. 
Ongoing  enhancements  include  the  development  of  an  improved  interface 
to  allow  users  to  submit  (and  track)  their  publications,  invited  talks,  and 
awards  that  are  based  on  data  obtained  at  the  NCNR  to  the  User  Office 
through  the  NCNR-IMS  application  for  timely  processing.  In  addition,  the 
NCNR  is  in  the  process  of  incorporating  the  health  physics  training  into  the 
NCNR-IMS.  This  will  allow  users  to  take  the  computer-based  training  at 
their  leisure  using  their  own  computers  at  their  home  institutions  (users  will 
still  have  to  take  the  mandatory  tour  upon  arrival  at  the  NCNR  however). 

User  Laboratories 

The  NCNR  strives  to  provide  safe  and  functional  laboratories  for  the  use  of 
our  visitors  and  staff.  In  2007  we  renovated  a  laboratory  to  support  sample 
handing  for  the  wide  variety  of  scientific  programs  in  hard  condensed 
matter  science.  The  laboratory  is  now  equipped  with  a  specialized  glove 
box  for  the  loading  and  unloading  of  sample  canisters,  a  high  temperature 
tube  fijrnace,  and  other  basic  laboratory  equipment.  One  of  the  NCNR 
laboratories  was  equipped  with  a  specialized  fume  hood  for  working  with 
acids.  The  fume  hood  contains  three  different  tanks  for  acid  work  and  is 
connected  to  the  facility's  acid  neutralization  tank.  In  addition  to  laboratory 
renovations,  new  laboratory  equipment  has  been  purchased  to  improve 
the  research  capabilities  of  the  NCNR.  Some  examples  include  three  new 
high  temperature  furnaces  for  solid  state  synthesis  of  materials,  a  dry  bath,  a 
microscope,  and  a  vacuum  oven. 

In  addition  to  improving  the  laboratory  facilities,  we  have  continued  to 
improve  training  in  the  safe  use  of  NCNR  laboratories  and  equipment. 
To  enhance  and  standardize  training  for  NCNR  staff,  a  PowerPoint 
presentation  and  test  was  developed.  This  format  ensures  that  NCNR  staff 
is  well-versed  in  their  roles  and  responsibilities  for  safety  in  the  NCNR 
laboratories.  In  addition,  we  have  developed  a  new  user  manual  and 
additional  training  on  the  proper  handling  and  disposal  of  acids. 
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'he  mission  of  the  NIST  Center  for  Neutron  Research  is  to 
_L  assure  the  availabihty  of  neutron  measurement  capabihties 
to  meet  the  needs  of  U.S.  researchers  from  industry,  academia  and 
from  other  U.S.  government  agencies.  To  carry  out  this  mission, 
the  NCNR  uses  several  different  mechanisms  to  work  with 
participants  from  outside  NIST,  including  a  competitive  proposal 
process,  instrument  partnerships,  and  collaborative  research  with 
NIST. 

PROPOSAL  SYSTEM 

Most  of  the  time  on  NCNR  instruments  is  made  available 
through  a  competitive,  peer-review  proposal  process.  The 
NCNR  issues  calls  for  proposals  approximately  twice  a  year. 
Proposals  are  reviewed  at  several  different  levels.  First,  expert 
external  referees  evaluate  each  proposal  on  merit  and  provide 
us  with  written  comments  and  ratings.  This  is  a  very  thorough 
process  where  several  different  referees  review  each  proposal. 
Second,  the  proposals  are  evaluated  on  technical  feasibility 
and  safety  by  NCNR  staff.  Third,  we  convene  our  Beam  Time 
Allocation  Committee  (BTAC)  to  assess  the  reviews  and  to 
allocate  the  available  instrument  time.  Using  the  results  of  the 
external  peer  review  and  their  own  judgment,  the  BTAC  makes 
recommendations  to  the  NCNR  Director  on  the  amount  of 
beam  time  to  allocate  to  each  approved  experiment.  Approved 
experiments  are  scheduled  by  NCNR  staff  members  in 
consultation  with  the  experimenters. 

The  current  BTAC  members  are: 

Ramanan  Krishnamoorti  (chair.  University  of  Houston) 
Andrew  Allen  (NIST  Ceramics  Division) 
Collin  Broholm  (Johns  Hopkins  University) 
Ken  Chen  (Sandia  National  Laboratories) 
Michael  Kent  (Sandia  National  Laboratories) 
Young  Lee  (Massachusetts  Institute  of  Technology) 
Raul  Lobo  (University  of  Delaware) 
Robert  Leheny  (Johns  Hopkins  University) 
Lee  Magid  (University  of  Tennessee,  Knoxville) 
Janna  Maranas  (The  Pennsylvania  State  University) 
Stephan  Rosenkranz  (Argonne  National  Laboratory) 
Dawn  Sumner  (University  of  California  Davis) 
Michael  Toney  (Stanford  Synchrotron  Radiation  Laboratory) 
Lynn  Walker  (Carnegie-Mellon  University) 
David  Worcester  (University  of  Missouri) 

PARTNERSHIPS 

The  NCNR  may  form  partnerships  with  other  institutions  to 
fund  the  development  and  operation  of  selected  instruments. 
These  partnerships,  or  "Participating  Research  Teams",  may 
have  access  to  as  much  as  75  %  of  the  available  beam  time  on 
the  instrument  depending  on  the  share  of  total  costs  borne  by 
the  team.  A  minimum  of  25  %  of  the  available  beam  time  is 
always  made  available  through  the  NCNR  proposal  program  to 
aU  users.  Partnerships  are  negotiated  for  a  fixed  period  (usually 
three  years)  and  may  be  renewed  if  there  is  mutual  interest  and 
a  continued  need.  These  partnerships  have  proven  to  be  an 


important  and  effective  way  to  expand  the  research  community's  access 
to  NCNR  capabilities  and  have  been  very  successful  in  developing  new 
instruments. 

COLLABORATION  WITH  IMIST 

Some  time  on  all  instruments  is  available  to  NIST  staff  in  support 
of  our  mission.  This  time  is  used  to  work  on  NIST  research  needs, 
instrument  development,  and  promoting  the  widespread  use  of 
neutron  measurements  in  important  research  areas,  particularly  by  new 
users.  As  a  result  of  these  objectives,  a  significant  fraction  of  the  time 
available  to  NIST  staff  is  used  collaboratively  by  external  users,  who 
often  take  the  lead  in  the  research.  Access  through  such  collaborations 
is  managed  through  written  beam  time  requests.  In  contrast  to 
proposals,  beam  time  requests  are  reviewed  and  approved  internally  by 
NCNR  staff.  We  encourage  users  interested  in  exploring  collaborative 
research  opportunities  to  contact  an  appropriate  NCNR  staff  member. 

RESEARCH  PARTICIPATION  AT  THE  NCNR 

The  NCNR  continued  its  strong  record  of  serving  the  U.S.  research 
community  this  year.   Over  the  2008  fiscal  year,  2245  research 
participants  benefited  from  use  of  the  NCNR.  (Research  participants 
include  users  who  come  to  the  NCNR  to  use  the  facility  as  well  as 
active  collaborators,  including  co-proposers  of  approved  experiments, 
and  co-authors  of  publications  resulting  from  work  performed  at  the 
NCNR.)  Research  participants  represented  222  external  institutions, 
including  144  U.S.  universities,  38  national  laboratories,  and  40  U.S. 
industries.  The  NCNR  is  a  national  facility,  with  participants  from  40 
states,  DC,  and  Puerto  Rico. 
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In  2008  Research  Participants  were  from: 
18  NtST  Divisions  and  Offices 
38  U.S.  Government  Laboratories 
144  U.S.  Universities 
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FIGURE  1;  Research  participants  at  the  NCNR  1986-  2008 
2008  NCNR  PROPOSAL  PROGRAM 

Two  calls  in  the  past  year  resulted  in  the  submission  of  652  proposals, 
a  record  number,  of  which  414  were  approved  and  received  beam  time. 
The  oversubscription,  i.e.,  the  ratio  of  days  requested  on  all  proposals  to 
the  days  available,  was  2.2  on  the  average,  but  as  high  as  6.8  for  specific 
instruments  in  one  call.  Proposal  demand  has  grown  constantly  since 
the  NCNR  first  began  accepting  proposals  in  1991,  and  has  doubled 
in  the  past  six  years.  The  following  table  shows  the  data  for  several 
instrument  classes. 
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Instrument  class 

Proposals 

□ays  requested 

□ays  allocated 

SANS  and  USANS 

232 

883 

451 

94 

750 

404 

270 

2092 

781 

Diffraction 

23 

91 

70 

Imaging 

23 

137 

70 

Total 

652 

3853 

1776 

USER  GROUP  HOLDS  SECOND  ELECTION 

The  NCNR  Users  Group  (NUG)  provides  an  independent  forum 
for  all  facility  users  to  raise  issues  to  NCNR  management,  working 
through  its  executive  officers  to  carry  out  this  function.  In  early  2007, 
it  conducted  a  comprehensive  online  survey  of  user  satisfaction  and 
concerns,  eliciting  260  responses  with  many  useful  comments  on  areas 
such  as  instrument  quaUty  and  problems,  sample  environment,  data 
analysis  software,  and  user  laboratories.  The  survey  input  is  being  used 
as  a  guide  to  make  some  significant  improvements  in  user  support. 
In  December  2007,  the  NUG  conducted  its  second  online  election 
of  officers,  so  that  now  all  of  the  Executive  Committee  members 
have  been  elected  by  the  general  NUG  membership.  The  Committee 
members  are  Stephan  Rosenkranz  (Argonne  National  Laboratory, 
chair),  Mark  Dadmun  (University  of  Tennesee,  Knoxville),  Kalina 
Hristova  (Johns  Hopkins  University),  Moon  Jeong  Park  (University 
of  California,  Berkeley),  Tonya  Kuhl  (University  of  California,  Davis), 
Robert  Leheny  (Johns  Hopkins  University),  and  Lynn  Walker 
(Carnegie-MeUon  University). 

PANEL  OF  ASSESSMENT 

The  major  organizational  components  of  NIST  are  evaluated 
periodically  for  quality  and  effectiveness  by  the  National  Research 
Council  (NRC),  the  principal  operating  agency  of  both  the  National 
Academy  of  Sciences  and  the  National  Academy  of  Engineering. 
A  panel  appointed  by  the  NRC  reported  on  the  NIST  Center  of 
Neutron  Research  in  early  2007  and  again  in  early  2008.  Their  findings 
from  the  2007  assessment  are  summarized  in  a  document  that  may 
be  viewed  online  at  http://books.nap.edu/openbook.phpPrecord_ 
id=120108cpage=Rl 

The  panel  members  included  Philip  Pincus  of  the  University  of 
California,  Santa  Barbara  (chair),  Stuart  Bush,  Rohm  and  Haas 
Corporation,  C.  W.  Chu,  University  of  Houston,  Sebastian  Doniach, 
Stanford  University,  Christopher  Gould,  North  Carolina  State 
University,  Tonya  Kuhl,  University  of  California,  Davis,  Herbert 
Mook,  Oak  Ridge  National  Laboratory,  V.  Adrian  Parsegian,  National 
Institute  of  Child  Health  and  Human  Development, ,  Kenneth  C. 
Rogers  (retired),  Nuclear  Regulatory  Commission,  and  Barbara 
Wyslouzil,  the  Ohio  State  University. 

THE  CENTER  FOR  HIGH  RESOLUTION  NEUTRON  SCATTERING 
(CHRNS) 

The  CHRNS  is  a  national  user  facility  that  is  jointly  funded  by  the 
National  Science  Foundation  and  the  NCNR.  Its  primary  goal  is  to 
maximize  access  to  state-of-the-art  neutron  scattering  instrumentation 
for  the  research  community.  It  operates  six  neutron  spectrometers  and 
diffractometers  at  the  NCNR,  enabling  users  from  around  the  nation 
and  the  world  to  observe  dynamical  phenomena  involving  energies 
from  «  30  neV  to  ~  100  meV,  and  to  obtain  structural  information 
on  length  scales  from  1  nm  to     10  |im.  A  more  detailed  account  of 
CHRNS  activities  may  be  found  on  p.  59  of  this  report. 


PARTNERSHIPS  FOR  SPECIFIC  INSTRUMENTS 
NG-7  SANS  Consortium 

A  consortium  that  includes  NIST,  the  Industrial  Partnership  for 
Research  in  Interfacial  and  Materials  Engineering  (IPRIME)  led 
by  the  University  of  Minnesota,  and  the  ExxonMobil  Research  and 
Engineering  Company,  operates,  maintains,  and  conducts  research 
at  the  NG-7  30  m  SANS  instrument.  The  consortium  uses  65  % 
of  the  beam  time  on  this  instrument,  with  the  remaining  35  % 
allocated  to  the  general  scientific  community  through  the  NCNR's 
proposal  system.  Consortium  members  conduct  independent  research 
programs  primarily  in  the  area  of  large-scale  structure  in  soft  matter. 
For  example,  ExxonMobil  has  used  this  instrument  to  deepen  their 
understanding  of  the  underlying  nature  of  ExxonMobil's  products 
and  processes,  especially  in  the  fields  of  polymers,  complex  fluids,  and 
petroleum  mixtures. 

NIST  /  General  Motors  -  Neutron  Imaging 
An  ongoing  partnership  between  General  Motors  and  NIST  has 
resulted  in  an  exciting  collaboration  employing  neutron  imaging  to 
visualize  the  operation  of  fuel  cells  for  automotive  vehicle  applications. 
Neutron  imaging  is  an  ideal  method  for  observing  the  movement  and 
behavior  of  water  in  proton  exchange  membrane  fuel  cells,  a  key  to 
developing  hydrogen-powered  cars  that  employ  such  cells  instead  of 
internal  combustion  engines  as  their  central  power  source.  In  the  past 
year,  the  spatial  resolution  of  the  technique  has  been  improved  ten-fold. 
Experiments  were  performed  to  simulate  the  operation  of  stacks  of 
cells,  a  necessary  advance  towards  development  of  practical  fuel  cell 
vehicles. 

INTERAGENCY  COLLABORATIONS 

The  Smithsonian  Institution's  Nuclear  Laboratory  for  Archeological 
Research  is  part  of  the  Anthropology  Department  at  the  National 
Museum  of  Natural  History.  It  has  had  a  productive  31  y  partnership 
with  the  NCNR,  during  which  time  it  has  chemically  analyzed  over 
37,000  archaeological  artifacts  by  Instrumental  Neutron  Activation 
Analysis  (INAA),  drawing  extensively  on  the  collections  of  the 
Smithsonian,  as  well  as  on  those  of  many  other  institutions  in  this 
country  and  abroad.  Such  chemical  analyses  provide  a  means  of 
linking  these  diverse  collections  together  in  order  to  study  continuity 
and  change  involved  in  the  production  of  ceramic  and  other  artifacts. 

The  Center  for  Food  Safety  and  Applied  Nutrition,  U.S.  Food  and 
Drug  Administration  (FDA),  directs  and  maintains  a  facility  at  the 
NCNR  that  provides  agency-wide  analytical  support  for  food  safety 
and  food  defense  programs.  Neutron  activation  (instrumental,  neutron- 
capture  prompt-gamma,  and  radiochemical),  x-ray  fluorescence 
spectrometry,  and  low-level  gamma-ray  detection  techniques  provide 
diverse  multi-element  and  radiological  information  about  foods 
and  related  materials.  Ongoing  work  includes  studies  of  cryogenic 
homogenization  of  foods,  preparation  of  an  in-house  reference 
material,  improved  detection  capabilities  for  toxic  and  nutritional 
elements  in  foods,  rapid  screening  of  food  products  for  dangerous  levels 
of  cadmium,  lead,  and  mercury,  and,  most  recently,  a  study  of  residual 
bromine  in  baked  goods  and  beverages. 
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CHRNS  is  a  national  user  facility  that  is  jointly  fUnded  by 
the  National  Science  Foundation,  through  its  Division  of 
Materials  Research  (grant  number  DMS-0454672),  and  by  the 
NCNR.  The  primary  purpose  of  this  partnership  is  to  maximize 
access  to  state-of-the-art  neutron  scattering  instrumentation  for 
the  research  community  through  the  NCNR's  proposal  system. 
Proposals  to  use  the  CHRNS  facilities  are  critically  reviewed  on  the 
basis  of  scientific  merit  and/or  technological  importance. 

The  core  mission  of  CHRNS  is  fourfold:  (i)  to  develop  and  operate 
neutron  scattering  instrumentation,  with  broad  application  in 
materials  research,  for  use  by  the  general  scientific  community,  (ii) 
to  promote  the  effective  use  of  the  CHRNS  instruments  by  having 
an  identifiable  staff  whose  primary  function  is  to  assist  users;  (iii) 
to  conduct  research  that  advances  the  capabilities  and  utilization 
of  CHRNS  facilities;  and  (iv)  to  contribute  to  the  development  of 
human  resources  through  educational  and  outreach  efforts. 

Scattering  Instruments  and  Research 

CHRNS  currently  operates: 

1.  the  30  m  Small  Angle  Neutron  Scattering  (SANS)  instrument 
at  NG-3, 

2.  the  Ultra-SmaU  Angle  Neutron  Scattering  (USANS) 
instrument  at  BT-5, 

3.  the  Spin-Polarized  Inelastic  Neutron  Scattering  (SPINS) 
spectrometer  at  NG-5, 

4.  the  Disk  Chopper  Spectrometer  (DCS)  at  NG-4, 

5.  the  High  Flux  Backscattering  Spectrometer  (HFBS)  at  NG-2, 
and 

6.  the  Neutron  Spin-Echo  (NSE)  spectrometer  at  NG-5. 

The  small  angle  scattering  instruments  supported  by  CHRNS 
provide  structural  information  over  length  scales  from  =  1  nm  to 
=  10  \im.  The  spectrometers  collectively  yield  dynamical  information 
over  time  scales  from  =  3x10""  s  to  =  10"^  s  (energy  scales  from 


=  100  meV  to  =  30  neV).  These  wide  ranges  of  accessible  distances  and 
times  support  a  very  diverse  scientific  program,  allowing  researchers 
in  materials  science,  chemistry,  biology,  and  condensed  matter  physics 
to  investigate  materials  such  as  polymers,  metals,  ceramics,  magnetic 
materials,  porous  media,  fluids  and  gels,  and  biological  molecules. 

In  the  most  recent  Call  for  Proposals,  208  were  received  for  the  six 
CHRNS-operated  instruments,  and  131  were  awarded  beam  time. 
Of  the  1158  days  requested,  the  Beam  Time  Allocation  Committee 
approved  530.  Roughly  half  of  the  users  of  neutron-scattering  techniques 
at  the  NCNR  use  CHRNS-fijnded  instruments,  and  about  one  third  of 
NCNR  publications  (see  the  "Pubhcations"  section  on  p.  62)  over  the 
current  one-year  period  are  based  on  research  performed  using  these 
instruments.  In  2007  more  than  24  Master's  and  Ph.D.  theses  were 
completed  using  results  from  CHRNS-sponsored  instruments. 

Scientific  Support  Services 

An  important  ingredient  of  the  CHRNS  operation  is  its  support  of 
staff  whose  responsibility  is  to  provide  services  that  are  essential  to  the 
success  of  neutron  scattering  investigations.  One  such  service  is  the 
provision  and  operation  of  an  expanding  range  of  sample  environment 
equipment.  The  sample  environment  team  often  takes  on  projects  to 
expand  the  range  of  existing  equipment,  enhance  its  safety,  improve  its 
reliability,  and  even  develop  entirely  new  capabilities.  Examples  from  this 
past  year  include  a  new  SANS  high  pressure  cell  with  windows  (shown 
below),  higher  reliability  SANS  windows  for  the  9  T  horizontal  field 
magnet,  and  a  prototype  high  temperature  closed-cycle  refrigerator  with 
considerably  improved  performance. 

Other  services  include  the  maintenance  of  several  well-equipped  user 
laboratories,  and  the  development,  documentation  and  maintenance  of 
data  reduction,  visualization,  and  analysis  software  through  the  DAVE 
project.  DAVE  is  a  suite  of  programs  that  enables  users  to  rapidly  reduce, 
examine  and  display  their  experimental  data.  It  may  be  downloaded  free 
of  charge  from  the  website  http://www.ncnr.nist.gov/dave/,  either  as  a 
binary  executable  with  an  embedded  runtime  license  or  as  the  complete 


source  code.  Hie  DAVE  project  has  lately  been  expanded  to  include 
support  for  spectrometers  at  the  Paul  Scherrer  Institute  (Switzerland). 
A  major  redesign,  known  as  DAVE  2,  has  been  released  with  the  latest 
development  version. 

Education  and  Outreach 

One  of  the  missions  of  the  CHRNS  program  is  to  contribute  to  the 
development  of  human  resources  through  education  and  outreach, 
thereby  fostering  future  generations  of  scientists  and  neutron  scatterers. 

This  year,  as  in  previous  years,  CHRNS  and  the  NCNR  jointly 
sponsored  a  neutron  scattering  summer  school.  The  2008  school, 
focused  on  small  angle  scattering  and  reflectometry,  attracted  a  record 
number  of  applicants  (125)  for  the  36  available  places.  The  students' 
backgrounds  ran  the  gamut,  from  chemistry,  engineering,  physics,  and 
materials  science,  to  polymer  science,  complex  fluids,  and  biology.  The 
NCNR/CHRNS  summer  schools  have  developed  into  a  very  eifective 
way  of  educating  potential  new  users  of  the  facihties  at  the  NCNR.  For 
the  first  time  we  had  a  poster  session  where  students  presented  their 
research,  or  that  of  their  group.  Course  materials  have  been  placed  on 
the  Web,  at  http://www.ncnr.nist.gov/summerschool/ss08/. 


l|     FIGURE  2:  2008  SURF  students  accompanied  by  their  advisors 


CHRNS  also  participated  in  NIST's  Summer  Undergraduate  Research 
Fellowship  (SURF)  program,  hosting  nine  undergraduate  students 
who  worked  directly  with  staff  members  studying  materials  that  ranged 
from  magnetic  nanoparticles  to  biological  membranes.  The  students' 
work  encompassed  sample  preparation,  neutron  scattering  using  diverse 
instrumentation,  and  development  of  data  analysis  software.  The 
students  enthusiastically  agreed  that  their  summer  research  experiences 
were  productive,  expanding  their  understanding  of  the  scientific  process. 

2008  was  the  second  year  that  NIST  organized  a  "Summer  Institute" 
for  middle  school  science  teachers.  Sixteen  teachers  from  local  schools 
took  part.  At  the  NCNR  they  attended  lectures  and  demonstrations  of 
projects  suitable  for  their  students  and  toured  the  faciUty.  They  were  also 
given  the  opportunity  to  interact  with  users,  SURF  students  and  staff. 
The  program  should  give  teachers  a  greater  appreciation  of  what  can  be 
learned  using  neutron  scattering  and  how  that  knowledge  affects  the 
j     world  we  Uve  in. 


FIGURE  3:  Middle  school  teachers  and  NCIMR/CHRIMS  staff 


FIGURE  4:  DCS  tutorial  participants  and  CHRNS  staff  members 


The  increasing  number  of  summer  school  applicants  in  recent  years 
prompted  us  to  offer  a  series  of  tutorial  demonstrations  based  on  the 
experiments  and  lectures  presented  at  the  2007  Summer  School  on 
Methods  and  Applications  of  Neutron  Spectroscopy.  Twelve  participants 
from  local  universities  and  staff  members  from  other  NIST  divisions 
attended  the  tutorials. 

CHRNS  staff  participated  in  the  Joint  Annual  Conference  of  the 
National  Society  of  Black  Physicists  and  the  National  Society  of 
Hispanic  Physicists  (NSHP),  in  W^ashington,  DC.  The  conference 
gave  us  an  opportunity  to  interact  with  professors  and  students  from 
Historically  Black  Colleges  and  Universities  (HBCUs)  and  wdth 
members  of  the  NSHP.  Graduate  and  undergraduate  students  toured 
the  NCNR  facility. 

Other  education  and  outreach  activities  included  "Bring  your  Kids  to 
Work  Day",  "Adventures  in  Science",  the  Boy  Scouts  Nuclear  Science 
merit  badge,  and  tours  for  middle  school  and  high  school  students.  We 
also  had  our  first  high  school  intern. 
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Taner  Yildirim  of  the  NCNR  was  a 
recipient  of  a  (Washington,  D.C.) 
Downtown  Jaycees  2007  7\rthur  S. 
Flemming  Award  in  AppUed  Science 
for  his  "innovative  research  on  hydrogen 
storage  materials.  "The  Flemming  Awards 
honor  outstanding  federal  employees 
with  three  to  15  years  of  puWic  service  experience.  Yildirim 
joins  NCNR  Flemming  alumni  David  Jacobson  (2006)  and 
Muhammad  Arif  (2002)  in  this  honor. 


7\lan  K.  Thompson 
and  Muhammad  Arif 

of  the  Ionizing 
Radiation  Division, 
Physics  Laboratory 
and  working  at  the 
NCNR,  along  with 
Robert  E.  Vest  and 
Charles  W.  Clark  of  the  NIST  Electron  and  Optical  Physics 
Division;  and  Michael  A.  Coplan  of  the  Institute  for  Physical 
Science  and  Technology,  University  of  Maryland,  received  an 
R&D  100  Award  for  their  development  of  a  new  neutron 
detector  based  on  light  emissioji  from  the  excited  hydrogen  atoms 
that  are  a  product  of  the  neutrori  absorption  by  ^He.  R&D  100 
awards  are  given  for  "the  100  most  technologically  significant 
products  introduced  into  the  market"  during  the  previous 
year,  as  selected  by  an  independent  judging  panel  and  the 
editors  of  R&D  Magazine.  See  the  article  on  p.  49. 

Taner  YUdirim  of  the  NCNR 
was  awarded  the  Department 
of  Commerce  (DOC)  Gold 
Medal  for  "For  scientific 
breakthroughs  in  the  creation  and 
characterization  of  nanostructured 
materials  for  safe  and  efficient 


solid-state  hydrogen  storage. 


Robert  E.  Williams,  Paul  A. 
Kopetka,  and  Scott  J.  Slifer  of 

the  NCNR  were  awarded  the 
DOC  Silver  Medal  for  the 

a?zalysis,  design,  and  construction 
of  the  Liquid  Hydrogen  Cold 
Neutron  Source  at  the  NIST 


Peter  Gehring  of  the  NCNR  was  awarded 
the  DOC  Bronze  Medal  ''For  insightful 
neutron  scattering  studies  that  have  elucidated 
the  nanoscale  structure  and  dynamics  of  relaxor 
ferroelectrics. "  See  the  highlight  of  his  recent 
work  on  p.  12. 

Terrence  Udovic  ot  the  NCNR  was 
awarded  the  DOC  Bronze  Medal  "For 
leadership  of  neutron  scattering  studies  of  metal 
hydride  systems  for  practical  hydrogen  storage" 


JefF Zeigler  of  the  NCNR  was  awarded  the 
DOC  Bronze  Medal  "/vsr  innovative  design 
and  implejnentation  of  electronic  systems  for 
neutron  instrumentation''' 


Robert  Greenberg  of  the  NIST  Analytical 
Chemistry  Division  was  awarded  the  DOC 
Bronze  Medal  "For  leadership  in  the 
develop?7tent  and  application  of  neutron 
activation  ajialysis  as  a  primary  method  of 
chemical  analysis" 

Jason  Gardner  of  the  NCNR  and  Indiana 
University  has  been  appointed  Fellow  in  the 
Institute  of  Physics,  UK:  "For  his  outstanding 
contributioJis  to  our  understanding  of  the  spin 
dynamics  and  magnetic  order  in  complex 
magnetic  systems  by  means  of  neutron  scattering 
a?2d  contributions  to  the  lOP journals." 

Yun  Liu  of  the  NCNR  was  selected  by  the 
NIST  chapter  of  Sigma  Xi  for  a  2008 
outstanding  poster  award  for  his 

presentation:  "Increasing  hydrogen  uptake  by 
increasing  the  density  of  adsorbed  hydrogen. " 
The  work  is  highlighted  on  p.  30. 


Center  for  Neutron  Research. ' 


60 


Jae-Hyuk  Her,  SUNY  Stony  Brook,  and 
currently  at  the  NCNR  and  the  U.  of 
Maryland,  was  awarded  the  Margaret  C. 
Etter  Student  Lecturer  Award  by  the 

American  Crystallographic  Association  at  its 
July  2008  meeting  for  his  presentation: 

"'Maximum  Entropy  Crystal  Structure 
Reconstruction  using  X-ray  Powder  Diffraction  Data.'' 

Prof.  Sow-Hsin  Chen  of  MIT  is  the 
recipient  of  the  2008  Clifford  G.  ShuU  Prize 

of  the  Neutron  Scattering  Society  of  America 
with  the  citation:  "For  seminal  contributions  to 
understanding  the  dynamical  properties  of 
supercooled  and  interfacial  water  using  neutron 
scattering  techniques,  and for  an  exceptional 
record  of  training  young  scientists  in  the  use  of scattering  techniques  to 
solve  topical  interdisciplinary  problems  in  complex fluids  and  soft 
matter."  Many  of  Prof  Chen's  recent  measurements  were  carried 
out  at  the  NCNR. 

Prof.  Frank  Bates  (U.  Minnesota)  is  the 
recipient  of  the  2008  Sustained  Research 

Prize  of  the  Neutron  Scattering  Society  of 
America  with  the  citation:  "For  his  pioneering 
SANS  experiments  that  probe  the  structure  and 
thermodynamics  of  polymeric fluids  and  block 
copolymers. "  The  U.  of  M.  shares  in  a  PRT 
operating  the  NG-7  SANS  instrument  at  the  NCNR  where 
many  of  Prof.  Bates  experiments  were  carried  out. 

■■I^^^V'  Prof.  Seung-Hun  Lee  (U.  Va)  is  the  recipient 
V^^BP^^mH  2008  Science  Prize  of  the  Neutron 

^H^^P  Scattering  Society  of  America  with  the 

^^^^W  citation:  "For  his  innovative  and  insightful 

^^^^^^^^^^  neutron  scattering  studies  of  frustrated  magnetic 

^^^^^  ^^^1  systems. "  Much  of  this  work  was  done  when 

^^^Ui//tWKKM  Prof  Lee  spent  the  first  part  of  his  career  at 
the  NCNR  where  he  helped  develop  the  SPINS  instrument. 


Thomas  Trabold,  Jon  Owejan,  and  Jeffrey  Gagliardo  from 
General  Motors  Fuel  Cell  Activities  (Honeoye  Falls,  NY)  were 
recently  honored  with  GM's  Charles  L.  McCuen  Special 
Achievement  Award  in  recognition  oi" Extraordinary 
Accomplishment  in  the  area  of  Neutron  Imaging  Methods  for  Fuel  Cell 
Water  Visualization." ]tS Is  stationed  at  the  NCNR  where  the 
work  was  performed  at  the  Neutron  Imaging  Facility. 

The  Faculty  of  Science  of  Stockholm 
University  has  honored  Prof.  Stephen  M. 
White,  U.  California,  Irvine  with  a  conferral 
of  a  Doctorate  of  Philosophy  honoris  causa 

in  recognition  of  his  work  on  biological 
membranes.  Prof.  White's  team  has  made 
significant  advances  in  membrane 
measurements  using  neutron  reflectometry  at 
NCNR's  AND/R  instrument. 

Kevin  Yager,  NIST  Polymers  Division,  was 
selected  by  the  NIST  chapter  of  Sigma  Xi  for 
a  2008  outstanding  poster  award  for  his 

presentation:  'Neutron  Scattering 
Measurements  of  Nanoparticle  Directed 
Self-Assemply  of  Block-Copolymers.  "Yager's 
work  is  highlighted  on  p.  40. 

Matthew  E.  Helgeson,  Dept.  of  Chemical 
Engineering,  U.  of  Delaware,  a  graduate 
student  working  with  Prof  Norman  Wagner, 
has  been  awarded  the  NSSA  Prize  for 
Outstanding  Student  Research  for  his 
project  entitled:  ''Spatially-resolved  Structural 
Evolution  of  Wormlike  Micelles  through  the  Shear 
Banding  Transition."The  work  is  highlighted  in  this  report  on  p.  42. 
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(2007). 
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NCNR  2008  Accomplishments  and  Opportunities 


Instruments  and  Contacts 


Instruments  and  Contacts:  (name,  tel.  301-975-xxxx,  email) 

High  Resolution  Powder  Diffractometer  (BT-1): 

•  Judy  Stalick,  6223,  judith.stalick@nist.gov 

•  Qingzhen  Huang,  6164,  qing.huang@nist.gov 

•  Mark  Green,  4297,  mark.green@nist.gov 

•  Sarah  Pouhon,  8491,  sarah.poulton@nist.gov 

Residual  Stress  Diffractometer  (BT-8): 

•  Thomas  Gnaeupel-Herold,  5380,  tg-h@nist.gov 

30-m  SANS  Instrument  (NG-7): 

•  Mark  Laver,  4404,  mark.laver@nist.gov 

•  Paul  Butler,  2028,  paul.butler@nist.gov 

•  JefFKrzywon,  6650,  jkrzywfon@nist.gov 

30-m  SANS  instrument  (NG-3)  (CHRNS): 

•  Boualem  Hammouda,  3961,  hammouda@nist.gov 

•  Steve  Kline,  6243,  steven.kline@nist.gov 

•  Susan  Krueger,  6734,  susan.krueger@nist.gov 

•  Cedric  Gagnon,  2020,  cedric.gagnon@nist.gov 

USANS,  Perfect  Crystal  SANS  (BT-5)  (CHRNS): 

•  Andrew  Jackson,  3146,  andrew.jackson@nist.gov 

•  David  Mildner,  6366,  david.mildner@nist.gov 

•  Paul  Butler,  2028,  paul.butler@nist.gov 

•  John  Barker,  6732,  barker@nist.gov 

Reflectometer,  Vertical  Sample,  Polarized  Beam  Option 
(NG-1): 

•  Brian  Kirby,  8395,  brian.kirby@nist.gov 

•  Julie  Borchers,  6597,  julie.borchers@nist.gov 

•  Chuck  Majkrzak,  5251,  cmajkrzak@nist.gov 

AND/R,  Advanced  Neutron  Diffractometer/Reflectometer 

(NG-1): 

•  Joe  Dura,  6251,  joseph.dura@nist.gov 

•  Brian  Maranville,  6034,  brian.maranville@nist.gov 

Reflectometer,  Horizontal  Sample  (NG-7): 

•  Sushil  Satija,  5250,  satija@nist.gov 

•  Bulent  Akgun,  6469,  bulent.akgun@nist.gov 

Double-focusing  Triple  Axis  Spectrometer  (BT-7): 

•  Ying  Chen,  6442,  ying.chen@nist.gov 

•  Jeff  Lynn,  6246,  jeff.lynn@nist.gov 

Triple-axis  Spectrometer  (BT-9): 

•  William  Ratcliff,  4316,  wiUiam.ratcliff@nist.gov 

•  Ross  Erwin,  6245,  rerwin@nist.gov 

•  Jeff  Lynn,  6246,  jeff.lynn@nist.gov 

SPINS,  Spin-polarized  Triple-axis  Spectrometer  (NG-5) 
(CHRNS): 

•  Sung  Chang,  8369,  sung.chang@nist.gov 


FANS,  Filter-analyzer  Neutron  Spectrometer  (BT-4): 

•  Terry  Udovic,  6241,  udovic@nist.gov 

•  Juscelino  Leao,  8867,  juscelino.leao@nist.gov 

DCS,  Disk-chopper  Time-of-flight  Spectrometer  (NG-4) 
(CHRNS): 

•  John  Copley,  5133,  jcopley@nist.gov 

•  Yiming  Qiu,  3274,  yiming.qiu@nist.gov 

•  Craig  Brown,  5134,  craig.brown@nist.gov 

HFBS,  High-flux  Backscattering  Spectrometer  (NG-2)  (CHRNS): 

•  Madhusudan  Tyagi,  2046,  mtyagi@nist.gov 

•  Tim  Jenkins,  8396,  timothy.jenkins@nist.gov 

NSE,  Neutron  Spin  Echo  Spectrometer  (NG-5)  (CHRNS): 

•  Antonio  Faraone,  5254,  antonio.faraone@nist.gov 

•  Jason  Gardner,  8396,  jason.gardner@nist.gov 

•  Michi  Nagao,  5505,  michihiro.nagao@nist.gov 

Cold  Prompt-gamma  Neutron  Activation  Analysis  (NG-7): 

•  Rick  Paul,  6287,  rpaul@nist.gov 

Thermal  Prompt-gamma  Neutron  Activation  Analysis  (VT-5): 

•  Liz  Mackey,  5149,  liz.mackey@nist.gov 

Other  Activation  Analysis  Facilities: 

•  Greg  Downing,  3782,  gregory.downing@nist.gov 

Cold  Neutron  Depth  Profiling  (NG-1): 

•  Greg  Downing,  3782,  gregory.downing@nist.gov 

Neutron  ivriag<ng  Station  {BT-2'): 

•  David  Jacobson,  6207,  david.jacobson@nist.gov 

•  Dan  Hussey,  6465,  daniel.hussey@nist.gov 

•  Muhammad  Arif,  6303,  muhammad.arif@nist.gov 

Neutron  Interferometer  (NG-7): 

•  Muhammad  Arif,  6303,  muhammad.arif@nist.gov 

•  Dave  Jacobson,  6207,  david.jacobson@nist.gov 

•  Dan  Hussey,  6465,  daniel.hussey@nist.gov 

Fundamental  Neutron  Physics  Station  (NG-6): 

•  NG-6M:  Scott  Dewey,  4843,  mdewey@nist.gov 

•  NG-6U:  Pieter  Mumm,  8355,  pieter.mumm@nist.gov 

•  NG-6:  Jeff  Nico,  4663,  nico@nist.gov 

Theory  and  Modeling: 

•  Taner  Yildirim,  6228,  taner@nist.gov 

•  Joseph  Curtis,  3959,  joseph.curtis@nist.gov 

•  Norm  Berk,  6224,  nft)@nist.gov 

Sample  Environment: 

•  Julia  Scherschligt,  5328,  julia.scherschligt@nist.gov 

•  Juscelino  Leao,  8867,  juscelino.leao@nist.gov 

•  Dan  Dender,  6225,  dender@nist.gov 
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